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We have proposed that the material spalled 
from the surface of an exploding supernova is 
the source of cosmic radiation. Nuclei in the sur- 
face of the star may acquire many orders of mag- 
nitude more than the average energy per particle 
released in the explosion because of the large 
ratio of matter density between the core and the 
outer mantle. A shock from a sudden pressure 
increase in the core intensifies as it advances 
into lower-density material, thereby imparting 
extreme relativistic energies to the outermost 
layers. We are presenting a preliminary account 
of our investigations covering the mechanism of 
supernova explosions, the hydrodynamic motion 
following the explosion, the subsequent flow of 
radiation, and the structure of relativistic shock 
fronts, all of which are significant influences on 
cosmic radiation. 

Astronomical evidence indicates that a super- 
nova is a highly evolved star with a central tem- 
perature of 300-500 kev in which the nuclear fuel 
has burned to the most stable form, i.e., Fe. If 
we assume the total mass of such a star to be 
10M, (Mg =mass of the sun), the radial density 
and temperature distribution as determined from 
hydrostatic and radiative equilibrium conditions 
are those shown in Fig. 1. In the explosion mech- 
anism currently under study,* the energy source 
is the gravitational instability occurring when the 
adiabatic equation of state in the core corresponds 


TEMPERATURE IN ev 



















lev 10 io? 10° 10% 10° ~~ 108 
1 1 
io’ Fe a“ 
10 F 7 
10> 4 
oO 
o 
~ 
oo 
§ ot F 
~ DENSITY 
(vs RADIUS) 
2 
= ot 4 
> 
3 
” eh. 
2 10 5 
Ww 
Oo 
TEMPERATURE 
iO F (vs DENSITY ) “7 
ae ail 
10'e - 
RADIUS x 108 cm 
wt ¢ 3 2 > * 7 F F 2 2 4 
FIG. 1. Pre-supernova density and temperature 
distributions. The central density for 10 solar masses 


has been chosen from the estimates of Burbridge 


et al., ' with a density distribution corresponding to 


radiative equilibrium (p = T3- 4%). 
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to a ratio of specific heats of less than 4/3. This 
condition can arise for temperatures in the neigh- 
borhood of 500 kev, at which nuclei begin to dis- 
sociate and the resulting dissociation energy is a 
considerable part of the heat capacity. We hope 
ultimately to start the hydrodynamic calculation 
with the onset of instability? and follow in detail 
the gravitational collapse to provide the initial 
conditions for the subsequent explosion. 

Current hydrodynamic calculations are made 
with a one-dimensional hydrodynamic code.* The 
validity of the code in following shock propagation 
through a medium in which the density changes by 
many orders of magnitude has been verified by 
comparison with known analytic solutions‘ for 
special cases. 

Since astronomical evidence indicates that ap- 
proximately one-tenth of the mass (i.e., 1M,) is 
ejected in a supernova explosion, the calculation 
begins with an energy release in the core adjusted 
so that one solar mass escapes from the gravi- 
tational field. The required energy is about 5x10" 
ergs/g or 10°" ergs altogether which gives the 
necessary 2 10"* ergs/g to the outer solar mass. 
The resulting fluid velocity behind the shock is 
shown in Fig. 2. The velocity increases by ap- 
proximately a factor 10 to a value of 1x10"° cm/sec 
at the point where the density has fallen to 30 
g/cm*. The fluid velocity well after the passage 
of the shock, also shown in Fig. 2, has reached 
twice the corresponding fluid velocity behind the 
shock. 

At velocities greater than 10° cm/sec, rela- 


tivistic effects become important. Relativistic 
equations of motion (including shock continuity 
conditions), obtained by setting the divergence of 
the energy-momentum tensor equal to zero, have | 
been considerably studied.* In the present appli- 
cation, the internal energy behind the shock is 
almost entirely in the radiation field; consequently 
the temperature never exceeds 100 kev and the 
equation of state is that of a radiation gas. The | 
relativistic equations are being programmed for / 
the IBM 709 computer, so that the numerical cal- 
culations shown in Fig. 2 can be extended to the 
lower densities at which the energies of cosmic- 
ray magnitude appear. 

Meanwhile, we have found analytic solutions of 
considerable generality. Characteristic space- 
time curves exist for the one-dimensional rela- 
tivistic equations. They are the trajectories of 
sound waves traveling in the same or the opposite 
direction to the fluid. They differ from nonrela- 
tivistic characteristics only in that the fluid and 
sound velocities are added according to the rela- | 
tivistic composition law for velocities. Incase | 
the fluid is isentropic or in case the fluid satis- | 
fies the equation of state p=E /3, the character- 
istic equations are integrable. For the latter 
case, the quantities [(1+)/(1- 8) ]E¥*” and 
[(1+)/(1 - 8) JE-¥*” are constant on the forward 
and backward characteristic, respectively, cor- 
responding to the Riemann invariants of classical 
hydrodynamics. In these relations £ is the fluid 
velocity divided by the velocity of light c, p is 
the pressure in the fluid at rest, and E is the 
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FIG. 2. Fluid velocity. 
Curve I: the velocity of the 
fluid immediately behind the 
shock as a function of initial 
density at the shock position. 
Curve II: the final velocity of 
4 the fluid after expansion with 
gravitation as a function of its 
initial density. These calcu- 
lations were made for the star 
shown in Fig. 1 where the ini- 
tial energy inserted at the 
center of the star has been ad- 
= justed so that the total energy 
behind the shock corresponds 
to escape for all points beyond 
the solar mass cut. 
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proper energy density. The expansion into a 
vacuum can be treated analogously to the corre- 
sponding well-known problem of classical hydro- 
dynamics. If the initial energy density of the 
fluid is up,c*, where p, is the rest mass density 
of the constituent fluid particles, the kinetic en- 
ergy per gram of rest mass after expansion is 
4u’ for y>>1. The flow behind a shock can also 
be integrated because, in the case of interest, 
the preponderance of internal energy behind the 
shock is in the radiation field and p=E/3 accurate- 
ly represents the equation of state. Then the final 
kinetic energy per gram of material having an 
initial density p, is (3)'+ ¥3(E,/p,c?)®-5¢8-V3), 

where E,‘*? is the value of the invariant on the 
forward characteristic, assumed the same at 

the shock front and in the subsequent expansion. 
An approximate value of E, is 30c* ergs/cc, cor- 
responding to the density 30 g/cm® at which the 
shock becomes relativistic. A more accurate 
value of E, as well as the actual change in E, 
during the expansion will be obtained from the 
numerical calculations. 

The final kinetic energy of the expanded ma- 
terial gives directly the energy distribution of 
cosmic radiation. According to the polytrop 
solution® of the hydrostatic and radiative equili- 
brium conditions, the density p at any radius r 
is proportional to the temperature T at the same 
radius to the 3.25 power. In the outer layers of 
the star where the pressure is just the weight of 
material beyond r, the mass of this material is 
then proportional to pT~=p****. From the previous 
paragraph, y, the final energy of material orig- 
inally at density p, is proportional to p-*%, 
whence M(<p) <N(>y) «y~%-3¥0-6=4-2. It is ex- 
pected that the numerical calculations for more 
accurate values of E, referred to above will show 
a somewhat less rapid decrease in intensity at 
higher energies, in better agreement with the 
experimental law, N(>+y)a y “8,798 

Above 10* Bev the initial density estimated 
from the analytic solution is less than 10°° g/cm‘, 
and the mass in the remaining mantle of the star, 
less than 1 g/cm?, is insufficient to contain the 
radiation field. However, it appears that plasma 
oscillations (to be discussed below) provide a 
mechanism by which the shock can propagate to 
much lower densities. Such a shock would be 
dissipated when the Debye length becomes equal 
to the scale height which occurs at a density of 
10™ g/em*. The corresponding material energy 
is about 10° Bev per nucleon, which is about the 
upper limit of the observed cosmic-ray spectrum. 


The absolute intensity of cosmic radiation can 
be estimated from the frequency of supernova 
explosions, about 1 per 100 years, and the mean 
life of a cosmic ray within the galaxy, about 
5x10*® years, corresponding to an average den- 
sity of 0.1 particle/cm® in the galaxy.® Since the 
volume of the galaxy is = 5x10°* cm’, the radia- 
tion from a single explosion therefore occupies 
on the average a volume of 10% cm*. We estimate 
that the initial density corresponding to a final 
expansion energy of 10 Bev is about 1 g/cm’. 

The mass in the star’s mantle beyond this point 

is about 107° g or 6 x10** protons, which distrib- 
uted over a volume of 10% cm* gives a density of 
6x10°** per cm*. The corresponding total cosmic- 
ray flux is 0.02 per cm? sec, which agrees with 
experiment within the uncertainties in the esti- 
mate. As the cosmic-ray population at any in- 
stant represents the product of 510° supernova 
explosions diffused through the galaxy, density 
fluctuations should be very small. 

The distribution of elements in cosmic radiation 
agrees well with that expected in the outermost 
parts of the supernova star as represented by the 
outer mantle of the sun. However, it is probable 
that the distribution is altered in the relativistic 
shock front. Viewed from the shock frame of 
reference, material is entering the shock front 
with velocity c and, in equilibrium with the radia- 
tion field, is leaving the shock front with a veloc- 
ity c/3. Processes within the shock front must 
produce the created equilibrium radiation field 
and must also decelerate the ions to the velocity 
c/3. The radiation field acts directly on the 
electrons so that the width of the shock front is 
determined in part by the Compton scattering 
length and in part by the distance required for 
bremsstrahlung to effect radiative equilibrium. 
Within the shock width an electric field due to 
charge separation by the radiation pressures on 
the electrons decelerates the ions, predominantly 
hydrogen. Since the charge-to-mass ratio of the 
nonhydrogenic nuclei is one half that of hydrogen, 
these nuclei lose only one half their energy in 
traversing the electric field. Proton collisions, 
which are in part responsible for bringing these 
nuclei to rest, then become sufficiently disruptive 
to destroy the nuclear species for final energies 
larger than a few Bev. 

The shock structure can be profoundly modified 
by plasma oscillations. The relative velocity of 
electrons and ions due to the charge separation 
described above produces a two-stream instability 
by which plasma oscillations may grow. Conse- 
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quently, it is likely that the internal energy in the tConsultant to the Lawrence Radiation Laboratory. 
first part of the shock is entirely carried in large- ‘E, M. Burbridge, G. R. Burbridge, W. A. Fowler, 
amplitude oscillations.‘ The majority of heavy and F. Hoyle, Revs. Modern Phys. 29, 547 (1957). 
ions may then be brought to relative rest by virtue *W. H. Grasberger is currently performing these 








calculations. 
of their phase with respect to these oscillations 3R. H. White has performed these calculations using 
with very little destruction of the nuclear species. a hydrodynamic code which has been modified to include 
As mentioned above, these oscillations also pro- gravitation. 
vide a mechanism for the transmission of the ‘We are indebted to J. M. Burgers who brought to 
shock to much lower densities. The damping of our attention one such solution in plane parallel geome- 
such oscillations by wave and by particle inter- try; W. P. Robbertse and J. M. Burgers, Proc. 
actions and behavior of the material after the aan” Akad. Wetenschap. Amsterdam 52, 968, 
oscillations break up are subjects currently under 51. D. Landau and E. M. Lifshitz, Fluid Mechanics 
study. (Pergamon Press, New York, 1959). 

A final problem under investigation is the dis- 6M. Schwarzschild, Structure and Evolution of the 
position of the blown-off material with energies Stars (Princeton University Press, Princeton, New 
below a few Bev. The amount of such material Jersey, 1958). 
is at least one order of magnitude too large to be 'S. F. Singer, in Progress in Elementary Particle 





and Cosmic-Ray Physics, edited by J. G. Wiison and 
S. A. Wouthuysen (Interscience Publishers, Inc., New 
York, 1958), Vol. IV. 





contained as cosmic radiation by the estimated 
galactic magnetic fields. We are examining mech- 


anisms whereby such material is confined by *L,. Biermann has pointed out that the cosmic-ray 
local magnetic fields to the neighborhood of the spectrum is known considerably more accurately than 
explosion long enough to lose its energy by syn- the errors in the equipment might indicate. However, 
chrotron radiation, as appears to be the present our present calculations are uncertain by more than 
situation in the Crab nebula. the indicated errors; L. Biermann and L. Davis, Z. 


It is a pleasure to acknowledge discussions with Naturforsch. 13a, 909 (1958). ” . 
*The average density includes the “galactic halo’ of 


and suggestions from many of our colleagues. ; 

We are particularly indebted in this regard to trapped cosmic rays. An extensive review is given by 
V. L. Ginzburg, in Progress in Elementary Particle 

Dr. Edward Teller. and Cosmic-Ray Physics, edited by J. G. Wilson and 
S. A. Wouthuysen (Interscience Publishers, Inc., New 
York, 1958), Vol. IV. 

*Work done under the auspices of the U. S. Atomic 10©, N. Parker and D. A. Tidman, Phys. Fluids 3, 
Energy Commission. 369 (1960). 














PLASMA RESONANCE IN A RADIO-FREQUENCY PROBE 


K. Takayama and H. Ikegami 
Electrical Communication Laboratory, Musashino-shi, Tokyo, Japan 


and 


S. Miyazaki 
Radio Research Laboratory, Kokubunji, Tokyo, Japan 
(Received July 21, 1960) 


As the frequency of a low-voltage rf signal where N, is the electron density. 
superposed on a probe is swept within an adequate The schematic diagram of the apparatus used 
frequency range, the dc component’ of the elec- is shown in Fig. 1(a). The rf probe does not 
tron current due to the nonlinearity of the sheath necessarily have to be a cathode-probe system, 
impedance is measured. It has been found that a where the rf voltage is applied between the cathode 
resonant increase appears in the current at the and the probe. The experimental results are qual- 
electron plasma frequency, itatively unaffected by the shape or dimensions 
of the probe. In this experiment a plane probe 
“p =(4nNe*/ my, (1) consisting of a nickel disk 2 cm in diameter is 
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FIG. 1. (A) Characteristic curves of the rf probe 
with a constant discharge current of 5 ma, correspond- 
ing to 5.6x10° electrons per cc, for various rf voltages 
ranging from 0.2 to 0.8 volt. The electron tempera- 
ture is 1800°K. (B) Characteristic curves for various 
discharge currents corresponding to 1.2x10°, 6.6x10*, 
and 7.9x10® electrons per cc. (C) Resonance frequency 
versus electron density determined by the Langmuir 
probe method. The solid line denotes Eq. (1). 


used. It is biased relative to the cathode so as 
to be always negative with respect to the plasma 
potential, even when the rf voltage is superposed 
on the probe. 

The discharge tube” is a Pyrex glass sphere 
30 cm in diameter filled with mercury vapor at 
a pressure of 10°‘ mm Hg. Two anode-hot cathode 
systems are arranged in such a way as to pro- 
duce a very low density and homogeneous plasma, 


and the discharge currents can control the elec- 
tron density. The plasma in the tube being of a 
diffusive type, there is no beam present. 

The relation between the frequency of the rf 
voltage and the corresponding dc component of 
electron current to the probe is shown in Fig. 1(a) 
for various amplitudes of rf voltage with a con- 
stant discharge current of 5 ma, and with the 
probe set 0.5 volt below the plasma potential. In 
Fig. 1(b), however, the probe is set at the wall 
potential. The curves in this case are for various 
discharge currents with a constant rf potential 
difference of 0.4 volt. 

Every characteristic curve in Fig. 1 consists 
of the following three frequency ranges: the first 
range where the electron current maintains a 
constant value independent of frequency, the 
second range where a resonance peak appears, 
and the third range where the superposed rf field 
has no effect on the dc current to the probe. 

In the first frequency range, the electron cur- 
rent density at the probe surface for a super- 
posed rf voltage V coswt is given by 


j=j,(exp(-eV coswt/kT)) ay 
=jyJo(ieV/kT), 


where the average is taken over a long time. 
J,(iz) is the zeroth order Bessel function for a 
pure imaginary argument, and 


jo=N e(kRT/21m)” exp(-eV,/kT). 


(2) 


T is the electron temperature, and V, is the 
potential difference between the plasma and the 
probe. 

jo corresponds to the constant current in the 
third frequency range, since the Boltzmann dis- 
tribution of electron density for the oscillating 
electric field in Eq. (2) is no longer valid when 
the frequency is higher than the electron plasma 
frequency. Accordingly Eq. (2) holds when the 
field changes slowly in a time interval of the 
order of wp™’. The ratio of j to j, gives the elec- 
tron temperature through the argument of the 
Bessel function. 

In the second frequency range, it will be noted 
from Fig. 1(a) that the positions of the peaks 
along the frequency axis coincide with each other 
whenever the discharge current is the same. In 
Fig. 1(c) the values of resonance frequency are 
plotted against the electron density determined 
by the well-known Langmuir probe method. The 
solid line was calculated from Eq. (1). It is seen 
from the figure that the resonance occurs at the 
electron plasma frequency. 
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The height of a resonance peak depends not onlv 
on the magnitude of the rf voltage, but on the col- 
lision frequency of electrons with atoms. Thus 


this height decreases as the pressure is increased. 


The mechanism of the resonance is interpreted 
in terms of the interaction between the electron 
beam and the plasma medium.*~* 

In the experiment, the oscillating electric field 
applied between the plasma and the probe drives 
off some electrons from the positive ion sheath 
around the probe into the plasma during one half- 
cycle, and during the next half-cycle, because 
of the shielding effect of the plasma, the field 
cannot restore these electrons which were driven 
off. These electrons which are driven off can be 
regarded as an electron beam with its density 
modulated by the frequency of the oscillating 
field. 

An electron beam of this kind is known to ex- 
cite growing plasma waves‘ when the kinetic 
energy of the beam exceeds the thermal energy 
of the electrons of the medium. 

This interpretation, in terms of a beam-plasma 
interaction, seems to be supported by the exper- 
iments shown in Fig. 2, where the resonance 
peak nearly vanishes when the energy supplied 
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FIG. 2. The height of resonance peaks inya versus 
rf potential difference in volts, for several discharge 
currents. The electron thermal energy is 0.17 ev 
which is shown by a short line on the abscissa. 
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by the field to the electron beam becomes com- 
parable to the thermal energy of the electrons. 

The growing wave gives a density fluctuation 
N, of the following form: 


N, =A exp(yd +i¢ -iwt) 
=N exp(-iwt), (3) 


where d is a suitable distance; ¢ is a phase dif- 
ference, expected to be 7, given, in terms of 
the beam velocity u, by 


g=wd/u=n, 


and y is a growth rate which is a function of w 
having a sharp maximum at w =wp and rapidly 
decreasing to zero for w different from wp.*»* 

The dc electron current in this case is given 
by 


4=FoJo(is)[1 -iNJ,(is)/N,J(is)], (4) 


where s=eV/kT, and J, is the first order Bessel 
function. In the first frequency range, N is 
small enough that (4) reduces to (2). 

At a frequency nearly equal to the electron 
plasma frequency, y in (3) gives a maximum 
value of N, and j shows a resonant increase. 
Hence the resonance can be fully explained by 
the behavior of N. 

The authors would like to thank Professor 
S. Kojima for his advice and encouragement, 
and also Dr. M. Sumi and many members of our 
Laboratory for their interesting discussions. 
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INFRARED MAGNETOREFLECTION IN BISMUTH. I. HIGH FIELDS 


B. Lax, J. G. Mavroides, H. J. Zeiger, and R. J. Keyes 
Lincoln Laboratory, * Massachusetts Institute of Technology, Lexington, Massachusetts 
(Received August 4, 1960) 


The pulsed-field infrared experiments on bis- 

muth reported by Keyes et al.’ at wavelengths 
between 10 and 22 » and at room temperature, 
were originally interpreted as transitions cor- 
responding to cyclotron resonance. Well-defined 
dispersion curves were obtained by reflection 
techniques. The value of the magnetic field at 
the center of each line and the measured wave- 
length were utilized to obtain an apparent effec- 
tive mass from the well-known cyclotron equation 
m*=cH/w,c. Three distinct sets of mass values 
obtained from the strongest absorption data’ 
yielded apparent effective masses for electrons 
which showed a large variation with magnetic 
field, with mass values above and below those 
deduced from the de Haas-van Alphen effect? and 
microwave cyclotron resonance.* This suggested 
that the large change in the effective masses was 
due to a rapid change in the curvature of the 
energy bands. To explain this a theoretical ex- 
pression was developed for the energy-momentum 
relation, assuming only two sets of interacting 
energy bands which are located near the edge of 
the zone. The derivation uses the k+p perturba - 
tion method, following the work of Kane* on InSb 
and that of Cohen and Blount® on bismuth. We 
showed that the nonparabolic energy dependence 
of the bands could be represented by the expres- 


sion® 
€ > a- 
(1+ Z)-F G5 (1) 
& 


where €g is the energy gap and a is the inverse 
mass tensor at the bottom of the conduction band 
and (with a minus) at the top of the valence band. 
When a magnetic field is applied, the Hamilton- 
ian including spin becomes 


> a> Bo > 1 
%2[ > +50 g£ i Tree,” (2) 


where ll=p+eA/c, A is the magnetic vector 
potential, 8 is the Bohr magneton, and g is the 
effective spectroscopic splitting tensor® at the 
bottom of the bands. The solution of Eq. (2) with 
the magnetic field taken in the z direction yields 
an equation for the eigenvalues for the valence 
and conduction bands, 


€ hk? 
—_— = 1 
(te ) (n+ gw + = 





Her, (3) 





Wk 2 v2 
cya Heg|idegh for dite oe 2m. + $8 at" ; 
(4) 


w.=eH/m*c, cyclotron mass m* =[a,a,0,mz]"™, 
Mz =, ~*A,? +@,~1A,7 +571”, and Sore =[2,7A," 
+£"A," +83 A3']”; the a; are the components of 
the inverse mass tensor a and the A; are the 
directional cosines of the magnetic field relative 
to the principal ellipsoidal axes. The plus sign 
in front of the radical in Eq. (4) gives the energy 
of the conduction band and the minus sign, that 

of the valence band. 

In attempting to fit the experimental curves for 
bismuth with Eq. (4), assuming cyclotron reso- 
nance, it was necessary to choose very small 
values for the energy gap ég (~0.01 ev) and also 
for the effective masses; at room temperature 
the data yielded masses at the Fermi level ap- 
proximately one quarter of those found at low 
temperatures.”** In addition, two other dilemmas 
were present. The first involved the existence of 
well-defined “resonant” reflection traces, such 
as shown in Fig. 1. This was in contradiction 
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for bismuth obtained by reflection using pulsed magnetic 
fields at infrared wavelength, A=13.0 microns. 
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with the expected magnetoplasma effect for bis- 
muth at room temperature, where the plasma 
wavelength occurs at ~10 uw for a carrier con- 
centration of m= 10**/cm* as derived from gal- 
vanomagnetic data.” From theoretical consider - 
ations, a resonant reflection trace should not 
have been observed over at least part of the wave- 
length range (i.e., near 20 u). The second puzzle 
involved the appearance of two resonances with 
the magnetic field in the binary direction [1120], 
where only one was expected. The intensity of 
the extra line appeared greater than that of the 
line which agreed best with known values of 
masses. 

These difficulties all seem to be resolved if it 
is assumed that the observed transitions do not 
correspond to intraband cyclotron resonance but 
are interband transitions between magnetic levels 
of the valence and conduction bands near the edge 
of the Brillouin zone. Such direct transitions 
are known to have high absorption coefficients 
from similar studies in semiconductors.* The 
estimated penetration depth for this process in 
bismuth is of the order ~0.1 4, whereas that for 
free carrier penetration due to the presence of a 
plasma (for wrt > 1) would be one order of magni- 
tude larger, i.e., ~1 yu. 

The explanation of the first dilemma can be 
considered in a more general context as suggested 
by Dresselhaus.® The argument is similar to that 
of Anderson’® for the direct observation of cyclo- 
tron resonance of minority carriers in the pres- 
ence of majority carriers even at a frequency 
well below the plasma frequency. If the complex 
conductivity of the medium consists of a large 
part o(w,H) which varies slowly with magnetic 
field and a small part Ao(w, H) which varies 
“resonantly” with magnetic field, then the “reso- 
nance” line will be observable, modified in shape 
only slightly by the background absorption. Thus 
a trace as in Fig. 1 of reflection versus magnetic 
field will represent the “resonance” due to 
A4o(w,H). In the particular case of bismuth the 
dominant contributions to o(w, H) are the non- 
resonant background interband transitions and 
Ao(w, H) is due to the small “resonant” interband 
transitions between Landau levels. The analysis 
for such a situation in the absence of a plasma 
has been treated for semiconductors." The im- 
portant conclusion for other metals is that the 
presence of a large plasma contribution o(w, H) 
due to the conduction electrons should not inter - 
fere with the observation of “resonant” interband 
transitions between Landau levels, although the 
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frequency at which these observations are made 
is below the plasma frequency. However, the 
intensity of the magnetoreflection lines in a more 
highly conducting metal would probably be smaller 
than that in bismuth. 

In order to test the hypothesis that interband 
transitions have been observed, the pulse data 
reported were replotted in terms of photon energy 
as a function of magnetic field. The three sets of 
lines converged to a value at zero field of 0.040 
+0.005 ev, corresponding to the energy gap near 
the edge of the zone at room temperature. Had 
these been cyclotron resonance transitions the 
lines would have converged to the origin. Further- 
more, upon re-examining the original photographic 
traces, we found additional weak lines correspond- 
ing to what we believe to be transitions between 
higher quantum states. With this in mind, it is 
now possible to interpret the data from a model 
of the energy level diagram shown in Fig. 2, in- 
cluding spin splitting.* However, since the scat- 
ter of the present data was large, only a semi- 
quantitative interpretation was possible. It appears 
that some of the points in the [1120] direction can 
be explained as transitions between the large 
mass levels corresponding to the resonance mass 
value m* ~0.2m,. These appear as reasonably 
well-defined lines. 

We again attempted to interpret the data, this 
time assuming interband transitions and using 
the theory for nonparabolic bands since on theo- 
retical grounds Eq. (4) should apply; however, 
the pulsed data and those in the following Letter 
indicate that the deviation from parabolicity is 
not as large as that predicted by this equation. 

It is possible that an analysis similar to the pres- 


FIG. 2. Model for bis- 
muth of the energy levels 
near the edge of the Brill- 
ouin zone in a magnetic 
field. For the low-mass 
directions, the spin split- €g 
ting of the Landau levels 
is determined from 
S etg™2™,/m*. The 
arrows indicate j states 
rather than pure spin 
states. 
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ent one, but modified to take into account the 
effects of other bands, would provide a better 
model for the quantitative treatment of the data. 
The effective Hamiltonian given in Eq. (2) can be 
readily extended to the case of a band interacting 
with a number of other bands, but the location 
and character of other bands near the zone edge 
in bismuth are not at present known. 

This experiment indicates the power of the 
method of interband magnetcreflection measure- 
ments in exploring the energy band structure of 
metals. Masses and g values for both conduction 
and valence bands can be studied. In addition, 
bands can be studied not only at the Fermi sur- 
face (as with de Haas-van Alphen effect and cyclo- 
tron resonance) but at energies above the Fermi 
energy; and information about lower lying bands 
can be obtained. Furthermore, the temperature 
dependence of the band parameters can be studied 
up to room temperature and possibly higher. It 
is interesting that the detection of transitions 
due to a mass as large as ~0.2m, is possible at 
room temperature and with transient pulse tech- 
niques. This lends encouragement to experiments 
in other metals with high carrier masses at low 
temperatures, with steady magnetic fields. By 
using magnetic fields of the order of 100 kilo- 
gauss and higher, this method could be applicable 
to the investigation of a number of metals. 

The pulse experiments logically suggested the 
possibility that this phenomenon could be observed 


at lower fields and at low temperatures with photon 


energies well above those used in the pulse exper - 


iments, thereby permitting the observation of 
interband transitions to higher quantum numbers. 
Such experiments were indeed successfully 
realized and are discussed in the following Letter. 
We wish to thank Dr. L. M. Roth and Professor 
G. F. Koster for valuable discussion in connec- 
tion with the theory for two interacting bands in 
a magnetic field. 
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INTERBAND MAGNETOREFLECTION IN BISMUTH. Il. LOW FIELDS 


R. N. Brown, J. G. Mavroides, M. S. Dresselhaus, and B. Lax 
Lincoln Laboratory, * Massachusetts Institute of Technology, Lexington, Massachusetts 
(Received August 4, 1960) 


Magnetoreflection measurements in bismuth 
have been carried out at liquid air temperatures 
and below with static magnetic fields up to 38 
kilogauss. The measurements were made using 
a magnetospectrometer similar to that of Zwerd- 
ling et al.,* except that infrared reflection tech- 
niques were used. Best results were obtained by 
sweeping the magnetic field at a fixed wavelength. 
This was done at wavelengths between 6 and 14 u, 
using an NaCl prism. 

Single-crystal specimens were grown from Bi 
of “99.999” % purity. Surfaces were obtained by 


cleaving at nitrogen temperature along a trigonal 
plane. A diamond saw was used to cut two sam- 
ples 2 mmx 3 mmx 25 mm from the crystal. The 
cuts were made so that the magnetic field parallel 
to the exposed trigonal surface was in one case 
along the bisectrix axis [1010], and in the other 
case along the binary axis [1120]. The surface 
was then polished with Buehler No. 3 grit to ob- 
tain a good reflecting surface. The surface 
strains were removed by electropolishing.? Back- 
reflection Laue patterns gave the crystal orien- 
tation and revealed minor strains. The samples 
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FIG. 1. Variation of power reflection with magnetic 
field on bismuth for photon energy Ae =0.1041 ev. The 
arrows indicate the field for which interband transitions 
of this photon energy occur. The magnetic field is 
along [1120] and the inner Dewar is at 4. 2°K. 














(kilogauss) 


were freely mounted with Teflon tape in a copper 
sample holder attached by a copper rod to the 
bottom of the inner Dewar which contained either 
liquid helium or liquid air.* 

A typical trace of the variation of power reflec- 
tion with magnetic field at constant photon energy 
is shown in Fig. 1. The line shape is character- 
istic of dispersion curves. The center of the 
straight portions of trace, indicated by arrows, 
is taken as the value of the magnetic field for 
which the infrared photon induces an interband 
transition between a magnetic level of the valence 
band to one of the conduction band. The resolu- 
tion was sufficient to identify lines with peak to 
peak intensities from 3% to 0.05%. The intensity 
of a given line was found to increase with increas- 
ing magnetic field and to decrease with increas- 
ing quantum number. 

Each interband transition appearing in the mag- 
netic field trace (as in Fig. 1) was identified, and 
the results plotted in terms of photon energy Ace 
vs magnetic field H, as in Fig. 2 [1120] and Fig. 3 
[1010]. These measurements were obtained near 
liquid helium temperature. Quantitatively, the 
results at liquid air temperature were essentially 
the same except that the lines sharpened and in- 
creased in intensity as the temperature was low- 
ered. For H along [1120], only the light-mass 
(~0.01m,) transitions are seen. For H along 
[1010], both light (~0.01m,) and heavy (~0.02m,) 
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masses are observed. Experimentally, the light 
masses are distinguished from the heavy masses 
by their relatively greater intensity and by the 
disappearance of some of the light-mass transi- 
tions when polarized radiation is used. 

For both orientations, the points identified 
with a given transition follow straight lines, par- 
ticularly for transitions with low quantum number. 
As in the pulsed-field experiments,‘ the results of 
Figs. 2 and 3 show a convergence of the straight 
lines for different transitions to an energy gap 
€g= 0.047 + 0.003 ev.° The value for €g was deter- 
mined in each case by drawing the best straight 
line through the experimental points associated 
with low quantum number. The lines drawn 
through the experimental points for high quantum 
number were joined to this value of the gap. The 
deviation from linearity for the higher transitions 
may be an indication of nonquadratic energy sur- 
faces at these high energies. 

In order to interpret the interband transitions, 
the model illustrated by Fig. 2, of the preceding 
paper was adopted. The energy of an interband 
transition is written as 


c c 
Ae= 1 i 
€ E+ (a, + Mw + 2B 8 ote H 


a wt v 
+ (nm, + 2)hur, t 2B Bac H, (1) 


MAGNETIC FIELD (kilogauss) 


in which c and v refer to the conduction and va- 
lence bands, respectively, geff is the spectro- 
scopic splitting factor, and 8 is the Bohr mag- 
neton. Using this expression, one can in prin- 
ciple determine from the four lowest transitions 
the masses and g values for both the valence and 
conduction bands. From Eq. (1), we find that the 
experimental results are consistent with approxi- 
mately equal masses for the two bands and g 
factors approximately as given by Cohen and 
Blount, ° 


Bagg 2M p/m". (2) 
For H in the [1120] direction, the value me" /m, 
= 0.010 is obtained for the mass of the conduction 
band. The effective mass of the valence band 
my*/m, differs from m,*/m, by less than 10%. 
For H in the [1010] direction, we find the effec- 
tive mass of the light electron to be m,*/m, 

= 0.0080+ 0.0002 with | (m¢*- my*)/m,| < 0.0003. 
For the heavy electrons in this orientation there 
are two possibilities. Either the valence and con- 
duction masses are approximately equal, with 
me*/m,=0.019 and |(me*-m,,*)/m,| <0.002; or, 
mc*/m, and my*/m, are quite different, m,*/m, 
= 0.019 and m,*/m,=0.011. This ambiguity can 
be resolved by using higher fields to look for the 
lowest observable transition (010!). The conduc- 
tion band masses are in agreement with those 
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deduced from Shoenberg’s’ de Haas-van Alphen 
data and Galt et al.* (cyclotron resonance). 

These results definitely indicate that the theo- 
retical expression, Eq. (4) of the preceding paper, 
does not give a good fit. The curvature for the 
theoretical curves of energy versus magnetic 
field is much sharper than that indicated in Figs. 2 
and 3. Extension of these measurements to higher 
magnetic fields will permit the study of the de- 
parture from parabolicity and the added resolu- 
tion will allow the independent determination of 
the masses for the valence and conduction bands 
together with the g values associated with them. 

We are grateful of Dr. G. B. Wright for his 
generous loan of the apparatus for carrying out 
these experiments, Mr. E. Warekois for x-ray 
studies of the samples, Mr. D. F. Kolesar for 
his capable assistance with the experimental 
work, and Dr. G. Dresselhaus and Dr. H. J. 
Zeiger for many fruitful discussions. 
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A series of measurements of the heat capacity 
of annealed niobium in the normal and supercon- 
ducting phases, carried out in this laboratory, 
was reported at the 1958 International Confer- 
ence on the Electronic Properties of Metals at 
Low Temperatures.’ One of the principal ob- 
jectives of the investigation was the determina- 
tion of the temperature dependence of the super - 
conducting electronic specific heat C,,, and its 
comparison with the detailed predictions of the 
Bardeen, Cooper, Schrieffer? (BCS) theory. 
Although an exponential dependence of Ces with 
temperature was observed as required by the 
theory, it was not possible to make measure- 
ments below 0.27, and hence the predicted 
change to the lower temperature exponential® 
could not be investigated. A new cryostat was 
therefore constructed for this and other investi- 
gations and the research on niobium was re- 
sumed. 

According to current ideas, the determination 
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DIFFERENCE IN LATTICE SPECIFIC HEATS IN THE NORMAL 
AND SUPERCONDUCTING PHASES* 


H. A. Boorse, A. T. Hirshfeld, and H. Leupold 
Pupin Physics Laboratories, Columbia University, New York, New York 
(Received August 18, 1960) 





of C,, requires a measurement of C,,, the spe- 
cific heat of the metal in the normal phase, and 
of C;, the specific heat in the superconducting 
phase. Values of C, were found to follow the 
relation: 


C =yT+A(T/e)5, 


with y = 7.62 millijoules/mole deg” and © = 231°K. 
Values of the normal lattice term, A(7T/6)°, 
were subtracted from C, at the same tempera- 
ture to yield values of C,, which, down to 1.7°K, 
were found to satisfy the following relation: 


C,,=7T [10 exp(-1.637 /7)], 











with 7, =9.09°K. It should be pointed out that 

this procedure assumes that the superconducting 
specific heat consists of independent lattice and 
electron contributions and that the lattice speci- 
fic heat is identical in the normal and supercon- 
ducting phases. 
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The validity of this method for determining 
Ces was discussed some time ago by Chester* 
who pointed out that it could be proved only for 
those superconductors which obey the similarity 
rule for critical magnetic fields. Critical field 
studies of isotopic mixtures are limited, but 
where they have been made it has been found that 
the similarity principle holds. On the basis of 
this evidence it seems to have been assumed that 
the principle holds generally and that therefore 
the phonon spectrum in the normal and super- 
conducting phases, at least for superconductors 
of several isotopes, is the same. Further, the 
current theoretical point of view as provided by 
the BCS theory presents no contradiction to these 
ideas. 

The recent specific heat measurements of 
Bryant and Keesom® on indium now indicate that 
the phonon spectrum in that superconductor at 
least is not the same in the normal and super- 
conducting phases. 


The purpose of this Letter is to present evi- 
dence that in niobium also the lattice specific 
heat is different in the two phases. When our 
measurements were taken at temperatures less 
than 1.7°K the total superconducting specific heat 
was found to be less than the corresponding lat- 
tice heat in the normal state as shown in Fig. 1. 
After another independent set of measurements 
had confirmed this result, a complete remeas- 
urement of the superconducting and normal heat 
capacities was made using as a thermometer a 
1-watt Allen-Bradley carbon resistor. This re- 
sistor was carefully recalibrated by comparing 
the resistance with (1) the vapor pressure of 
helium from 1.02°K to 4.2°K, (2) the transition 
temperature of a sample of very pure lead taken 
as 7.17°K,® and (3) the triple point of normal hy- 
drogen, 13.96°K. The results gave a complete 
confirmation of the 1958 measurements and also 
confirmed that C, at the lowest temperatures 
was smaller than the normal lattice heat in the 
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same range. Figure 1 shows the results of the 
four separate measurements below 7*=5.0. If 
less weight is given to the first and less precise 
set of data represented by the open circles, then 
below 7*=2.2 (T=1.48°K=0.1637,) the individual 
determinations fall along a straight line which 
passes through the origin. A possible interpre- 
tation would be that the electronic specific heat 
at this temperature is inconsiderable in com- 
parison with the superconducting lattice specific 
heat and that the latter can be represented by a 
T® term with an altered value of © which may be 
evaluated from the data as 243°K. 

We regard this interpretation as provisional 
since a number of small corrections must be 
made to these data to account for the thermal 
capacities of the heater wires wound around the 
sample, for the Allen-Bradley carbon resistor, 
and for the glyptal adhesive. These corrections, 
now being investigated, would all be subtracted 


and their influence, if significant, could only in- 
crease the discrepancy with the normal line. A 
detailed report of the measurements will be sub- 
mitted .n the near future. 

We wish to thank Mrs. Claire Metz, Mrs. Sa- 
rina Hirshfeld, and Dr. Truly C. Hardy for their 
kind assistance in taking these measurements. 
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CHANGES IN SUPERCONDUCTING CRITICAL TEMPERATURE PRODUCED 
BY ELECTROSTATIC CHARGING* 


R. E. Glover, III and M. D. Sherrill 
University of North Carolina, Chapel Hill, North Carolina 
(Received June 16, 1960; revised manuscript received August 29, 1960) 


In the process of charging a parallel plate con- 
denser, electrons are added to one plate and re- 
moved from the other one. Using a film 100A 


thick as one of the condenser plates, it is possible 


to add or subtract about one electron for each 10* 
atoms. The maximum charge that can be added 


is limited by the breakdown strength of dielectrics, 


typically on the order of 10° v/em. The effect of 
electrostatic charging on the conductivity of a 
metal film has been investigated repeatedly over 
the last fifty years.'~* Recent measurements by 
Deubner and Rambke* and by Bonfiglioli, Coen, 
and Malvano®® have shown that the normal-state 
conductivity of metals can be modified by electro- 
static charging (field effect). Changes of relative 
conductivity on the order of one in 10° have been 
produced. The conductivity of the films can be 
increased or lowered depending on the sign of 
the added charge. 

We have made measurements of the effect of 
charging on the normal-state conductivity of 
films of gold and have examined the effect on 
both the superconducting transition temperature 
and the normal-state conductivity of tin and in- 
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dium. Charging resulted in a change in the normal- 


state conductivity of all metals tried. It causeda 
change in the transition temperature of the super- 
conducting materials. 

Ruby (muscovite) mica substrates about 5 u 
thick were used for the experiments. A gold 
coating was evaporated onto the back side to serve 
as one plate of the condenser, and gold potential 
and current contact patches were placed on the 
front. The substrates were then placed in a va- 
cuum system and the film to be studied was 
evaporated so as to form the second plate of the 
condenser. Charging voltages of 150 volts were 
applied giving a maximum electric field of 3 10° 
v/em. The conductivity of the film was measured 
using a 1000-cycle bridge. Film thickness was 
estimated from the temperature-dependent part 
of the resistance. 

Of the metals examined, gold would be expected 
to have the simplest conductivity mechanism. Our 
gold films were condensed at room temperature. 
Negative charging resulted in an increase in con- 
ductivity while positive charging caused a de- 
crease. The relative change in conductivity was 
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of the same order of magnitude as the ratio of 
added or subtracted electrons to the number of 
atoms in the film. Measurements made at room 
temperature and at 77°K gave essentially the 
same results. The magnitude of the change in 
conductivity was the same for positive and nega- 
tive charging. 

Tin and indium films were prepared by evapora- 
tion onto mica substrates held at liquid nitrogen 
temperature. The films were annealed to about 
200°K. Measurements were made of the normal- 
state field effect as the films were cooled down 
to liquid helium temperature. The effect in both 
indium and tin was opposite in direction to that 
found for gold. Negative charging decreased the 
conductivity while positive charging increased it. 
For tin the effect was very much larger than fo 
indium. . 

For the superconducting measurements the 
film temperature was slowly lowered through the 
transition region. Resistance vs temperature 
curves for an indium and a tin film are shown 
by the dashed curves in Figs. 1 and 2. The next 
step was to stabilize the temperature at various 
points in the transition region and measure the 
change of resistance caused by charging. The 
results are shown by the solid curves of Figs. 1 
and 2. The displacement of the transition curve 
necessary to cause the observed resistance 
changes was calculated from the slope of the 
R vs T traces. 
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FIG. 1. Transition curve and change of resistance 
caused by negative charging of an indium film. Posi- 
tive charging causes a change in resistance of the same 
size but in the opposite direction. 


Qualitative measurements made in the super- 
conducting transition region for five indium films 
showed in every case an increase in resistance 
with negative charging and a corresponding de- 
crease with positive charging. This effect is the 
same in sign but much larger than that found for 
normally conducting indium. The measured re- 
sistance changes correspond to a decrease in 
transition temperature with negative charging. 
Quantitative measurements were made on three 
indium films. Thicknesses ranged from 60 to 
120 A. The results of Fig. 1 are typical. A field 
of 2.6x10° v/cm produced a shift in transition 
temperature on the order of 10°K. Approxi- 
mately 3x10 electron were added for each atom 
of the film. The amount of added charge was ob- 
tained from capacity measurements and the 
charging voltage. The number of atoms was cal- 
culated from the thickness estimate using the 
normal density of indium. To an accuracy of 10% 
the decrease in transition temperature caused 
by negative charging was equal to the increase 
caused by adding positive charge. 

The effect of charging on the transition tem- 
perature of tin was opposite to that found for 
indium. Qualitative measurements on seven films 
all showed an increase in transition temperature 
with negative charging. Quantitative results were 
obtained with four tin films having thicknesses 
from 70 to 110 A. Results typical of three of 
these are shown in Fig. 2. Here again the mag- 
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FIG. 2. Transition curve and change of resistance 
caused by negative charging of a tin film. Positive 
charging produces a change in resistance of the same 
size but in the opposite direction. 
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nitudes of the shifts of transition temperature 
with positive and negative charging were equal. 

It should be noted that the effect of charging on 
the resistance of tin films is opposite in sign 

in the superconducting transition region to that 
found for normally conducting material. The 
fourth tin film showed a shift in transition tem- 
perature an order of magnitude larger. For this 
film, the size of the effect was not symmetrical 
with respect to charge, the increase in transition 
temperature being about six times larger than 
the decrease. However, the width of the transi- 
tion region was also an order of magnitude larger 
than for the other films, perhaps indicating a 
lack of uniformity. 

After the first experiments on tin had been 
made but before measurements on indium were 
started, Donald Seraphim pointed out to us that 
an effect similar to the one observed might re- 
sult if the substrate were piezoelectric. Strains 
which would depend on the sign of the applied 
charge could cause changes in transition temper- 
ature. The space group Co9,(2/m) usually quoted 
for muscovite’ has a center of symmetry, a 
property that rules out the possibility of piezo- 
electric behavior. Unfortunately, a closer look 
at the literature indicates that the details of the 
muscovite structure are still in doubt since the 
proposed structure does not account for all of 
the x-ray diffraction lines.* Accordingly, ex- 
periments were made to try and detect piezo- 
electric behavior in muscovite. None was found. 
Strips of mica 3 mm by 19 mm were cut from 
material similar to that used for the field-effect 
measurements. Thin gold films were evaporated 
onto both sides of the strips. The specimens 
were then clamped at one end in a special vise 
which also made electrical contact to the gold 
films. The free end of the mica was observed in 
a microscope as fields of 3 x10° v/cm were ap- 
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plied first in one direction and then in the other. 
A piezoelectric behavior capable of causing a 
charge-dependent raising and lowering of the 
superconducting transition temperature would 
have shown up as a deflection dependent on the 
direction of the applied field. No such effect was 
observed. Any such strain (A///) was less than 
1x10. Measurements of the effect of pressure 
on the transition temperature of tin give AT. 
=30.0AV/V.° Assuming 3A//1=AV/V, the change 
in transition temperature that could be produced 
by a still undetected piezoelectric effect would 
be less than 1x107°°K. This is two orders of 
magnitude smaller than the shifts in transition 
temperature reported here. A further reason for 
believing that the effects are not caused by strains 
in the substrate is that the shifts are in opposite 
directions for tin and indium. This is in con- ) 
trast to measurements® which show that com- 
pression lowers the transition temperatures of 
both materials. 
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MAGNITUDE AND SIGN OF THE CONDUCTION ELECTRON 
POLARIZATION IN RARE-EARTH METALS 


V. Jaccarino, B. T. Matthias, M. Peter, H. Suhl, and J. H. Wernick 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received August 26, 1960) 


From nuclear magnetic resonance (NMR) and 
electron paramagnetic resonance (EPR) meas- 
urements on rare-earth intermetallic compounds, 
we have determined the magnitude and sign of 
the conduction electron polarization. These re- 
sults may be interpreted in terms of a negative 
exchange interaction of ~0.1 ev between the local- 
ized f electrons and conduction electrons. 

The NMR of Al’ in the rare-earth intermetallic 
compounds having the cubic Laves structure’ 
XAl, X =La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, 
Ho, Er, Tm, Yb, Lu, and as well, U)? was 
investigated as a function of field (8-14 koe) and 
temperature (4.2° - 300°K) using a Varian NMR 
spectrometer. The EPR of EuAl, and GdAl, was 
studied at 55 kMc/sec between 100° - 300°K. 

Large, temperature-dependent, Knight shifts 
(k) were observed for all of the 4f magnetic com- 
pounds CeAl, through TmAl, and for UAl,. The 
sign of k 4) obeyed the following simple rule: If 
we denote by S, L, and J the spin, orbital, and 
total angular momentum of the ground state of 
the rare-earth ion, then 


k ,,>0 for J=L-§, 


k <0 for J=L+8. (1) 
For the nonmagnetic compound LaAl,, ka) =+0.13 
+0.02%. The Knight shift results from the con- 
tact hyperfine interaction AI-S between the Al?” 
nuclear spin I and the s conduction electron spin 
s. 

An explanation for the magnitude, sign, and 
temperature dependence of k,, in the magnetic 
compounds may be found in the enhanced conduc- 
tion electron polarization that results from the 
exchange interaction between the localized f elec- 
trons and the conduction electrons; i.e., 


Hog =-5gf5°5: (2) 
By virtue of $5¢ an additional hyperfine inter- 
action, A’I-8, exists, which we assume for sim- 
plicity to be isotropic. The fact that J is a good 
quantum number necessitates projecting § on to 
J; hence A’I-§—A’I-J (§-J)\/J(J +1). We may 


then express the Knight shift at a temperature 7, 





Isf% GD 
k=k,| 1+ 5 , (3) 
8,8,"48 J(J +1) 


where k, is the Knight shift in the absence of 
f-electron polarization effects and Xy/n is the 
f-electron susceptibility (at T) per rare earth 
ion. 

Since 


(-3)<0 for 


(-H)>0 for +S, (4) 


the observed relations (1) for k are consistent 
with a negative 9,7. The values of 9, derived 
from the observed k and measured values’ of x 
using (2) are given in Fig. 1. The value of 
k,)(LaAl,) was assumed for kp. 

For PrAl, the temperature independence of 
sf was established by measurements between 
77° and 300°K. Since ka) (PrAl,:300°K) is but 
4 times k,, it was possible to verify that Rk is 
not proportional to x, as a function of 7, but is 
given by Eq. (3). For all but Ce, Pr, Sm, and 
Gd, the magnitude and field dependence of the 
linewidth at 300°K indicated that the Knight shift 
was anisotropic and could not be accounted for 
by electronic dipolar broadening. The large 
linewidths result from either an anisotropic ex- 
change interaction or partial quenching of the 
orbital angular momentum (anisotropic g) or 
both. 

We briefly mention certain other results of 
these experiments. Large Al*’ quadrupole inter- 
actions were observed, the magnitudes of which 
are given in Table I. The La‘*® NMR in LaAl, 
was studied as a function of temperature and ky, 
was found to be weakly temperature dependent 
(see Table I). The relative broadening of the 
La’*® and Al?” NMR vs either Gd or Dy impurity 
concentration in LaAl, was studied and found to 
be consistent with the Ruderman-Kittel- Yosida*»® 
conduction electron polarization description. 

Yosida® has shown that the conduction electron 
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FIG. 1. The magnitude and sign of the exchange 
interaction in the rare earth-aluminum intermetallic 
compounds of the cubic Laves structure XAl, as deter- 
mined from the Knight shift of the Al?’ NMR and the 
measured susceptibilities. The number of 4/ electrons 
and the corresponding valence is that which is usually 
given for the free ion. [C. Kittel, Introduction to Solid- 
State Physics (John Wiley & Sons, Inc., New York, 
1956).] X-ray data and susceptibility measurements 
support the assignment given in the figure for these 
quantities. The relative size of the anisotropic inter- 
action observed is indicated by the vertical line. 
Neither the inconstancy of §;¢ nor its magnitude and 
sign are understood at present. The anomalously 
small value of 3, -(SmAl,) derived probably results 
from the assumption that the state °H;» for the Sm** 
ion is pure. Due to the close spacing of the fine struc- 
ture levels and the importance of crystal field and ex- 
change interactions in admixing states for which J 4, 
this assumption is questionable. 





polarization gives rise to an electronic g shift 
of the EPR of the magnetic ions in metals. The 
fractional change in the g value expected is 


Ag/g=3n I op/2n BE > (5) 


Here, 2n/n, is the number of conduction elec- 
trons per atom, and ¢€¢ is the Fermi energy. 

The observed shift AH =Hye, -H, was only a 
fraction of the half-power linewidth (6H ~900 oe 
in GdAl,). The resonance line shows a Lorentz 
shape as would be expected in the case of strong 
exchange narrowing.® The observed line profile 
indicated a mixture of the real and imaginary 
parts of the rf susceptibility." The relative admix- 
ture was determined for each measurement. The 
field for resonance showed a temperature depend- 
ence which results from demagnetizing effects. 
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Table I. The Al’ nuclear electric quadrupole inter- 
action constants erqQ*" as determined from the field 
dependence of the broadening of the 1/2+—~ -1/2 tran- 
sition. The temperature dependence of the Knight shift 
of the La‘*® NMR in LaAl, is given in the lower portion 
of Table I. No measurable temperature dependence 
was found for the Al*' Knight shift in LaAl,: 2a) =0.13 
+ 0.02%. (Reference solutions: LaNO; and KAISO,.) 








Compound Valence 2qQ”" (Mc/sec) 
LaAl, 3+ 4.63 
CeAl, 3+ 4.54 
PrAl, 3+ 4.54 
YbAl, 2+ 2.30 
UAl, ? 4.31 
Temp(°K) 4.2 20 295 
Rya(LaAl,),% 0.65 0.64 0.60 








This dependence was eliminated by linear extra- 
polation to zero magnetization. 
Using these procedures, we obtained the values 


GdAl,: g= 1.9824 0.003, Ag= 0.010, 
EuAl,: g= 1.9944 0.003, Ag=0, 


where Ag is obtained from a comparison with the 
ionic value for Gd of g=1.992 estimated by La- 
croix.* A comparison with published measure- 
ments on ionic crystals® indicates that crystalline 
field effects on the g value are negligible. The 
measured Ag/g in GdAl, leads to a value of 9, 
=-0.03 ev if we assume €-¢=5 ev and n/n ¢=1. 
No polarization was detected within our accuracy 
in EuAl,. 

It remains in interpreting the NMR results to 
justify the statement [Eq. (3)] that the conduction 
electron polarization at the aluminum sites due 
to the exchange interaction with the rare earth 
ions is proportional to Iof as though the rare 
earth ion spins imposed a uniform magnetic field 
on the conduction electrons. In actual fact the 
s -f interaction is highly localized. Suppose that 
x(F -F,) is the conduction electron magnetization 
at r due to a localized “field” of unit strength at 
fF; (in this case the exchange field). Then the 
total magnetization at f is proportional to 


- - 


> a. > i(r anal r -)e 
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where the x(q)’s denote the Fourier coefficients 

in the expansion of x, assuming the usual bound- 
ary conditions. The sum over fr; gives zero unless 
q@ is zero or one of the reciprocal lattice vectors 
K. Thus the magnetization is 


> 


x(0)+ Dex (K)e""* 


The leading term is, of course, the magnetiza- 

tion due to a spatially uniform field of unit 

strength. Assuming that x is simply the Ruderman- 
Kittel susceptibility, x(q) decreases from x(0) 

to zero at g=2k;. Therefore, unless 2k; is com- 
parable with one of the K’s the sum will be neg- 
ligible, and Eq. (3) is justified. The internal 
consistency of the data suggests that no such 

special conditions prevail. 

Only for a naive model would one expect a de- 
tailed agreement between 9, as derived from the 
NMR and EPR experiments since the former 
gives the polarization at the Al site while the 
latter that at the Gd site. The fact that both 
measurements show a negative polarization is 
strong evidence for assuming 9,,< 0. 

To our knowledge this is the first determination 
of the sign of conduction electron polarization in 
a magnetic metal.'° 

We are particularly indebted to H. J. Williams 
and R. C. Sherwood for the measurements of the 
susceptibilities, to A. M. Clogston for many 


critical discussions, and to Mrs. V. B. Compton 
for x-ray analyses. 





‘The cubic Laves structure AB, belongs to the space 
group O;,"=Fd3m; the point symmetry for the A site is 
43m and the B site 3m. An axial quadrupole interaction 
is possible for the B site. All NMR measurements 
were made above the Curie point. 

*The compounds were prepared by induction heating 
of stoichiometric amounts of the pure metals in silica 
or alumina crucibles in an argon atmosphere. J. H. 
Wernick and S. Geller, Trans. Am. Inst. Mining, Met. 
Petrol. Engrs. (to be published). S. E. Haszko, Trans. 
Am. Inst. Mining, Met. Petrol. Engrs. (to be published). 
For the NMR and EPR experiments, the samples were 
pulverized and sieved (400 mesh). 

3H. J. Williams and R. C. Sherwood (private com- 
munication). 

‘m. A. Ruderman and C. Kittel, Phys. Rev. 96, 99 
(1954). 

°K. Yosida, Phys. Rev. 106, 893 (1957). 

Sp, W. Anderson and P. R, Weiss, Revs. Modern 
Phys. 25, 269 (1953). Assuming dipolar broadening, 
this theory predicts 6H = 500 oe in GdAl). 

™. Peter and B. T. Matthias, Phys. Rev. Letters 
4, 449 (1960). 

®R. Lacroix, Helv. Phys. Acta 30, 374 (1957). 

°J. W. Orton, Reports on Progress in Physics (The 
Physical Society, London, 1959), Vol. 22, p. 204. 

103, Owen, M. E. Browne, V. Arp, and A. F. Kip, 
J. Phys. Chem. Solids 2, 85 (1957), have found no 
shift for the EPR of Mn** in copper but have detected 
a broadening of the Cu NMR. 








OPTICAL OBSERVATION OF SPIN-ORBIT INTERACTION IN GERMANIUM 


Jan Tauc and Emil Antoncik 
Institute of Technical Physics, Czechoslovak Academy of Sciences, Praha, Czechoslovakia 
(Received August 12, 1960) 


Philipp and Taft’ have recently determined the 
optical constants of germanium from the meas- 
urement of reflectivity in a broad wavelength 
range. In the reflection spectrum they observed 
near 2 ev a characteristic peak which Phillips? 
from an analysis of the absorption curve ascribed 
to the transition L,:~ L, (see also Roth and Lax’). 
Our measurements of the reflection spectrum on 
etched single crystals of germanium have shown 
that this peak is split into two peaks, a at 2.1 ev 
and 6 at 2.3 ev (see Fig. 1). Using the optical 
constants as determined by Philipp and Taft! and 
considering the splitting a - b as a small pertur- 
bation of the reflectivity reproduced in reference 
1, it is possible to show that in the absorption 
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FIG. 1. Reflection spectrum of germanium (at 300°K) 
near 2 ev. 
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spectrum a similar splitting appears (naturally 
with the peaks somewhat shifted).* 

The double peak may be explained by the split- 
ting of the state L, by spin-orbit interaction. 
The observed value of about 0.2 ev agrees well 
with the theoretical estimate by Roth and Lax.* 

The correctness of this interpretation is sup- 
ported by measurements on Ge - Si alloys: the 
splitting diminishes with increasing content of 
Si (e.g., in an alloy with 22% atomic percent of 
Si the energy interval of the peaks a, 6 in the 
reflection curve is only 0.15 ev; in Si it was not 
observable at all).* Similar double peaks were 
also found’ in the reflection spectra of GaAs (the 
peak a is at 2.94, b at 3.20 ev), GaSb (2.00; 2.48), 
InAs (2.53; 2.82), InSb (1.82; 2.38); the general 
features of our interpretation are here in good 
accord with theoretical expectations as was kindly 
pointed to us by Phillips.® 

It may be expected that the corresponding spin- 
orbit splitting in the T point should also be ob- 
servable optically; however, in the absorption 
spectrum of Ge determined by Dash and Newman’ 


this is not apparent. This point deserves further 
experimental investigations. 

For help in our work and valuable discussions 
we thank A. Abraham, F. Lukes, J. C. Phillips, 
Z. Trousil, and B. Velicky. 





'H. R. Philipp and E. A. Taft, Phys. Rev. 113, 
1002 (1959). 

23. C. Phillips, J. Phys. Chem. Solids 12, 208 
(1960), 

3L, M. Roth and B, Lax, Phys. Rev. Letters 3, 217 
(1959). 

‘at our suggestion F. Lukes from the University of 
Brno determined the absorption index k of germanium 
by the method of D. G. Avery [Proc. Phys. Soc. 
(London) B65, 425 (1952)]; his preliminary results con- 
firm our observations. He found the peaks in the k 


curve at about 2.11 and 2.32 ev (private communication). 


5J, Tauc and A. Abrahim, Proceedings of the Con- 
ference on Semiconductor Physics, Praha, 1960 (to 
be published). 

83. C. Phillips (private communication). 

'w. C. Dash and R. Newman, Phys. Rev. 99, 1151 
(1955). 





PROTON HELICITY FROM A DECAY 


Robert W. Birge and William B. Fowler* 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received August 22, 1960) 


When A hyperons are produced with K° mesons 
in 7” -p reactions there is a large up-down asym- 
metry of the decay products with respect to the 
production plane.’ The angular distribution of the 
decay pion from a completely polarized hyperon 
at rest can be written as? 


dN =(1/4n)(1+a@ cosyx)dQ, 


where d& is the solid angle of the pion momentum 
vector P,, and x the angle between P_ and the 
spin of the hyperon. The constant a is given by 


a =2 Re(A*B)/(|A|?+1B1?), 


and characterizes the degree of mixing of parities 
in the decay. A and B are the amplitudes for de- 


cay into s,, and p,,, final states of the pion-nucleon 


system. The quantity aP, which has the possible 
values 0< |aP| <1, is a measure of the up-down 
asymmetry and has been experimentally shown* 
to be > 0.73+0.14. This large asymmetry can 
exist only if the A’s are highly polarized in the 
production process and if there is nonconserva- 
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tion of both parity and charge conjugation in the 
decay process. 

Another necessary consequence of parity non- 
conservation in the decay process is a longitudinal 
polarization of the decay proton from unpolarized 
A’s decaying at rest. It can be shown that this 
longitudinal polarization equals -a.‘* 

Fortunately, this longitudinal polarization of 
the proton, referred to the center of mass of the 
A, appears as a partial transverse polarization 
in the laboratory system when a A decays in 
flight and, hence, can be measured by a suitable 
scattering experiment. In this way the helicity 
of the proton can be obtained, whereas in the 
aP experiments only the lower limit to the mag- 
nitude can be determined. The sign of a was 
determined in an experiment of Boldt et al.,° who 
found a =+0.8579-}5, based on 54 selected events 
from a total of 257 in a multiplate cloud chamber. 

In the course of an experiment designed to 
produce = particles® from a high-momentum 
(1.1-Bev/c) K~ beam’ impinging on the Berkeley 
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30-inch propane chamber, about 20000 A’s were 
produced. A fraction of these (about 800 events) 
were observed to decay and have a subsequent 
scattering of the proton within the liquid of the 
chamber. This constituted a considerably larger 
sample of events than had been obtained by Boldt 
et al.; hence it seemed worth while to repeat the 
experiment because of its fundamental nature. 
Preliminary and incomplete results of this ex- 
periment have already appeared elsewhere.® 

All events visually identified as A’s were meas- 
ured, and the relevant data were then calculated 
on an IBM-650 computer. A constraint program 
was used to find the best fit to the data, taking 
into account transverse momentum balance and 
coplanarity of the A’s with respect to the produc- 
tion origin as well as the “Q” of the decay. In 
addition, the ionization of all tracks was visually 
checked for consistency with measured momen- 
tum and particle assignments. Above about 800 
to 900 Mev/c there was difficulty separating 6°’s 
from A’s, and those events were deleted. 

A major difficulty was encountered because the 
proton track, prior to scattering, was often too 
short to measure the momentum accurately. Be- 
cause very large errors were assigned to these 
tracks, the constraint program could not readily 
compute reliable values for the momenta, yet 
these values were needed to obtain the asym- 
metry parameter in the scattering process. For 
this reason, we assumed an elastic scattering 
and used the momentum obtained from the scat- 
tered prong instead. The events were accepted 
only when the proton momenta calculated by this 
method (after appropriate corrections for energy 
loss by ionization) were consistent with the values 
needed to give the right @ for the A. In most 
cases the scattered proton stopped in the chamber 
and a very accurate momentum determination 
could be made; unfortunately, however, these 
low-momentum events are just those where the 
asymmetry from scattering by carbon is quite 
small. Therefore these events do not help to 
measure the initial polarization of the proton. 

In order to eliminate possible inelastic scatter - 
ing events an acceptance cutoff was made in the 
scattering angle at a point where the elastic scat- 
tering cross section still dominates the inelastic 
processes. Even in those events in which both 

the incoming and outgoing momenta could be well 
measured, it was impossible to detect excitation 
into the low-lying levels of carbon; therefore, 

the cutoff procedure was used for all events. 

This procedure was necessary because the asym- 


metry parameter for inelastic scattering can be 
of opposite sign to that of elastic scattering in 
certain angular regions. Those events with two 
prongs that appeared to be hydrogen scatterings 
were subjected to an additional constraint pro- 
gram to determine whether they were indeed 
elastic hydrogen scatterings. 

As a result of rejecting @-like events, inelastic 
scattering, poorly measured events, events with 
large kinks in the tracks, and events in which all 
tracks were too short to be measured accurately 
or the momentum was obviously too low for analy- 
sis, only 212 events remained for further analy- 
sis. Finally, all single A’s without visible pro- 
duction origins were eliminated from the sample, 
leaving 183 events. Of these, 36 were cases in 
which the proton was scattered by hydrogen in 
the propane and the remainder were elastic car- 
bon scatterings. 

It is difficult to see how any bias can creep into 
the selection of these events even with such a high 
filtering factor, because the sign of the polariza- 
tion is not an obvious quantity when the event is 
viewed in a bubble chamber picture. 

The following procedure was then used to analyze 
the selected events. First, for each event the 
angle between the spin of the proton and the direc- 
tion of motion of the proton in the laboratory sys- 
tem was computed in a manner prescribed by 
Stapp.® The sine of this angle multiplied by the 
magnitude of the original longitudinal polariza- 
tion (-a@) is the magnitude of the transverse 
polarization. Second, correction was made for 
the precession of the proton spin in the magnetic 
field of the bubble chamber prior to the scattering 
event. The rate of precession of the spin was 
computed according to the equations derived by 
Ford and set forth in the paper by Nelson et al.,’° 
namely, 


eB-v 


= [1 +(4¢- ty]-2- 1)(3g -1)— — 


a 
s 








where B coma field, V=velocity of the pro- 
ton, and $g=2.79275. Third, because of spin- 
orbit forces, the scattered intensity is propor - 
tional to 1 +(-aP,)- P,(6), where -aP, is the 
transverse polarization of the incoming proton 
at the position of the scattering (as determined 
by the above transformation) and P,(@) is the 
asymmetry parameter in the scattering process. 
The direction of P. 2(8) | is along the normal to the 
scattering plane, 7), =kxk’/lkxk’|. The magni- 
tude of P,(@) is a function of incoming momentum 
and scattering angle and has been determined 
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experimentally elsewhere. For example, with 
an initial transverse polarization downwards, 
P,(@) is determined as 


P,(0)=(N,-N AND +N), 


where Np and N; are the number of protons 
scattered to the right and to the left. (The angle 

@ between P, and P, is 0 deg and 180 deg, respec- 
tively.) In. Figs. 1 and 2 the values of P,(@) are 
summarized as a function of laboratory-system 
momentum and scattering angle for protons 
scattered by carbon and hydrogen, obtained from 
references available in the literature and by pri- 
vate communication. A list of references from 
which these data were taken is available on re- 
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FIG. 1. Curves of constant asymmetry for elastic 
scattering of protons from carbon as a function of 
laboratory-system momentum and scattering angle. 
Dashed line indicates cutoff where elastic-scattering 
cross section still dominates inelastic scattering. 
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FIG. 2. Curves of constant asymmetry for proton- 


proton scattering as a function of laboratory-system 
momentum and scattering angle. 
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quest. Data from these charts were put into the 
memory of the computer and intermediate values 
obtained by interpolation. A dashed line in Fig. 1 
indicates the elastic cutoff. Below 300 Mev/c 
the asymmetry parameter P, was set equal to 0, 
even though at very low momenta P, again be- 
comes finite. This condition eliminated another 
63 events and left only 120 events of significance. 

Finally, the probability that a proton scatter to 
the right is Np/(Np+Nyz)= 3(1 - aP;), where 
P;=P,P, cos¢ for the ith event. The product of 
all of these independent probabilities is the like- 
lihood function 


L=J].(1-aP.), 
1 z 


which may be plotted as a function of a to obtain 
the most likely value for the magnitude and sign 
of a. 

In Fig. 3 the natural logarithm of L is plotted 
for 120 events. 

Although these data indicate that a = -0.45+ 0.4," 
the relevant question to be answered here is the 
sign of a and not the magnitude (which is already 
known to be > 0.734 0.14). Hence one may state 
that our value is three standard deviations away 
from a =+0.73. The ratio of the likelihood func- 
tion at a =-0.73 to that at a =+0.73 is 65:1; for 
a@ =¥0.45 the ratio is 12:1. The negative sign for 
a@ implies positive helicity for the proton in con- 
tradiction to theoretical speculations” based on 
the universal Fermi interaction, and to the ex- 
perimental findings by Boldt et al.,° whose re- 
sults are plotted in the same graph as a dashed 
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FIG. 3. Natural logarithm of the likelihood function 
L as a function of a. The longitudinal polarization of 
the proton from A° decay is -a. Solid curve is this 
work. Dashed line is that of Boldt et al.® This ex- 
periment indicates a=-0.45 +0.4. 
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line, and contrary to our own preliminary re- 
sults.® 

A second way of determining the sign of a is 
to measure the right-left asymmetry directly. 
The result obtained in this manner is (Np -N)/ 
(Np+Nz)= 20/120 =0.16+0.091. The computed 
value for this ratio, (-2a(P,P,)ay/m), based on 
the average of the product of the input polariza- 
tion -@P, and the asymmetry parameter P,(@), 
is 0.16 for a=1.0, and 0.07 for a=-0.45. The 
experimental asymmetry is in agreement with 
both of these values. 

Finally, we have computed the quantity D 
=-)) ;P; cos¢; as suggested by Bowen et al.,’* 
so that these data may be combined with other 
results. The value obtained here is D = -2.79, 
to which may be added Bowen’s -0.142. 

Every effort has been made to detect errors in 
the experiment, to the extent that the signs and 
magnitude of the polarization of the individual 
events making up the most significant part of the 
data have been checked numerous times. Although 
these data disagree with previous measurements 
of the sign of a, we have no other choice but to 
present the results as they now stand. 

It is with pleasure that we acknowledge the 
work done by Howard White and his group in pro- 
gramming the computer and processing all the 
data. Dr. Cyril Henderson collaborated with us 
on a preliminary version of this experiment. We 
are indebted to Dr. Wilson Powell and to all the 
members of the 30-inch propane chamber group, 
and to the staff of the Bevatron for their part in 
making a successful run of the chamber. 

We thank the many scanners and technicians 
for their help in the analysis. Discussions with 
Dr. Francis Muller and Dr. Oreste Piccioni were 
most helpful, as well as discussions with many 
other people too numerous to list. 





*Now at Brookhaven National Laboratory, Upton, 
New York. 
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Several authors have proposed the existence of 
new particles and 7-7 resonances.'~*® Along these 
lines we are studying the reactions 


He*+1n° (1a) 
Sade He*+w°® (1b) 
He*+ a+ +77 (1c) 
He*® + 7°+ 7° (1d) 
H°+2* (2a) 
p+d—< H8+wt (2b) 
H°+2*++7°, (2c) 


where w may be a particle of mass intermediate 
between that of a 7 meson and aK meson. Our 
first experiment consists of measuring at a fixed 
laboratory-system angle of 11.7 deg the He*® mo- 
mentum spectrum from 1.0 Bev/c to 1.6 Bev/c 
for incident proton energies ranging from 624 
Mev to 743 Mev. 

Figure 1 is a schematic drawing of the experi- 
mental arrangement. Protons extracted from the 
184-inch cyclotron impinged upon a Y-shaped 
gaseous deuterium target operated at a pressure 
of 320 psi and at liquid-nitrogen temperature. He* 
nuclei produced at an angle of 11.7+ 0.2 deg passed 
through the collimators C, through C, into the 
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magnet spectrometer and were detected by scintil- 
lation counters S, through S,. Q,, an 8-inch triplet 
quadrupole magnet, focused the particles initially 
at f,, and Q, (a lens similar to Q,) focusedthem 
finally at S,. The slit at f, was 0.5 in. wide at 
full proton energy of 743 Mev and 2.0 in at re- 
duced proton energies, thereby yielding momen- 
tum intervals of + 0.6% and + 2.5%, respectively. 

Identification of He* was made by requirements 
of momentum, time of flight, range, and dE/dx 
in S,. For He* momenta greater than 1200 Mev/c, 
S,S,S,8,5,5,5, coincidences were required, while 
below 1100 Mev/c, S,S,S,S,S, were required. The 
momentum interval between 1100 and 1200 Mev/c 
was measured with both arrangements to check 
that they agreed. Backgrounds were measured 
by using H, gas of the same stopping power as 
the D, in the target. The backgrounds thus meas- 
ured agreed with target-empty measurements 
when corrections were made for energy losses. 
At the full proton energy the backgrounds were 
always less than 1% of the maximum counting 
rates. However, at reduced energies they a- 
mounted to as much as 20% of the peak rate be- 
cause the greater dimensions of the incident 
proton beam due to scattering in the energy de- 
grader increased the production of He® in the 
walls of the target. 

Figure 2 shows the He* momentum spectra 
corrected for energy loss and for the change in 
momentum interval, Ap, with momentum p. The 
assigned errors are statistical and constitute 
the major part of the error. Each spectrum has 
a high-momentum peak due to single 7° produc- 
tion and a lower momentum continuum associated 
with double meson production. The solid curves 
drawn for the 7° peaks are the calculated mo- 
mentum- resolution functions of the system, nor- 
malized to the area under the experimental points. 

For comparing the continua with theory, we 
have chosen a simple phase-space calculation as 
a first approximation. The expression used for 
the phase-space volume element ¢, is 


¢ -€p__b;* (t-4” (3) 
S dp,dQ ws, r) F 








where p is the relativistically invariant three- 
body phase-space volume, w, and p, are the total 
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FIG. 2. He*® momentum spectra for various incident 
proton energies. (See text for description of curves.) 
Ordinate scale is correct to within a factor of two. 


energy and momentum of the He*, respectively, 
and ¢ is the square of the total energy of the two 
pions in their own c.m. system in units of the 
pion mass. The dashed curves drawn in the con- 
tinua are these calculations fitted by the method 


of least squares to all of the experimental points. 


It is clear that these calculations alone cannot 
reproduce the peaks that occur in the continua. 
We have investigated some of the conventional 
mechanisms for extending these calculations: 
for example, the effects of Bose statistics for 
the two pions, pion-nucleon interaction, and 
final-state wave function for the He® nucleus. 
Our crude calculations of these effects do not 
give quantitative agreement with the data. We 
have therefore attempted to fit only that part of 
the data outside the peaks with ¢,; the results 
are shown as the solid curves in Fig. 2. In 
Fig. 3, A, the differences between the data and 
the solid curves of Fig. 2, have been plotted as 
a function of He’ momentum. The curves are the 
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FIG. 3. Difference between He’ data and solid curves 
of Fig. 2 for various incident proton energies. (See 
text for description of curves.) - 


calculated resolution functions for a particle or 
a resonance of zero width and located according 
to kinematics for a mass or total energy of 310 
Mev. The areas under the curves have been 
normalized to the areas under the points. 

Some caution must be exercised before firm 
conclusions are drawn about the rather broad 
observed peak at full energy, as contrasted to 
the more narrow distribution expected. First, 
since the peak resides so close to the sharply 
changing edge of the continuum, the subtraction 
is sensitive to errors in momentum settings and 
the exact shape of ¢,. Second, the data have been 
combined over many runs, so that errors in set- 
tings from run to run tend to smear out any narrow 
peaks. We think that it is unlikely that these two 
possibilities are the entire explanation for the 
width observed. We therefore fitted the widths of 
the peaks with a simple Breit-Wigner one-level 
formula with experimental resolution folded in. 
The results of these calculations gave, for the 
natural linewidth, law =10+6 Mev. No definite 
conclusion on the linewidth can be drawn at this 
time. 

Information on the isotopic spin assignments 
for the particle or resonance may be obtained 
from a study of Reactions (2a), (2b), and (2c), 
which proceed via pure J/=1 production. Figure 4 
shows our attempts to measure these reactions 
at the highest proton energy. Because our spec- 
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FIG. 4. H*® momentum spectrum for incident proton 
energy of 743 Mev. Ordinate scale is arbitrary. 


trometer was designed for the He* measure- 
ments, we had to degrade the H® with Be toa 
momentum of 800 Mev/c. The small number of 
points and the large errors assigned are results 
of the large backgrounds and low counting rates. 
Experimentally, we find that the ratio of the 
cross sections of (2a)/|(2b) +(2c)] is about 2.3:1. 
We expect from charge independence the same 
ratio for the /=1 part of (1a)/[(1b)+(1c)] or (1a)/ 
peaks if the He* peaks are J/=1 and the remainder 
of the continuum is J/=0. From our fits, we find 
the ratio (la)/peaks to be 1.9:1. Thus, within 
the large errors involved, an J/=1 assignment 
for the peak is possible. 

Upon the suggestion of Chew we attempted to 
fit the full-energy He* data with a combination 
of an S-wave (J=0) ¢, and a P-wave (J=1) reso- 
nant bp given by 


$5 =¢4% -4)/t\IF 1’, (4) 


IF |? =[(¢ - t ? +(e -4)T?/t}", 


where (t,)'* is the position of the resonance in 
units of the meson mass and [I is a parameter 
appearing in the expression for the pion form 
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factor F, of Frazer and Fulco.* We obtained 
reasonable agreement for c, =3, ¢,=5.0+ 0.2, 
and ['=1.0 to 2.0. These values of ¢, and I do 
not agree with those of Frazer and Fulco, who 
predict ¢,, between 12 and 16 and I =0.4 in fitting 
the nucleon form factors and magnetic moments. 

In concluding, we can summarize the follow- 
ing: 

1. The data are inconsistent with the relativis- 
tically invariant phase space assumed. The dis- 
crepancy appears in the form of a narrow peak 
which kinematically behaves like a system with 
a mass or total energy of 310+10 Mev. The 
natural linewidth of this system is not more than 
about 16 Mev. 

2. Plausible explanations of the line are the 
existence of a new neutral particle or a resonant 
m-n system. The possible isotopic spin assign- 
ments are J/=1 or J/=0. The H® data crudely 
support the J/=1 assignment. 

3. Alternate explanations may be possible al- 
though we have not found any in quantitative 
agreement with the data. Further experimental 
work, particularly on the H® reactions, is clearly 
essential. Our present conclusions therefore 
must be regarded as tentative. 

A more complete discussion of the experiment 
and any new results will be forthcoming shortly 
in a more extensive article. 

We would like to thank Professor Geoffrey 
Chew, Professor Kenneth Watson, and Professor 
Emilio Segré for several enlightening discussions 
of the experiment and the results. We also wish 
to thank Mr. Morris Pripstein for his assistance 
during the early phases of the experiment, and 
Mr. James Vale and the cyclotron crew for their 
cooperation. 





*Work carried out under the auspices of the U. S. 
Atomic Energy Commission. 
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We have recently put into operation a new large 
double-focusing magnetic spectrometer capable 
of analyzing electrons or other singly-charged 






















energy of the scattered electron is high. 
We have begun to carry out a series of exper- 
iments on the proton at various energies and 










































































a particles up to a momentum value of 1000 Mev/c. angles with these spectrometers. In both instru- 
This spectrometer is of the 180° type previously ments we used Cerenkov counters as detectors. 
' used in this laboratory and has a mean radius of In. Fig. 2(a) we show a typical electron-scattering 
curvature of 72 in. We have employed this spec- a : , 
an trometer in electron-scattering work on the proton 
between incident electron energies of 650 Mev and = soo} ELECTRON SCATTERING | | 
900 Mev and between scattering angles of 45° and 7 AT 900 MEV | | , —/ 
nt 145°. A typical value of the solid angle employed 5 75° . | } | 
a in detecting electrons is 5.6x10~° steradian. The & 400 = oM"spectroweten | | 
| higher energies (> 650 Mev) have been realized $=CH, | al | 
by virtue of the recent extension in the length of < sop — ttt a Ws nin 
the Stanford linear accelerator. z | irs | 
The 72-in. spectrometer forms part of a two- 5 a L 4 aap | 
1 magnet system; the second part is the 36-in. z =n 3 
rly spectrometer described previously.’ Figure 1 8 | ro Pt 7 
is a schematic drawing of the two spectrometers a ae r | 4 
in a position in which they are 120° apart. Both 500 MEV 525 MEV 550 MEV 
nt spectrometers are arranged so that they can be Ws 108 08 OA Ot 
y rotated independently about a common scattering POTENTIOMETER READING 
center. We have taken data simultaneously with (a) 
both spectrometers, usually employing the 36-in. 
magnet at large angles, e.g., up to 145°. The = | | am 
cd 36-in. spectrometer can handle scattered elec- ae ’ 
sh trons only up to 500 Mev/c without excessive = : si 
oa deterioration of focusing and we have used the 3 — a asiamunl = | 
i 36-in. spectrometer in these experiments only > SPECTROMETER a ce 
a in the very safe region below 370 Mev/c. Such & mn / | 
a procedure limits, and has limited in the past, = 120 f+ | 
the ability of this spectrometer in obtaining = 
scattering data at high energies and small ang- = se aimee 
les,” i.e., in those circumstances where the 5 es 
° / + aa 
g / 
| SLIT BOx 1000 MEV CARRIAGE 2 } 
SPECTROMETER © aol I. >. 
ers - SeapPp os po o—- T An a ee 
sie ay MEV 
310 MEV 318 MEV 
| a 0 1 
foe — 132 134 136 138 14 142 
: ieee oneal POTENTIOMETER READING 
SCREEN |’ (b) 
sts FIG. 2. (a) Electron-scattering peak obtained with 
a 0.237-in. polyethylene target at an incident energy 
* of 900 Mev and a scattering angle of 75°. The target 
was at 45° with respect to the incident beam. (b) Same 
FIG. 1. The double spectrometer system used in except parameters are 850 Mev, 145°, 36-in. spec- 
these experiments. trometer, and target at 0°. 
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peak obtained with the 72-in. spectrometer at Tobie 1. Sammaxy of reseite. | 
an incident energy of 900 Mev and a scattering 

angle of 75°. The target material was polyethy- E r) g (do /dQ) x 108 

lene and we have subtracted the carbon back- (Mev) (deg) (fermi~*) (cm?/sr) 

ground in analyzing the data to obtain the proton- 

scattering peak. The carbon data were obtained 597 60° 6.96 67.9 | 
in the same run with a separate graphite target. 597 90° 11.20 8. 76 | 
In Fig. 2(b) we show a corresponding peak ob- road a rege “go 
— in the 36-in. spectrometer at 850 Mev and 700 60° 9. 16 38.0 

145°. 

The cross sections obtained in this work are 700 75° 12.01 11.2 
absolute cross sections and are based on the oa on me a 
readings of the Faraday cup shown in Fig. 1 750 45° 6. 86 136 
which gives the number of electrons in the main 750 75° 13.45 10.0 
beam. The secondary emission monitor was an , 
adjunct used as a monitor for the 72-in. spec- Ln Bd nos ed 
trometer and was calibrated during each run 750 141.5° 21.2 4 0.735 
with reference to the Faraday cup. The Faraday 750 145° 21.42 0.580 
cup was the same unit used in previous experi- 775 135° 21.86 0.644 
ments.” In our present experiments the vacuum 800 45° 7.70 104 
pipe leading to the Faraday cup was not used. 800 60° 11.53 23.5 

It is difficult at the present time to estimate 800 75° 14.93 8.19 
the errors associated with our absolute cross 800 90° 17.75 2.99 
sections. The main uncertainty is the detection 800 120° 21.64 0. 888 
efficiency of our Cerenkov counters. These 800 135° 22. 86 0.605 
counters have been tested carefully but it is still 800 145° 23.44 0. 360 
possible that a 10% or smaller error remains 850 45° 8.59 80.6 
in our efficiency measurements. Our efficiency 850 60° 12.77 17.8 
measurements differ very little from 100% and 850 75° 16.46 5.62 
we have therefore used this figure as the detector 850 120° 23.60 0.711 
efficiency. Statistical errors are usually small 850 135° 24.88 0.509 
compared with possible systematic errors and 850 145° 25.50 0.371 
perhaps a measure of our errors should be the 875 40° 7.56 127 
differences between individual measurements 67S 45° 9.05 82.0 
made under the same conditions of energy and 875 60° 13.42 15.7 
angle. These differences are always less than 875 145° 26.54 0.376 
+10%. Another small source of uncertainty is see 45° 9.61 61.5 
the value of the dispersion constant of the 72-in. =e pan 46.08 16.3 

: 900 75 18.03 5.35 
spectrometer. We have measured the dispersion 
constant to an accuracy of about 5%. 900 90° 21.24 2.09 

Our results are presented in Table I and in- 900 145° 27.58 0.347 
clude standard radiative corrections.’ In addi- 
tion to the cross sections in the Table we give 
the values of the square of the momentum trans- Unfortunately we have no firm data of that group 
fer. In many cases we have given a single value to compare with our results. 
as an average of several measured cross sec- We wish to acknowledge most gratefully the 
tions. As a general rule we believe our results cheerful cooperation and assistance we have re- 
are accurate to about + 10% because of the afore- ceived from D. Aitken, P. Auvil, C. Buchanan, 
mentioned possible errors. Relatively speaking, G. Burleson, B. Chambers, H. Collard, E. Dally, 
errors are probably less than + 5%. Our inter- P. Gram, T. Janssens, M. Ryneveld, and 
pretation of the above results will be given in an W. Wadensweiler. 
accompanying paper.® 

We understand that the Cornell group* is in- *This work was supported in part by the Office of 
vestigating the same electron-scattering problem. Naval Research and the U. S. Atomic Energy Com- 
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mission, and by the U. S. Air Force, through the 
Office of Scientific Research of the Air Research and 
Development Command. 

tVisitor from University of Strasbourg, Strasbourg, 
France. 

‘g, E. Chambers and R. Hofstadter, Phys. Rev. 103, 
1454 (1956). 

*R. Herman and R. Hofstadter, High-Energy Elec- 





tron Scattering Tables (Stanford University Press, 





Stanford, California, 1960), and the work of F. Bumil- 
ler and R. Hofstadter shown in Fig. 8 of this refer- 
ence. 
3R. Hofstadter, F. Bumiller, and M. Croissiaux, 
following Letter (Phys. Rev. Letters 5, 263 (1960)]. 
‘kK. Berkelman and J. Cassels (private communica- 
tion). 





SPLITTING OF THE PROTON FORM FACTORS AND DIFFRACTION IN THE PROTON* 


R. Hofstadter, F. Bumiller, and M. Croissiauxt 
Department of Physics and High-Energy Physics Laboratory, Stanford University, Stanford, California 
(Received August 25, 1960) 


Electron- scattering studies of the proton ob- 
tained in the last few years have been summa- 
rized recently. The measurements showed that 
the proton form factors (F,, F,) were less than 
unity, implying a finite structure, and lay ina 
region in which they were approximately equal 
to each other at momentum transfers (q) as high 
as q?=9.3 in- units of squared inverse fermis. 

At this value of the momentum transfer the 
measured ratio was F,/F,=1.23+0.20.? The ex- 
periments were confined to angles larger than 
60° at the highest energies then obtainable (650 
Mev) because of the limitation imposed by the 
energy-handling ability of the 36-in. spectrom- 
eter. It was therefore not possible to solve for 
F, and F, separately at values of q?>9.3. Sev- 
eral independent experiments”: * indicated that 
the F, values were slightly greater than the F, 
values at the same momentum transfer, but for 
simplicity and ease of calculation, in the past, 
the ratio of form factors was usually taken to be 


unity. 

We have now succeeded in splitting apart the 
two proton form factors. Because of the great 
interest in the proton form factors and because 
our data appear to be internally consistent, we 
wish to present in this paper some conclusions 
drawn from the experimental results given in 
the accompanying paper.‘ 

Our procedure has been to solve for the sepa- 
rate form factors (F,, F,) at conditions lying be- 
tween 7.7 <q? < 25 by choosing a pair of experi- 
mentally measured cross sections at the same 
value of q? but at different correlated values of 
energy and angle. We have used the method of 
intersecting ellipses® to find the form factors. 

Table I shows the values selected and the form 
factors found by combining the results. In a few 
cases, indicated by asterisks, we have used 
older data and combined the older values with 
the newly- measured cross section at the same 
value of g?. In two cases (866 Mev, 75°; 675 


Table I. Form factors F, and F,. 








q* (f-%) E, (Mev) 6; (deg) (do/dQ),(cm?/sr) E,(Mev) 6, (deg) (do/dQ),(cm*/sr) Fy, F, 
7.70 800 45° 1.04 107%! | 400 124° *1.06x 10" 0.520 0.490 
9.16 700 60° 3.80 x 10-2 464 135° *6.26x10-% 0.500 0.420 

11.50 800 60° 2.35 x 107% 500 135° *4.18x10-3 0.451 0.341 

14. 06 900 60° 1.43 x 107%? 597 120° 2.65x10-% 0.423 0.214 

16.97 866 75° 5.56 x 107% 650 135° 1.51x10-3 0.430 0.160 

18. 03 900 75° 5.35 x 107% 675 135° 1.23x10~3 0.451 0.108 

21.24 900 90° 2.09 x 10733 750 141.5° 7.35104 0.405 0.087 
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PROTON FORM FACTORS 


FORM FACTOR 


EXPONENTIAL 
0.80 f 





FIG. i. The proton form factors obtained in Table I, 
plotted against g?. F, may be approaching a diffraction 
zero. 


Mev, 135°) we have interpolated between two 
newly- measured results in order to obtain pro- 
perly matched pairs of cross sections. 

The form factor results now show the behavior 
plotted in Fig. 1. The dashed line is the form 
factor corresponding to the exponential model 
and F,=F,. Apparently our new F,, which is 
seen to approach zero, indicates qualitatively 
that the Pauli magnetic moment cloud is a “soft,” 
spread-out distribution. On the other hand, the 
constancy of F, suggests qualitatively that the 
Dirac electric/magnetic cloud has a small, per- 
haps point-like, core. 

The form factors found in the above manner 
were then put back into the well-known Rosen- 
bluth Eq. (40) of reference 1: 
do/dQ = Ong l@41F 1°+4 12" 1" 9 +a,.F,"} , () 
where the values of the coefficients a,,, a,,, and 
a,, are taken from the tables’ at the appropriate 
energies and angles. When this is done we ob- 
tain the results shown in Fig. 2. Notice that in 
Fig. 2(c) the cross section appears to be going 
through a diffraction dip, so characteristic of 
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electron-scattering studies on heavier nuclei. 
The experimental data indeed show this diffrac- 
tion dip and we believe that this is the first time 
diffraction has been observed in the proton. 

Within experimental error the new experi- 
mental results appear to be in agreement with 
the split form factor curves. It is very interest- 
ing to observe that the new form factors account 
for an increase of the cross section above the 
exponential case at small angles, merge approxi- 
mately with the exponential case at 120°, and 
drift below the exponential case at the large 
angles 135° and 145°. This is what the experi- 
ments appear to indicate and the result is a ra- 
ther complicated pattern of cross sections which 
the form factors must satisfy. The experimental 
data appear to fit the calculated curves for sepa- 
rate form factors absolutely as well as relatively. 

The data are in excellent agreement with the 
earlier experimental results. The measure- 
ments of a proton root-mean-square radius ap- 
pear to remain undisturbed because those meas- 
urements were made at low q values. However, 
we are aware that at higher values of q? the con- 
clusions about the neutron’s form factors may 
be influenced slightly.*»* This question is now 
under investigation by R. Herman and the authors. 
It may be pointed out that the inelastic electron- 
scattering studies on the deuteron should perhaps 
yield new information on the F, 1orm factor of 
the neutron when combined with these results. 

It is interesting to speculate on whether the pro- 
ton’s F, factor rises again after approaching 
zero at about g*=24 or whether it becomes nega- 
tive at that point. In our analysis we assumed 
F,=0 at q?>24. 

By use of these results new information on F, 
of the neutron should result from a study of the 
deuteron’s elastic scattering at large angles. 

We wish to thank Mr. Francis Lewis for his 
help in making some of the calculations with in- 
tersecting ellipses. 








*This work was supported in part by the Office of 
Naval Research and the U. S. Atomic Energy Commis- 
sion, and by the U. S. Air Force, through the Office 
of Scientific Research of the Air Research and Devel- 
opment Command. 

TVisitor from University of Strasbourg, Strasbourg, 
France. 

'R. Herman and R. Hofstadter, High-Energy Elec- 
tron Scattering Tables (Stanford University Press, 
Stanford, California, 1960). 

2F. Bumiller and R. Hofstadter; see Fig. 8, p. 28, 
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CROSS SECTION IN CM*/ STERADIAN 


90° 


SEPARATE Fs 
-—--- EXPONENTIAL 


600 700 800 900 
ENERGY IN MEV 


FIG. 2. Comparison of ob- 
served and calculated cross sec- 
tions. The experimental points 
are shown by hollow circles. The 
dashed line refers to the case F, 
=F, and corresponds to the form 
factors deduced from the old ex- 
ponential model.! The solid line 
is obtained from Eq. (1) and the 
newly-obtained form factors of 
Table I and Fig. 1. 


of reference 1. 

‘E. E. Chambers and R. Hofstadter, Phys. Rev. 
103, 454 (1956). 

F,. Bumiller, M. Croissiaux,and R. Hofstadter, 
preceding Letter [Phys. Rev. Letters 5, 261(1960)]. 
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5R. Hofstadter, Ninth Annual International Confer- 
ence on High-Energy Physics, Kiev, July, 1959 (un- 
published). See also reference 1, pp. 30-32. 
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BOUNDARY CONDITION MODEL FOR PROTON-PROTON SCATTERING 


D. P. Saylor, R. A. Bryan, and R. E. Marshak 
Physics Department, University of Rochester, Rochester, New York 
(Received July 18, 1960) 


An interesting feature of the potential model 
devised by Bryan’ for proton-proton scattering 
is that in the triplet states there is a region of 
great strength between 0.6%/uc (yu is the pion 


mass) and the radius of the infinite core, 0.38%/yuCc, 


and that outside this region the potential is com- 
paratively weak, and closely resembles most 
one- and two-pion theoretic calculations. These 
features made it seem likely that a boundary con- 
dition model with meson-theoretic potential tails 
could give a good description of nucleon-nucleon 
scattering in the 40- to 310-Mev energy range. 

The boundary condition model without potential 
tails was first suggested and applied to S-wave 
scattering by Breit and Bouricius.* Later, Fesh- 
bach and Lomon® extended the model to higher 
waves, but their predicted phase shifts resembled 
the No. 6 solution-of Stapp, Ypsilantis, and Metro- 
polis* which has since been ruled out. The pos- 
sibilities of the boundary condition model with 
potential tails were then explored by Tubis and 
Feshbach,° though the analysis has not been car- 
ried out to the extent of securing really good fits 
to the experimental data as yet. We have found 
that such a boundary condition model, using the 
Taketani-Machida-Ohnuma (TMO) potential,® in- 
deed yields a very good description of the p-p 
data up to 310 Mev consistent with the MacGregor, 
Moravcsik, and Stapp phase shift solution No. 1 
(MMS No. 1).” 

We have applied the boundary conditions at the 
same radius for all states of a given spin angular 
momentum S. Preliminary values of the boundary 
conditions were obtained by inward stepwise cal- 
culation of the wave functions for each state, 


Table I. Energy-independent boundary conditions 
used in conjunction with the TMO potential. The nota- 
tion is that of Feshbach and Lomon. # 




















Singlet 
jj Yy Pifuc) 
1S, 0.716 0.53 
'D, 1.52 0.53 
Triplet uncoupled 
j 1%, Fi/uc) 
3p, 19.2 0.56 
3p, 4.19 0.56 
oe -0.313 0.56 
Triplet coupled t 
Gj Gj G7 Moe) 
‘p, -*F, 1.32 4.71 4.54 0.56 
F, -5H, -0.644 22.5 -0.0815 0.56 





4See reference 3. 


starting in the asymptotic region with values cor- 
responding to the MMS No. 1 phase shifts at 310 
Mev. The energy-independent boundary conditions 
used in the model are given in Table I. Some of 
the resulting phase shifts appear in Table II and 
the predictions of this model are compared with 
experiment® in Fig. 1. It is seen that r, differs 
only slightly for the singlet (vr, =0.53%/yc) and 
triplet (ry, =0.56%/uc) states. Applying the bound- 
ary conditions at the same radius (7, =0.53//uc) 
for all states regardless of spin gives a slightly 
less good fit. 








Table II. Phase shifts predicted by TMO+ boundary condition model. The phase shifts are the nuclear bar phase 
shifts of reference 4 in degrees. 
E (Mev) 1S, 'D, 5P, 3p, ‘Pp, €o ¥, FF, ¥, & 5H, 
40 43.1 0.8 11.2 -7.3 3.8 -1.2 0.1 -0.2 0.0 0.0 0.0 
150 15.9 5.9 5.8 -17.5 14.6 -3.2 1.2 -2.3 0.7 -0.8 0.1 
210 5.3 7.6 -0.6 -21.6 17.0 -2.8 1.7 -3.0 1.4 -1.3 0.2 
310 -9.1 11.9 -11.5 -27.8 17.0 -1.6 2.2 -3.63 3.5 -1.5 0.4 
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vith 
"4 Because the computer program used in this difference in cross sections exhibited in Fig. 2(a). 
- d- work only takes the Coulomb force into account Figure 2(b) shows the extra strength of the TMO 
0) in an approximate way, this model is not com- potential in singlet states which brings about this 
* pared with low-energy experiments (less than larger phase shift. A defect of this model is the 
40 Mev) or with experimental data inthe Coulomb __ tendency of the cross section to dip too low in 
interference region at the higher energies. It the forward direction. As can be seen from the 
should be noted, however, that the boundary con- __ trend of the curves shown in Fig. 2, one can im- 
dition model will probably also give too large a prove the fit if one uses a potential which has 
hase scattering length (and a fortiori a negative shape more long-range attraction for the ‘D, state than 


\ 





parameter) a feature which has plagued most 
models to date.® 

The boundary condition method was also tried 
with the one-pion exchange potential (OPEP) in 
the outer region. Figure 2(a) shows the 150-Mev 
cross section with OPEP and TMO. The ’D, 
phase shift is 4° for OPEP and 5.1° for TMO 
when the 310-Mev MMS No. 1 value is matched. 
This difference is largely responsible for the 


is given by OPEP or TMO. 

Beyond the boundary radius 7,, the wave func- 
tions predicted by our model are very similar to 
those given by the Bryan potential. Evidently, 

a boundary condition model with static potential 
tails may simulate a more general velocity- 
dependent interaction than a potential model based 
on a combination of central, tensor, and spin- 
orbit forces. Furthermore, the present boundary 
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FIG. 2. (a) Plot of the do/dQ as predicted by OPEP and TMO models. (b) The 
singlet, even-parity OPEP and TMO potentials. 


condition model has the advantage, apart from 
its simplicity, of not requiring a knowledge of 
the meson-theoretic two-nucleon interaction be- 
low one-half of a pion Compton wavelength. This 
advantage may be considerable from the point of 
view of existing dispersion techniques. 





'R. A. Bryan, Nuovo cimento (to be published). 

’G. Breit and W. G. Bouricius, Phys. Rev. 79, 
1029 (1949). 

3H. Feshbach and E. Lomon, Phys. Rev. 102, 891 
(1956). 

‘H. P. Stapp, T. J. Ypsilantis, and N. Metropolis, 
Phys. Rev. 105, 302 (1957). 

5A. Tubis (private communication). See also A. M. 
Saperstein and L. Durand, III, Phys. Rev. 104, 1102 
(1956). 

®M. Taketani, S. Machida, and S. Ohnuma, Progr. 
Theoret. Phys. (Kyoto) 6, 638 (1951). 

™. H. MacGregor, M. J. Moravesik, and H. P. 
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RAPID REDUCTION OF COSMIC-RADIATION INTENSITY MEASURED 
IN INTERPLANETARY SPACE* 


C. Y. Fan,f Peter Meyer, and J. A. Simpson 
Enrico Fermi Institute for Nuclear Studies, University of Chicago, Chicago, Illinois 
(Received August 10, 1960) 


In a recent Letter’ we showed that the mecha- 
nism responsible for the rapid decreases of 
galactic cosmic-ray intensity was located beyond 
8 Earth radii, and that this result excluded all 
existing geocentric theories for this phenomenon. 
We have extended these experiments out to dis- 
tances of 8x 10° km and now can prove that there 
is no possible association of this phenomenon 
with the earth and its magnetic field. We also 
find strong evidence that there is a bulk motion 
outward from the sun of conducting solar plasma 
ejected in association with a solar flare which 
either carries magnetic fields within itself, or 
manipulates an interplanetary magnetic field to 
remove convectively the galactic cosmic radia- 
tion from a limited volume of the inner solar 
system. 

The experimental apparatus was carried in the 
space probe Pioneer V launched March 11, 1960 
on a trajectory inside the orbit of the earth and 
towards the orbit of Venus. The apparatus is 
identical with equipment on Explorer VI; namely, 
a wide-angle, triple coincidence counter tele- 
scope surrounded by 5 g/cm? lead which detects 
protons with energies in excess of 75 Mev and 
electrons of more than 13 Mev. Bremsstrahlung 
due to electrons in the energy range 0.2 - 13 Mev 
is measured independently. The details of the 
experiment have been described in earlier com- 
munications.‘»? To compare simultaneously the 
changes of cosmic-ray intensity at the earth and 
at the position of the vehicle in interplanetary 
Space, the recordings of nucleonic component 
intensity monitors were used. The telemetry 
limitations for Pioneer V confined the data re- 
cording to relatively short intervals several 
times each 24-hour period. Therefore, the cor- 
relation method used previously for Explorer VI 
data’ cannot be applied to Pioneer V results. 
Consequently, averaged, simultaneous observa- 
tions in the space probe and at the earth are used 
to determine the relative change in galactic in- 
tensity at the probe and at the earth both before 
the intensity decrease, and after the decrease 
reached its full extent. 

Although several intensity decreases have been 
detected in the Pioneer V data,* the largest and 


most readily analyzed event occurred on March 
31, 1960, at which time the probe was nearly on 
the Sun-Earth line and 5.2 x10° km inside the orbit 
of Earth. The changes in the nucleonic component 
at the earth—representing changes in the primary 
beam of particles with magnetic rigidities > 2.4 
Bv-are shown in Fig. 1. The corresponding data 
for Pioneer V are shown in Fig. 2(a). 

It is almost certain that the ejection of plasma 
from the solar flare* of importance 2 at 1455 - 
1858 U. T. on March 30 led to the commencement 
of the geomagnetic storm and the beginning of 
the cosmic radiation intensity decrease at about 
1200 U. T. on March 31. From this time delay 
we compute an average radial velocity of 2000 
km/sec for the advancing “front” of conducting 
plasma. The expected time difference of ~ 40 - 

50 minutes between the arrival of the plasma at 
Pioneer V and at the earth could not be observed 
due to the intermittent reception of data. 

After extrapolating the neutron intensity moni- 
tor data of Fig. 1 to the top of the atmosphere,’ 
we compute the ratio of the relative intensity 
change at the space probe and at the earth as 


Relative change in Pioneer V detector _ 1.34 0.15 
Relative change at the earth 


on April 1, 1960. 
The larger decrease measured by the Pioneer 
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FIG. 1. The changes in galactic cosmic-ray intensity 
for particles above 2.4-Bv magnetic rigidity. The 
neutron monitor is located at Climax, Colorado. 
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FIG. 2. Telemetered data from the space probe 
Pioneer V at distances 4 -5.5x 10° kilometers from 
the earth. The errors for each bar are less than 1%. 
The time for the first arrival of solar protons April 1 


was determined by the onset of enhanced ionization from 


protons at the polar cap.* The magnetometer meas- 


urements in Pioneer V are published by Coleman et al.’ 


V detector arises from the removal of cosmic 
radiation with magnetic rigidities below the cut- 
off for the detector at the earth (2.4 Bv) but not 
below the threshold of the probe detector. In- 
deed, this removal of low-energy particles per- 
sisted for more than 30 days after April 1 al- 
though the relativistic particle flux by that time 
had returned to the intensity prevailing before 
the decrease.® From analysis of our neutron 
monitor data extending to the geomagnetic equa- 
tor where the cutoff reaches ~15 Bv, we find 
independent evidence that the intensity decrease 
was more strongly dependent upon particle mag- 


netic rigidity than for the event of August 19 - 20, 


1959 previously reported.* 

These results prove that the rapid reduction 
in galactic cosmic radiation is a phenomenon of 
solar origin taking place in interplanetary space 
and is not related to the presence of the earth or 
its magnetic field. Since the observed intensity 
decreases only follow solar flares on the visible 
side of the sun, the reduction of cosmic-ray 
particle intensity is confined to a limited volume 
of the inner solar system. We also eliminate all 
electrostatic field deceleration models for this 
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effect. 

From the simultaneous measurements in Pio- 
neer V and at the earth we may also deduce some 
of the properties of these modulating magnetic 
fields: 

(a) It has been known for many years that the 
rate of decrease of cosmic-ray intensity may be 
as high as 5% per hour for some events. On the 
basis of our Explorer VI studies’ and the present 
results, we now conclude that this high rate ap- 
plies to the removal rate of particles and, hence, 
to the transient magnetic field parameters in the 
interplanetary medium. Since telemetered data 
from Pioneer V on March 31, 1960 did not pro- 
vide details on the rate of intensity change, we 
refer to the neutron data in Fig. 1 for these 
measurements. Approximately 28% of the par- 
ticle flux detected by Pioneer V disappeared in 
less than 20 hours, with half of this decrease 
occurring in less than six hours. From these 
data and the computed radial velocity of 2000 
km/sec for the advancing “front” of solar plasma, 
it is clear that a penetration to ~0.3 a.u. behind 
the front is required to reduce the intensity by 
14%-—half the full intensity decrease. For other 
events of this kind where the rate of decrease 
may be as high as 5 to 6% per hour the corre- 
sponding depth of penetration may be smaller by 
a factor of 3. This requires the rapid appear- 
ance of enhanced interplanetary magnetic field 
intensities; the magnitude of the changes and the 
degree of disordering of the required fields de- 
pends upon initial conditions and specific models 
which are not discussed here. Whereas our re- 
sults describe large-scale transient magnetic 
fields over great distances from Pioneer V, the 
magnetometer in Pioneer V registers field 
changes at the position of the vehicle perpendi- 
cular to its spin axis.*»’ These data’ are shown 
in Fig. 2(b) where magnetic field increases of a 
factor 10 - 20 times the quiescent field conditions 
occur at the time for which we deduce the exist- 
ence of large-scale magnetic fields from the 
changes in cosmic-ray flux. Both kinds of ob- 
servations show that magnetic fields are being 
moved or generated in interplanetary space as 4 
consequence of the solar flare on March 30. The 
only known way by which these transient fields 
could be established, or existing fields manipu- 
lated, is by moving, conducting plasma of solar 
flare origin. Therefore, we believe these Pio- 
neer V results provide the most direct evidence 
to date for the existence of conducting gas ejected 
at high velocity from solar flares—a concept 
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strongly supported already by many solar and 
terrestrial observations. Convective removal of 
galactic cosmic radiation by particle collisions 
with advancing magnetic field irregularities, 

such as a shock front, is the most likely explana- 
tion at this time. We do not discuss here the 
many models which may fulfill some of these 
conditions. 

(ob) By 0600 April 1 the advancing region which 
produced the full decrease of intensity had passed 
outward beyond Pioneer V and the earth. What 
were the magnetic field conditions remaining im- 
mediately thereafter in space between the sun 
and Pioneer V, i.e., behind the advancing “front”? 
This question may be investigated by examining 
the propagation of solar flare protons on April 1 
which penetrated the volume of space between 
the sun and Pioneer V only four hours after the 
full reduction of galactic cosmic-ray intensity. 
The solar flare was observed to begin at 0845 
and extend through 1222 U. T., with the first 
solar particles arriving to produce ionization in 
the polar ionosphere at 0945 U. T.* Thus the 
upper limit for the transit time is one hour. The 
solar particle increase is shown in Fig. 2(a) for 
protons > 75 Mev (but less than 1 Bev since no 
increases were observed by neutron monitors at 
the earth). Maximum intensity was reached 
within 50 minutes.* The subsequent decline of 
intensity with time ¢ closely follows a power law 
1/t. 

These data show that the interplanetary mag- 
netic field conditions behind the advancing front 
are either smooth and radial or weak, irregular 
fields (B,..,<5*10-° gauss), otherwise the low- 
energy solar protons would not penetrate to the 
orbit of the earth in less than one hour. 

Bremsstrahlung-producing radiation was also 
observed by our single, lead-shielded counter in 
Pioneer V° and by the University of Minnesota 
detectors® beginning about two days before and 
during the Forbush decrease. The measurement 


of detector efficiency is under way and these re- 
sults will be reported at a later date. 

For their contributions to the engineering and 
construction of our apparatus we wish to thank 
Mr. J. Jezewski, Mr. J. Lamport, Mr. L. Pe- 
traitis, and Mr. R. Takaki of the Laboratories - 
for Applied Science, University of Chicago. We 
also appreciate the over-all engineering and 
management of the Pioneer V project by Dr. A. 
Thiel, Dr. J. Lindner, and the staff of Space 
Technology Laboratories. The assistance of 
Dr. J. Lindsay of the N.A.S.A. for the Pioneer V 
project and the use of the Jodrell Banks radio 
telescope by Professor A. C. B. Lovell enhanced 
greatly the results of these experiments. 
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EXPERIMENTS ON THE ELEVEN-YEAR CHANGES OF COSMIC-RAY 
INTENSITY USING A SPACE PROBE* 


C. Y. Fan,? Peter Meyer, and J. A. Simpson 
Enrico Fermi Institute for Nuclear Studies, University of Chicago, Chicago, Illinois 
(Received August 10, 1960) 


The total intensity of cosmic radiation reach- 
ing the earth changes with the 1l-year solar ac- 
tivity cycle." The energy spectrum of relativi- 
stic particles up to at least 20 Bev changes ap- 
proximately from E~-**? (1954) to E~?*° at maxi- 
mum solar activity (1957-8) with the intensity 
being reduced at solar maximum by a factor 2 
for the relativistic particle flux,? and by a factor 
in excess of 4 for the integral flux down to ~ 100 
Mev for protons.* Figure 1 shows the change in 
nucleonic component intensity between 1954 and 
1960. Experimental evidence strongly supports 
the view that the spectrum at solar minimum 
approximates the full galactic spectrum outside 
the solar system, and that through an electro- 
magnetic modulation mechanism, initiated by 
solar activity, the galactic particle flux having 
access to the inner solar system, or the earth, 
is greatly reduced at solar activity maximum.” 
We have performed direct experiments which 
prove that any mechanism changing the cosmic- 
ray intensity over 11 years is centered about the 
sun, and the scale size is >1 astronomical unit 
(a.u.) for the volume of space in which the cos- 
mic-ray intensity is reduced at this period of 
the solar activity cycle. We also find that the 
gradient of cosmic-ray intensity near the orbit 
of Earth is so small as to suggest that modula- 
tion of galactic intensity occurs at distances 
much greater than the orbit of Earth. 
























OF | 
~~ 
2 | } | 
se 7 
= = Climax 
oe -10 L | | neutron intensity 
eo 
ae | 
oe L 
o” 6-15 | | 
8% | | M 
§ 8-20 | 
Oo | 
&  _25| Solar | | Solar | 
minimum | maximum 
4. 4. 4 \ 
1954 | 11956 | 1958 | 1960 


Years 


FIG. 1. Changes in the nucleonic component inten- 
sity between 1954 and 1960 at Climax, Colorado (~ 2.4 
Bv magnetic rigidity cutoff). 
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The experiment consists in carrying a cosmic 
radiation detector from the earth to great dis- 
tances in interplanetary space near the time of 
solar maximum. The detector measures protons 
above 75 Mev and, hence, integrates over that 
part of the cosmic-ray spectrum which changes 
by a factor of 2 to 4 in intensity between solar 
maximum and solar minimum. If the reduced 
intensity is restricted to the vicinity of the 
earth, the flux measured in interplanetary space 
will be greater by a factor of 2 to 4. If the vol- 
ume in which there is reduced intensity is helio- 
centric, then there will be an intensity gradient 
in the direction of the sun, the magnitude of 
which depends greatly upon the model adopted 
for the solar origin of the 11-year cosmic-ray 


intensity variations. 
Identical cosmic-ray detectors were carried 


by the satellite Explorer VI (launched August 7, 
1959) and the space probe Pioneer V (launched 
March 11, 1960), and are described elsewhere.*® 
Since this detector was not sensitive to brems- 
strahlung* from trapped electrons in the outer 
radiation belts of the earth, the cosmic-ray flux 
was measured from less than 10000 km to 48000 
km from Earth by Explorer VI. The data for 
August 9, 1959 are shown in Fig. 2. Between 
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FIG. 2. Counting rate of triple coincidence events 
in Explorer VI satellite (August 9, 1959) and Pioneer 
V space probe (March, 1960). The 13% correction in 
Explorer VI data arises from intensity changes known 
to take place outside the atmosphere between August, 
1959 and March, 1960 due to the 11-year cycle (see 
Fig. 1). Assigned errors on each interval of data are 
less than 1%. Pioneer IV (March, 1959) also observed 
constant intensity between 92000 and 658 000 km.° 
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this date and the launching of Pioneer V on 
March 11, 1960, the integrated cosmic-ray in- 
tensity at the earth increased by +13% as deduced 
from nucleonic component monitors. The August 
9, 1959 rates in Fig. 2 have, therefore, been in- 
creased by 13% for comparison with the Pioneer 
Vv data. We note that the cosmic-ray flux was 
approximately constant from inside the terres- 
trial field far into the interplanetary medium. 
Therefore the 11-year depression of cosmic-ray 
intensity is not geocentric. 

We know from analysis of nucleonic component 
data such as shown in Fig. 1 that the cosmic-ray 
intensity measured at the earth as the earth 
moves in its orbit about the sun is not signifi- 
cantly dependent upon solar system longitude. 
Thus, at solar maximum the cosmic-ray inten- 
sity is not only depressed within the inner solar 
system but the depression is roughly concentric 
about the sun in the ecliptic plane. 

Provided these intensity changes are brought 
about by solar modulation of the galactic flux, 
these results lead to the conclusion that there 
must be a cosmic-ray intensity gradient some- 
where between interstellar space and the inner 
solar system. This gradient will exist in the 
regions of space where the electromagnetic mo- 
dulation is taking place. 

On the one hand, if the modulating region lies 
beyond the orbit of Earth (>1 a.u.) then the gra- 
dient between Sun and Earth will be very small. 
A model in this category would be convective 
removal of particles by advancing magnetic field 
irregularities arising from instabilities of a 
solar wind.® On the other hand, if the earth lies 
inside the modulation region the cosmic-ray flux 
will decrease as one approaches the sun. Sever- 
al models, including a recent proposal of Elliot,’ 
require a large, negative gradient between Sun 
and Earth. 

One of the objectives for the space probe Pio- 
neer V was to examine this problem. At this 
time we have data over the initial 50 days repre- 
senting a change in solar radial distance of 0.1 
a.u. (10% of the distance to the sun) which sets 
upper limits on the radial gradient of cosmic-ray 
intensity in interplanetary space. To remove the 
time-dependent changes of cosmic-ray intensity, 
we determined the ratio of the intensities in Pio- 
neer V to the intensities at the earth represented 
by the nucleonic component monitor at Chicago. 
This ratio, normalized to unity at time of launch, 
March 11, is shown in Fig. 3 as a function of 
solar radial distance, and time. We have aver- 
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FIG. 3. The change of cosmic-ray intensity radially 
along the sun-earth line and relative to the intensity at 
the orbit of Earth. The errors assigned to these aver- 
aged intervals of data are less than 0.5%. The ar- 
rows indicate the existence of known perturbations in 
the low-energy portion of the cosmic-ray spectrum 
and data for these periods are not used in determining 
the intensity gradient. 


aged data over intervals of two to three days for 
this analysis. 

It is immediately apparent from Fig. 3 that 
there was a sudden loss of flux of nonrelativistic 
particles at the time of the rapid, Forbush inten- 
sity decrease of March 31- April 1.° This change 
in the primary spectrum persisted to the last’ 
available measurements 30 days later. The step- 
wise decrease in the very low energies was also 
detected in the Explorer VI data,* and will be 
discussed in a later publication. Smaller inten- 
sity decreases occurred at the times shown by 
arrows in Fig. 3. This effect is a change in 
spectrum which must be taken into account in 
measuring the gradient. There also exists the 
possibility that the negligible gradient observed 
from 0.95 a.u. to 0.90 a.u. (25th to 50th day) is 
the result of a space gradient plus a compensating 
gradual return of the low-energy particles. How- 
ever, such a coincidence of effects is not likely 
to continue over so long a period. 

We obtain a radial omnidirection intensity gra- 
dient of -(15+20)%/a.u. measured near the orbit 
of the earth and in the direction of the sun. This 
value is inconsistent with the existence of a posi- 
tive gradient such as required for an appreciable, 
steady solar production of the radiation. This 
small but negative gradient shows that (during 
the period of this measurement) any electromag- 
netic modulating mechanism required to account 
for the eleven-year intensity variations is located 
principally outside the orbit of Earth. Such re- 
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sults from Pioneer V place strong constraints 
upon acceptable heliocentric models for the 
eleven-year cosmic-ray intensity variation. 

We hope that additional data at later times in 
May may be extracted from the noisy data arising 
from the low signal level at these great distances. 

For their contributions to the engineering and 
construction of our apparatus we wish to thank 
Mr. J. Jezewski, Mr. J. Lamport, Mr. L. Pe- 
traitis, and Mr. R. Takaki of the Laboratories 
for Applied Science, University of Chicago. We 
also appreciate the over-all engineering and 
management of the Pioneer V project by Dr. A. 
Thiel, Dr. J. Lindner, and the staff of Space 
Technology Laboratories. The assistance of 
Dr. J. Lindsay of the N.A.S.A. for the Pioneer V 
project and the use of the Jodrell Banks radio 
telescope by Professor A.C.B. Lovell enhanced 
greatly the results of these experiments. 
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FORM OF THE ONE-PION EXCHANGE POTENTIAL* 


G. Breit, M. H. Hull, Jr., K. Lassila, and H. M. Ruppel 
Yale University, New Haven, Connecticut 
(Received August 18, 1960) 


The one-pion exchange potential (OPEP) has 
been used in calculations in nucleon-nucleon 
scattering by various authors employing the theo- 
retical expression as a whole.’ In connection 
with fits to data by means of adjustable phase- 
parameters it proved possible to make tests on 
the correctness of the mathematical form of the 
OPEP, and a brief report appears desirable. The 
tests were made by modifying the form of the 
OPEP through the introduction of additional terms 
multiplied by dimensionless parameters, g. Vary- 
ing one of the q together with the pion-nucleon 
coupling constant g, an adjustment was made to 
a minimum value of the mean weighted square 
deviation from experimental data, D, resulting 
in the determination of the most probable values 
of g, and of the g. Three types of variation have 
been studied. The first two, a and b, differ only 
in detail, both having the approximate meaning 
of changing the coefficient of the spin-spin term 
of the OPEP containing the factor (6,-6,) from 
the value 1 tol+g. For Gq this was literally the 
case. For q, the added term was (g,”/14)q, 
times the usual spin-spin term. In the units 
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used, g,*~14. The object in using two ways of 
introducing g was to provide a check. The two 

q are denoted below by g, when the distinction 
between a and 6b is immaterial. A third type of 
variation consists in the introduction of a central 
potential equal to q,. times (-1/3) of the singlet- 
even OPEP. The computation of D for various 
values of g, and of the q gave the following re- 
sults for the combined determination of g, and 
one of the g at a time. 

A preliminary form of p-p scattering data fit’ 
YLAM gave the pairs of values q, = -0.23+ 0.17, 
&" =14.0+1.3 and g- = -0.04+ 0.33, g,?=13.6+1.4. 
In this test the parameters not varied as OPEP 
phase parameters were all those for J <3 and in 
addition 67, 4» Pg, and K,. The p-p data and their 
analysis are more accurate than their n - p counter- 
parts and the above values may be supposed to be 
the more significant among those reported here. 

The n -p data fit? YLAN2M gave the pairs of 
values g,=0.23+0.19, g,?=12.841.8 and g,=-1.0 
+1.2, g,7=15.0+1.6. The one-pion parameters 
varied in this search were the same for T=1 as 
in YLAM while for T =0 the quantities not varied 
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were those with J <3 and in addition 6S, 9G, 
and p,. In this test the T=0 and T=1 parameters 
were varied together. 

The n-p data fit* YLAN3 gave the pairs of 
values g,=-0.15+0.18, g,’=13.341.8 and q, 
=-3.5+2.6, g,”=13.3+1.8 for variations of the 
OPEP with T=0 alone. The same fit was used 
also varying the OPEP parameters for T =0 and 
1 together. In this case the pairs of values g, 
=-0.12+0.13, g,?=14.2+1.8 and g, = -4.5+ 1.0, 

g,’ =11.8+ 0.8 were obtained. For test YLAN3 
the phase parameters which were not varied as 
OPEP quantities for T=0 were all those with 
J<3 and in addition those for *G,, *G,, and p, 
while for T =1 they were all those with J <3 and 
in addition *F, andp,. A later stage of YLAM 
was used for YLAN3 than for YLAN2M. In con- 
nection with the rather large q,. obtained in the 
last test, it is relevant that for it K, was included 
in the group of OPEP phase parameters sub- 
jected to the qg, g, variations. The possibility 
that the phase shift K, is not sufficiently well 
represented by the OPEP is suggested by the 
fact that the other tests described here do not 
give a definite indication of a nonvanishing q,. 

This explanation of the large q, is in agree- 
ment with the fact that the *G parameters change 
by amounts comparable to their OPEP values in 
the YLAN3 search. It should also be observed 
that K,, having the smallest L, gives the largest 
contributions to the amplitudes of the OPEP 
group. For this reason the apparent accuracy 
of the determinations of q,, is high in the searches 
including K,. The effect observed on q; is of 
approximately the magnitude expected from other 
evidence since -(1/3)(-4.5)=1.5 so that effectively 
the singlet-even long-range potential is roughly 
twice its OPEP value at the distances which mat- 
ter most in this work. 

A number of inaccuracies may be interfering 
with the validity of the numbers and error esti- 
mates given above. The data fits used are pre- 
liminary; there is an appreciable statistical 
error in the phase parameters of the non-OPEP 
group; the graphical determination of values of 
the q and of g, is not perfect and may be respon- 
sible for an inaccuracy of several units in the 
last digit of q,- Nevertheless the tests indicate 
the absence of large deviations from the OPEP 
form and furnish a confirmation of it. 

The results obtained are on the whole a cor- 
roboration of the tests of charge independence‘ 
which employed similar methods and of the ap- 
proximate estimates of the internucleon distances 


beyond which the OPEP may be assumed responsi- 
ble for the major part of the nucleon-nucleon in- 
teraction,® the OPEP groups of parameters being 
similar in the various sets of calculations. The 
procedure described above is capable of furnish- 
ing more definite information regarding the local 
ization of long-range deviations from the OPEP. 
One of the objects of the present note is to re- 
port the possibility of at least a partial empirical 
analysis of the form of the potential. In future 
work it may prove possible to test the range con- 
stants and through that perhaps to ascertain 
whether the usual g or f are the more fundamen- 
tal combinations. The tests themselves can be 
made to contain a more clean-cut spatial local- 
ization by including in the OPEP varied group 
additional parameters at the lower energies in 
accordance with the sensitivity curves calculated 
as in reference 5, and other extensions and im- 
provements are possible. It may be hoped, in 
this manner, to provide quantitative evidence 
regarding the two-pion exchange term which is 
furnished qualitatively in the present study. The 
absence of definite spin-spin and tensor force 
modifications and the presence of a central force 
addition indicated by the numbers reported on is 
in accordance with expectation for the two-pion 
exchange. 

The help of Mr. F. A. McDonald in performing 
some of the numerical reductions is gratefully 
acknowledged. 
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described in the last reference of reference 1. Further 
details regarding YLAM are contained in G. Breit, 
M. H. Hull, Jr., K. E. Lassila, and K. D. Pyatt, Jr., 
Phys. Rev. (to be published). A similar designation 
for n-p fits with YLAN as the first four letters has 
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K. Lassila, and K. D. Pyatt, Jr., Phys. Rev. Letters 
4, 79 (1960). 

3Fit YLAN3 is a development on fit YLAN2M with a 


major difference regarding the sign of p, through most 
of the energy range. It includes among other matters 
the employment of an improved mass treatment for the 
OPEP as described in the last reference of reference 2 
and difference in treatment of OPEP values for *F;, 
5F,, and p, at different stages of the search. 
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DEUTERON PRODUCTION IN HIGH-ENERGY COLLISIONS 


R. Hagedorn 
CERN, Geneva, Switzerland 
(Received August 29, 1960) 


The analysis of the beams of secondary particles 
resulting from 25-Gev protons hitting several dif- 
ferent kinds of targets shows the presence of high- 
energy deuterons.’ Deuterons have been observed 
also in cosmic-ray emulsion work’ and were in- 
terpreted as belonging to the high-energy tail of 
evaporation. 

Though the tail of evaporation and also the 
nuclear cascade may contribute to the total deu- 
teron yield, it is very probable that the CERN 
deuterons are mainly due to elementary nucleon- 
nucleon interactions. If that is so, then deuterons 
should be produced even in a pure proton-proton 
collision (and, by the way, not only deuterons 
but also H°, He*, and antideuterons and He‘, if 
namely one or two antinucleons are created’). 

It is the aim of the present note to explain this 
fact and give some numerical results. 

We consider a p->p collision and calculate the 
deuteron production by means of the statistical 
theory. It can be shown‘ that the statistical theory 
is derivable from S-matrix theory with the result 
that (n-particle end state) 


P,, =S(E,1,***N,) 


xR ,(E,T,T te? oe “+S ,m “7 (1) 


1 1 1 

is the probability of finding any one of the final 
states f contained ina set F. Here Ry is the 
total phase space of all feF including isospin, 
spin, and other statistical factors, whereas 
S(E,,***Np,) is the average with respect to F of 
the matrix element squared. E is the total energy 
(center-of-mass system) and n,:-~-n», are those 
parameters of the » particles which survive the 
integration over final states (e.g., coupling 
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strength, mass, isospin, structure) and there- 
fore influence still the averaged matrix element 
squared. Considering the inverse reaction leads 
to some knowledge of S: Namely, for the inverse 
reaction the individual matrix elements squared 
are the same, only the sum over final states of 
the reaction proper goes over into the average 
over exactly the same states for the inverse re- 
action, since initial and final states are now in- 
terchanged. The function S is by definition this 
average, and consequently 

Py nverS? Se, n ++, (2) 
On the other hand, the inverse reaction can take 
place only if the m particles (contained in a nor- 
malizing box of volume V) meet in a small volume 
of order 2 =(47/3)(1/u)°; 4: =meson mass (fh =c=1); 
hence S(E, m,+++N,) <(2/V)"~*. Inserting this into 
(1) and taking into account the Lorentz contraction 
leads to the usual formula for the statistical theory, 
which-—with certain refinements —seems to work 
quite well. 

Suppose now that a deuteron appears in the 
final state. Looking at the reaction proper it 
seems hard to believe that in 2 a deuteron (which 
is so much bigger) could be born and escape with- 
out being destroyed. But certainly it has its legal 
place in the possible sets of final states. Now 
things look better in the inverse reaction: a 
deuteron comes in and has to join inside © the 
other n-1 particles. The condition for this is 
twofold: It has to be there and, at the same time, 
it has to be as small as 2. Ina way this “melting 
in” is a measurement of the size of the deuteron 
(which of course destroys it—in fact it has to 
disappear as deuteron) and the probability of 
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finding it as small as 2 (or smaller) is given by 
the integral over the deuteron wave function: 


lp ,|?dVv. 
Ja 


Therefore, whereas each pion and nucleon is 
represented in S by a factor 2/V, the deuteron 
gives rise to a factor 


2, 2 
d_ 2 
ae ae dv. (3) 


Taking a Hulthén-type wave function with hard 
core (the latter has little influence) gives® 
lglvq\?dV = %-y. With this value for 27 one ! 
finds the results given in Tables Iand II. They 
show, though they do not apply directly to the 
experiments (no pp collisions have so far been 
analyzed with respect to deuterons), that this 
“elementary production” yields the correct orders 
of magnitude. The differences may be due to the 
presence of nuclear matter, to an anisotropy of 
nucleons in the center-of-mass frame (the trans- 
formation c.m..to lab involves the simplifying 
assumption of isotropy), and to “peripheral” col- 
lisions. 

It is satisfying that the above picture, which 
treats the deuteron as a quasi-elementary particle 
(at least makes no assumptions about how it is 
formed), can be supported by a kinematical con- 
sideration. It can be shown that the above formula 
can be made plausible also by asking the condi- 
tion that a neutron and a proton leave the inter- 
action region with a relative momentum, which 
is “acceptable” as a deuteron internal momentum 
(it belongs to the Fourier transform of yg). This 


Table I. Relative cross sections for deuteron pro- 
duction in pure p-p collisions. 








2.3 Gev 25 Gev 
af 
%W/mel, Sqr 0-0125 0.024 
~— 
Q,=7 9-006 0.012 
= a -5 
"den inet, Sg «9 0072 10 
at a +4 
Q =i 0.0036 10 
“atx “dean %q=8 
we + 1.73 not calculated 
d 





“This shows that at this primary energy a strong 
peak in the deuteron—as well as pion—c.m. spectrum 
should be observable. 


Table Il. The (deuteron/proton) ratio in the lab sys- 
tem at 15.9° in pp collisions with 25-Gev primary 


energy. 








Momentum of d Experiment 
and pinGev/e 2 ,=4 aa-t p-al* 
2 0.002% 0.001% approximately 
4 0.6 % 0.3 % constant, 2%,> 
6 2.8 % 1.4 % between 2.6 and 
8 7.8 % 3.9 % 5.5 Gev/c 
10 8.2 % 4.1 %& 





4See reference 1. 
Note that the experiment refers to p-Al collisions 
and the theory to pure p-p collisions. 


will be explained elsewhere. Even a detailed 
final-state interaction treatment® leads back to 
essentially the statistical formula with Qq as in 
(3). The phase-space integrals are computed 
rigorously (apart from the statistical errors of 

a Monte Carlo method’). All computing work was 
done on Ferranti-Mercury computers® (partly by 
the author but) mainly by Dr. W. Laskar, Univer- 
sity College, London. The author is very grateful 
to him for not only carrying out the machine runs 
but also writing up all data tapes. 

Many thanks are due to the University College, 
London, for offering computing time on their 
Mercury, to Dr. J. von Behr, CERN, for trans- 
forming the c.m. spectra to the lab system, and 
to many colleagues for discussions, in particular 
to J. von Behr, G. Bernardini, F. Cerulus, G. 
Cocconi, G. Fidecaro, A. W. Merrison, A. Pais, 
L. I. Schiff, and Y. Yamaguchi. The last-named 
was the first to obtain a numerical estimate on 
the basis of an argument involving restricted 
relative momenta for the outgoing nucleons. 





‘Vv. T. Cocconi, T. Fazzini, G. Fidecaro, M. Legros, 
N. H. Lipman, and A. W. Merrison, Phys. Rev. Let- 
ters 5, 19 (1960). 

2c. F. Powell, P. H. Fowler, and D. H. Perkins, 
Study of Elementary Particles by the Photographic 
Method (Pergamon Press, New York, 1959), p. 440. 

Calculations including these particles are in progress. 

‘R, Hagedorn, Nuovo cimento 15, 434 (1960). 

by, Yamaguchi (private communication). 

8... I. Schiff (private communication), and CERN 
Report No. 60-32 (unpublished). 

'F. Cerulus and R. Hagedorn, Suppl. Nuovo cimento 
9, 646 and 659 (1958). 

8R. Hagedorn, CERN Report No. 59-25 (unpublished). 
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PRODUCTION OF PARTICLE BEAMS AT VERY HIGH ENERGIES* 


S. D. Drell 
Institute of Theoretical Physics, Department of Physics, Stanford University, Stanford, California 





(Received August 25, 1960) 


We consider the production of beams of very- 
high-energy strongly interacting nuclear parti- 
cles. By studying the transition amplitudes near 
their poles corresponding to “real” one-particle 
intermediate states, we show that photons are 
very effective in initiating collimated beams of 
high-energy charged pions and K mesons, (anti-) 
nucleons, etc. A photon of energy & incident on 
a nuclear target produces a high-energy charged 
pion, say, of energy w,> $k (8 =c =1) in a for- 
ward cone of opening angle 6,.~/w, with a 
cross section that is reduced by roughly the fine 
structure constant, a =1/137, from geometric, 
~1/u?=20 mb. For nucleon-initiated processes, 
on the other hand, although one avoids the fine 
structure constant, the statistical model predicts 
that very high energy secondaries emerge in 
only a very small fraction of the collisions. 

This result is of significance for predicting 
and comparing yields from very high energy 
electron and proton accelerators. Experiments 
can be performed on existing machines to check 
the validity of this “pole analysis” which, as 
applied here to general inelastic processes, is 
an extension of the work of Chew, Goebel, and 
particularly of Chew and Low.’ 

Figure 1 illustrates the process being con- 
sidered, with X the nuclear target and (m) any- 
thing else produced in addition to the high-energy 
pion. Our basic approximation now, in calcula- 
ting the cross section when the pion is produced 
with energy w,~ and within the forward angular 
cone at an angle 6, <6,,.~ L/W, is to assume 
that the amplitude corresponding to Fig. 2 makes 
the major contribution to this process. In Fig. 2 
the photon produces a pair of charged pions, one 
of which emerges directly with W¢» 8 while the 





FIG. 1. Diagram for production amplitude y+X— 
m+(n), where (#2) denotes all particles produced in the 
final state, other than the observed 7, from the target 
: 


278 


other ploughs into X initiating reactions to various 
final states (m). The cross section of interest is 
the sum of squares of the amplitudes for all pos- | 
sible states (m). The above approximation is | 
made on the following grounds: 

(i) We are near the pole in the production amp- 
litude for Fig. 2 corresponding to a real pion in 
transit between (a2) and (6). The propagator for 
this pion is 


-[(k - q)? - u?]* 
= -kj7= - “1 
(2q-k) [2w e(1 8, cos? )] ; (1) 





where k,, = (k, k) and qp= (wg, q) are the photon 
and pion four-momenta, respectively. It is 
large for high energy Wg~k> ft and small angle 
6g<u /#q, equalling approximately 


{u°(e/w [1 + (0,0 /u ?}} =1/p?. 


This emphasizes the importance of this diagram 
relative to others for two or more pions or 
heavier particles propagating from (a) to (d). 
The next lightest state connecting (a) and (5) con- 
tains two pions and the continuous sequence of 
poles (the branch cut) for this state starts at a 
mass (2u)*, so that the propagator for the two- 
pion state then is less than 1/4u”. Another way 
of saying this is that the emitted high-energy 
pion is distributed over a cone of broader angle 
6’>2y /Wq for heavier intermediate states and 
so the flux per solid angle at Oq~U/wg is re- 
duced by a factor <1/4. 

(ii) We are led to a similar conclusion in con- 
sidering processes in which more than one pion 
emerges directly from the electromagnetic 








(n) 


FIG. 2. Diagram for one-pion contribution to the 
production amplitude y+X— 7+ (7). This dominates 
near the pole (k-g)* ~y’. 
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vertex (a). Integration over the undetected final- 
particle group (p) spreads the cone into which 
the observed pion emerges to a broad angle of 
width = [energy in particle group (p)]/wg > p/w¢. 

(iii) Lastly we may consider the amplitude cor- 
responding to one pion being exchanged between 
(a) and (b) but with the observed high-energy 
pion emerging from the nuclear part of the inter- 
action (b). We come to the same conclusion as 
above, namely that the cone into which the ob- 
served pion emerges is spread over a larger 
angle than p /wg. 

For these reasons we consider Fig. 2 a first 
approximation to Fig. 1 for w7~k>>y and 
998 u./wg. We can then compute the production 
cross section in terms of the pion electromag- 
netic vertex for a real photon multiplied by the 
total cross section for a pion of energy k - w 
incident on target X leading to all final states (m). 
This gives directly,” with Bg= 1g! /wq>1, 


& k,w ,0 
y, r#! q q 
a sin?6 dQ w (k-w )dw 
aed q qq q 4 
“22 (1-8 cos@ )? 4n RF 
n ( 5, q 





(k = )s 


*0,* +X, total (2) 


for the cross section that a 7* in energy interval 
dwg about w, and angle interval aQo about 9q is 
produced, along with anything else, by an inci- 
dent unpolarized quantum of energy & on a target 


Xx. 
Equation (2) leads to a sizable flux of pions of 


energy w,~k within a cone of angle @,,.= L/w¢. 
Let us compare it with other electromagnetic 

and nonelectromagnetic mechanisms for pro- 
ducing a 7 beam. First consider the electromag- 
netic production of 7+ pairs in the Coulomb field 
of a charge Ze (Fig. 3). The cross section ex- 
pressed in terms of the angle and energy distri- 
bution of one of the pions is® 








*~, at, a7 (, °q’ °” 
a [(u/w )*+sin*@ ] dQ k- 

alk (n/ q * q q w 4! wgdeg 
an (1-8) cosé |)" 4n k° 


x(Za)P 2% in(w /u). (3) 
q 
q 
In (3) we have neglected 1<In(wg/u) and taken 
the relativistic limit for the pions. Equations (2) 
and (3) are comparable in magnitude if in (2) we 


mone Aw 








FIG. 3. Diagrams for electromagnetic production 
of 7* pairs. 


uSe 07 +X, total=30 mb, as is appropriate for 
total pion-nucleon cross sections in the Bev 
range, while in (3) we insert Z =82. However, 

for any given target, process (2) will dominate 
for Og~u/w and Wg~k. Before leaving (3) we 
call attention to a difference in its derivation 
from that for (2). Both diagrams in Fig. 3 are 
necessary in order to arrive at the Pauli-Weiss- 
kopf cross section, whereas rescattering of the 
high-energy pion is not included in Fig. 2 and 

in the derivation of (2), on the basis of our earlier 
arguments. This is because the Coulomb scatter- 
ing of the virtual pion in Fig. 3 concentrates both 
final pions in the forward cone of 6,,.~pu/w and 
the two diagrams there are of comparable mag- 
nitude and opposite sign. 

Next consider production of pions by an incident 
proton beam on a hydrogen target. Hagedorn and 
von Behr* have calculated pion yields from such 
collisions on the basis of the statistical model; 
they predict the fraction of interactions leading 
to pions produced in given angular and momen- 
tum intervals. Their results are presented in 
graphical form in the CERN tables.* Since we 
are primarily concerned here with the produc- 
tion of high-energy collimated pion beams, we 
want to compare (2) and the statistical model for 
Oq~ L/W¢ and k~w,. Choosing k=25 Bev, w_=20 
Bev, and @g=1° to get an idea of the numbers, 
we find from (2) 


@o /d2Q dw =107? (5 Bev) 
Y,7 q q 


°" +X, total 


=0.3 mb/sr-Bev, 


for X a target proton. This is a factor of ten 
larger than the prediction of the statistical theory,° 
> ~4 = o- 

7x10 Op, total” Bev) =0.03 mb/sr-Bev. 
To be sure the predictions of the statistical model 
are least reliable at the high-energy tip of the 
spectrum of produced particles. It is also true 
that quantitative validity of (2) can be established 
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by experiment alone. Equation (2) may be pre- 
cise at the pole (1) in the nonphysical region 
cos6g = (1/8,,)>1, but we cannot theoretically 
anticipate whether or not there are appreciable 
corrections in the physical region at 6g~ L/w. 
At present the recent Cornell measurements® 
provide qualitative support for (2) in the energy- 
angle region of k~ 700-1000 Mev, Wg~ 500 Mev, 
and 64~ 15°. The peaking of their cross sections 
for mesons at forward angles is naturally ex- 
plained by the retardation denominator in (2). 
Also, the total yields as well as the ratio of the 
nm to the 1* yield can be approximately accounted 
for by known m+-nucleon cross sections occurring 
in (2). More detailed studies at larger Yq and 
smaller 6, are desirable.’ 

If the predictions of the statistical model for 
very high energy pion production in nucleon- 
nucleon collisions are verified, we conclude 
from (2) that a photon is more effective than a 
proton in initiating a high-energy beam of pions.*® 





aw (k-w )dw dQ 
q q q q 


Finally we turn to K-meson and baryon produc- 
tion. Equation (2) may be transcribed directly 
to the production of K* and K~ beams by changing 
the pion mass yu to ux for the K mesons, and 
the pion cross sections to ox+ ,y(k-Wg). Since 
K mesons are heavier than pions, several pions 
can accompany the K meson from (a) to (5) in 
Fig. 2 without appreciably broadening the cone 
of emission of the high-energy K*. In this ap- 
plication to K mesons then, Eq. (2) is to be given 
only qualitative significance in predicting the K 
beam intensity. It can be checked experimentally 
in terms of known total cross sections at en- 
ergies k-w,. The ratio of intensities of high- 
energy K* to K~ mesons is just the ratio of K~ 
to K* interaction cross sections. 

Going further and more recklessly afar from 
the pole, we can calculate also the production of 
nucleon and hyperon (anti-) particle beams by 
the mechanism of Fig. 2, obtaining in the high- 
energy limit 





#o, a: o”’ i Qn RS 


= © tall 2 oin2 ” 2) cin? 
(2 k/w )(a+n) sin oth sin 6 {2 (ew /2M ) sin 1 


dz °n+X, total” ~ “a 


(4) 





x 


where M is the mass of the particle n, a is its 
charge in units of e, and ) its Pauli anomalous 
magnetic moment in units of the Bohr magneton 
(a =+1 for the proton, = -1.91 for the neutron, 
etc.). We wish to point out in connection with (4) 
only that it predicts a considerable yield. For 
example, the cross section to produce anti- 
neutrons at 6 =M/w = 3°, with k =25 Bev and 
Wa = 20 Bev, fs : 


510756 


“1 -1 
N+X, total” Bev) sr-* Bev", 


and exceeds the prediction of the statistical 
theory calculation* for a proton-proton collision 
by a factor of =50. This is, however, a danger- 
ous and perhaps misleading application. 

Very valuable discussions, especially with D. 
Ritson and J. Ballam, have contributed to these 
remarks. 





*Supported in part by the U. S, Air Force through 
the Air Force Office of Scientific Research. 
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(1 - B, cosé )" 


, 





'G. F. Chew, Phys. Rev. 112, 1380 (1958); C. J. 
Goebel, Phys. Rev. Letters 1, 337 (1958); G. F. Chew 
and F, E. Low, Phys. Rev. 113, 1640 (1959). 

*The factors are arrived at as follows: 
$< s= (2c q)?= 81 AB g*wa? sin’, for the polariza- 
tion average of the pion electromagnetic vertex, times 
((q-k) -u"]~* = [2w (1-8, cos9g))~ for the square of 
the propagator, times d°q/ (21)*2w =B qwqdwgdQ,q/16r 
for the phase space of the observed pion, times 
k-\(k-wg)0 +X, total -g) for the absorption to all 
states wh of the intermediate pion, where the factor 
(k - wg)/k takes into account that the normalization 
factor for the incident photon is 1/2k, and 1/2(k- wg) 
for the pion is missing in extracting the cross section 
on +X, total -“g). 

3w. Pauli and ¥. Weisskopf, Helv. Phys. Acta 7, 
709 (1934); the formula in this reference is a factor 
of 4 larger than Eq (3) above. W. Pauli, Revs. Modern 
Phys. 13, 203 (1941); the result in this reference is 
a factor of 2 larger than Eq. (3) above. Equation (3) 
is derived for a point Coulomb field and is valid only 
for small angles @,<u/w, where the momentum trans- 
fers are ~y"/k. The Bethe-Heitler formula for the 
electromagnetic pair production of spin 1/2 Dirac par- 
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duc- ticles of the same mass p is larger by a factor ~6 at 8B. M. Chasan, G. Cocconi, V. T. Cocconi, R. M. 
ly bq™ Hiog due to the spin degree of freedom. Schechtman, and D. H. White, Phys. Rev. 119, 811 

i 4p. Hagedorn, Nuovo cimento 15, 434 (1960); J. von (1960). 
ing Behr and R. Hagedorn, CERN Report 60-20, 1960 TAn experiment is now in progress at Stanford Uni- 
: (unpublished) . versity, by L. Hand and W. K. H. Panofsky. 
— Cl 5pion and nucleon total cross sections are comparable ®Using the equivalent photon spectrum for a high- 
ons above a few Bev and both are ~ 30 -40 mb for target energy electron, one can directly compare the yields © 
1 protons. [See report of V. I. Veksler to Ninth Annual of pions from incident electron and proton beams. 
e High-Energy Physics Conference, Kiev, U.S.S.R., These calculations have been made by J. Ballam, 
= 1959 (unpublished); also private communication from Stanford University High-Energy Physics Laboratory 
iven J, Adams to B. Cork.] Report M200-8 (unpublished). 
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GENERALIZED BARDEEN-COOPER-SCHRIEFFER 
STATES AND ALIGNED ORBITAL ANGULAR MO- 
MENTUM IN THE PROPOSED LOW-TEMPERA- 
TURE PHASE OF LIQUID He*. P. W. Anderson 
and P. Morel [Phys. Rev. Letters 5, 136 (1960)]. 


In the paragraph after Eq. (8), we make the 
statements that “A, decreases with m, but that 
the total (binding) energy increases; m =I is 
always the best state.” More careful computa- 
tions show that this is not true. We find: 

(a) When normalized Yn! are used in (7), the 
binding energy and 4, vary in the same direc- 
tion. 

(b) The variation with m is more complicated, 

and apparently no very general rule holds. The 
best states are: Forl=1, m=1. Forl=2, 
m =1 and 2 have identical energies, lower than 
m=Q0. This agrees with reference 2. For /=3, 
m=2. One would guess, cautiously, that an in- 
termediate m is the best state in most cases. 

Since Fig. 1 is based on m =1, it, too, is not 
correct for 123. For 122 the density of states 
is linear in E, although the general appearance 
of curve 3 is maintained. 

The conclusions about aligned angular momen- 
tum remain valid. 


EVIDENCE OF THE NATURE OF ROTATING 
LIQUID HELIUM. J. R. Pellam [Phys. Rev. Let- 
ters 5, 189 (1960)]. 


In the first complete paragraph on page 190, 
the value of the apparent “rotation front” velocity 
radially inward should read “about 2 mm/min” 
instead of “about 2 mm/sec.” 
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~ In this section are printed the abstracts of Articles 
that have been forwarded to The American Institute of 
Physics for publication in THE PHYSICAL REVIEW. 
In quoting information obtained from this section be- 
fore the appearance of the corresponding Article, ref- 
erence should be made to “Physical Review (to be 
published)” rather than to this Journal. 











CARTESIAN TENSOR SCALAR PRODUCT AND 
SPHERICAL HARMONIC EXPANSIONS IN 
BOLTZMANN’S EQUATION. T. Wyatt Johnston, 
RCA Victor Company, Ltd., Research Labora- 
tories, Montreal, Canada (Received June 13, 
1960). 


The electron or ion distribution function, ex- 
panded in a sum of spherical harmonics, is 
shown to be equivalent to a tensor scalar product 
expansion. As such, it can be readily substituted 
into the Boltzmann equation to give transport 
equations integrated over angle and also the nec- 
essary equations for determination of the expan- 
sion coefficients. This has been done for terms 
up to order three, the order of the pressure 
transport tensor. 


NUCLEAR SPIN RELAXATION IN LIQUID HELI- 
UM 3. F. J. Low* and H. E. Rorschach, The 
Rice University, Houston, Texas (Received 
July 11, 1960) 


The longitudinal relaxation time of liquid He* 
has been measured as a function of temperature 
above 1°K and of magnetic field below 13 kgauss 
ina number of sample containers. At a tempera- 
ture of 2.0°K and in a magnetic field of 10 kgauss 
the longitudinal relaxation time, 7,, varied with 
the sample container from 60 seconds to 400 sec- 
onds. The transverse relaxation time, T,, was 
measured by a new method and was approximately 
30 seconds at a field of 10 kgauss in all sample 
containers. 7, was determined as a function of 
magnetic field at 2.0°K in a single sample con- 
tainer; the values increased from less than 50 
seconds in approximately zero field to 400 sec- 
onds at 13 kgauss. An impurity relaxation model 
is proposed to explain the T, results. By assum- 
ing both wall relaxation and a bulk relaxation 
given by the Bloembergen, Purcell, and Pound 
theory, the dependence of 7, on pressure and 
temperature can be quantitatively understood. 


The low values of T, are inconsistent with the 
Bloembergen, Purcell, and Pound theory and 
may be due to the presence of paramagnetic im- 
purities in suspension in the bulk liquid. 


*Now at Texas Instruments Inc., Dallas, Texas. 


RESONANCE ABSORPTION OF NUCLEAR 
GAMMA RAYS AND THE DYNAMICS OF ATOMIC 
MOTIONS. K. S. Singwi and A. Sjdjlander, Argonne 
National Laboratory, Argonne, Illinois (Received 
June 13, 1960). 


The theory of resonance absorption of nuclear 
y rays is generalized for an arbitrary system 
of interacting particles by expressing the rele- 
vant transition probability in terms of a space- 
time self-correlation function; and thus relating 
the resonance line shape to the incoherent differ - 
ential scattering cross section for slow neutrons. 
Two limiting cases, (i) a gas and (ii) a solid, 
have been considered. Discussion regarding the 
justification of the use of a classical self-diffusion 
function for a liquid is given and expressions for 
the broadening of the resonance line due to diffu- 
sive motions of the atoms of the interacting sys- 
tem are derived. It is suggested how a Mossbauer- 
type experiment could be used to give informa- 
tion regarding the diffusive motions of atoms in 
a solid and also, under more favorable circum- 
stances, in a liquid. 


MEASUREMENTS OF TEMPERATURES AND 
DENSITIES IN SHOCK-HEATED HYDROGEN AND 
HELIUM PLASMAS. Wolfgang Wiese,* H. F. 
Berg, and H. R. Griem, University of Maryland, 
College Park, Maryland (Received July 5, 1960; 
revised manuscript received August 8, 1960). 


Temperatures and densities in plasmas pro- 
duced by strong magnetically driven shock waves 
in a “T tube” were measured spectroscopically. 
Local thermal equilibrium exists within the ex- 
perimental accuracy for the conditions of this ex- 
periment. In both gases, the measured tempera- 
tures and densities are in disagreement with the 
values expected from the Rankine-Hugoniot re- 
lations, unless it is assumed that the conditions 
in the ambient gas are drastically affected by the 
initial discharge. Consistency between experi- 
ment and shock theory is achieved by considering 
magnetic pressure and dissociation and excitation 
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of the ambient gas by uv radiation from the dis- 
charge. 


*Now at the National Bureau of Standards, Washing- 
ton, D. C. 


ELECTRON PARAMAGNETIC RESONANCE OF 
EUROPIUM (II) IN SINGLE CRYSTALS OF 
CADMIUM SULFIDE. Paul B. Dorain, Brandeis 
University, Waltham, Massachusetts (Received 
July 5, 1960). 


Paramagnetic resonance absorption measure- 
ments have been made with single crystals of 
cadmium sulfide containing europium (II). The 
results can be fitted to the usual spin-Hamiltonian 
for *S state. The zero-field splitting is about 0.4 
cm™~. A comparison of the results with present 
theories for the crystal field splittings of °S is 
made. 


UNRESTRICTED HARTREE-FOCK METHOD: 
ELECTRON DENSITIES AND MAGNETIC FORM 
FACTORS FOR SPIN-POLARIZED Ni**. R. E. 
Watson,* Solid State and Molecular Theory Group, 
Massachusetts Institute of Technology, Cam- 
bridge, Massachusetts, and A. J. Freeman, 
Materials Research Laboratory, Ordnance Ma- 
terials Research Office, Watertown, Massachu- 
setts, and Solid State and Molecular Theory 
Group, Massachusetts Institute of Technology, 
Cambridge, Massachusetts (Received July 13, 
1960). 


The effect of relaxing the restrictions asso- 
ciated with the Hartree-Fock method are dis- 
cussed with particular emphasis on that constraint 
which requires common radial behavior for wave 
functions with all quantum numbers except ms 
(spin direction) in common. Results of such a 
“spin polarized” Hartree-Fock self-consistent 
field calculation are reported for the Ni*? ion 
and related to earlier calculations of Wood and 
Pratt, and Heine. Emphasis is placed on a con- 
sideration of the effects on the electron density 
and on x-ray and magnetic form factors. As is 
discussed, spin polarization of the 3d shell and 
the core results in an interesting effect on the 
magnetic form factor for this case. The calcula- 
tion suggests that one would obtain a magnetic 
form factor which is measurably expanded (hence 
a contracted charge distribution) in comparison 
with that appropriate for any single 3d electron. 
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Also presented are results of calculations of 
several hyperfine parameters which are in rough 
agreement with experiment. 


*Now at Avco, RAD, Wilmington, Massachusetts. 


CRYSTALLINE FIELD AND SPIN POLARIZA- 
TION EFFECTS ON ELECTRON DENSITIES AND 
MAGNETIC FORM FACTORS. R. E. Watson, * 
Solid State and Molecular Theory Group, Massa- 
chusetts Institute of Technology, Cambridge, 
Massachusetts, and A. J. Freeman, Materials 
Research Laboratory, Ordnance Materials Re- 
search Office, Watertown, Massachusetts, and 
Solid State and Molecular Theory Group, Massa- 
chusetts Institute of Technology, Cambridge, 
Massachusetts (Received July 13, 1960). 


The combined effects of spin (or exchange) 
polarization and an external crystalline field on 
charge densities, x-ray and magnetic form fac- 
tors, and hyperfine parameters are investigated 
following the analytic Hartree-Fock self-consistent 
field approach. The crystalline field was repre- 
sented by a crude cubic field arising from an 
octahedral array of point charges surrounding 
the central ion—in this case Ni**. In the strong- 
field approximation the atomic 3d electrons are 
“split” by the crystalline field and the spin polar- 
ization effect, resulting in a description of these 
electrons by a set of three distinct orbitals (each 
having different radial distributions and called 
togt, tagt ande,t). The ion’s spin density leads 
to a Fermi contact hyperfine term in better 
agreement with experiment than the value re- 
ported in an earlier spin-polarized calculation 
for the free Ni** ion and a magnetic form factor 
whose Fourier transform resembles none of the 
individual 3d charge distributions. 


*Now at Avco, RAD, Wilmington, Massachusetts. 


DYNAMICS OF RADIATION DAMAGE. J. B. 
Gibson, A. N. Goland, M. Milgram, andG. H. 
Vineyard, Brookhaven National Laboratory, 
Upton, New York (Received July 14, 1960). 


Radiation damage events at low moderate ener- 
gies (up to 400 ev) are studied by machine calcu- 
lations in a model representing copper. Orbits 
of knock-on atoms are found and the resulting 
damaged configurations are observed to consist 
of interstitials and vacancies. Thresholds for 
producing permanently displaced atoms (i.e., 
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interstitials) are about 25 ev in the (100) direc- 
tion, 25 to 30 ev in the (110) direction, and 
around 85 ev in the (111) direction. Collision 
chains in the directions (100) and (110) are prom- 
inent; at low energies the chains focus, at higher 
energies they defocus. Above threshold, the 
chains transport matter, as well as energy, and 
produce an interstitial at a distance. The range 
of (110) chains has been studied in detail. Local- 
ized vibrational modes associated with intersti- 
tials, agitations qualitatively like thermal spikes, 
ring annealing processes, and a higher energy 
process somewhat like a displacement spike have 
been observed. Replacements have been found to 
be very numerous. 

The configurations of various static defects 
have also been studied in this model. The inter- 
stitial is found to reside in a “split” configuration, 
sharing a lattice site with another atom. The 
crowdion is found not to be stable, and Frenkel 
pairs are stable only beyond minimum separations, 
which are found to be very much dependent on 
orientation. 


SPIN-WAVE SPECTRUM OF YTTRIUM IRON 
GARNET. R. L. Douglass, Department of Phys- 
ics, University of California, Berkeley, Cali- 
fornia (Received July 5, 1960). 


The spin-wave spectrum of yttrium iron garnet 
is treated using a Hamiltonian involving nearest- 
neighbor a-a, d-d, and a-d isotropic exchange 
interactions. Values of the exchange constants 
are estimated from the molecular field constants 
of Pauthenet. Anisotropy and magnetic dipole- 
dipole interactions are neglected. Twenty spin- 
wave modes are found, and their energies cal- 
culated at points of cubic symmetry in k space. 
The dispersion relation of the single “acoustical” 
spin-wave mode is found to agree with the value 
previously reported by Meyer and Harris. 


HALL EFFECT AND RESISTIVITY OF Ni-Pd 

ALLOYS. James A. Dreesen and Emerson M. 
Pugh, Carnegie Institute of Technology, Pitts- 
burgh, Pennsylvania (Received June 10, 1960). 


The two Hall coefficients and the resistivity of 
Ni-Pd alloys have been measured from 4°K to 
room temperature using fields up to 3.1 webers/ 
m*. The ordinary Hall coefficient is found to 
decrease in magnitude for small additions of Pd 
in Ni, but to increase as more Pd is added. It 


is also found that the ordinary Hall coefficient 
varies more slowly with composition for these 
alloys than it does for the Ni-Cu alloys. It is 
shown that a simple treatment successfully cor- 
relates the ordinary Hall coefficient, the resist- 
ivity, and the saturation magnetization of these’ 
alloys. The results indicate that the parallel 
half of the d band in pure Ni is not quite full at 
the absolute zero of temperature, and that the 
relaxation time is not a function only of the en- 
ergy for the thermal scattering in these alloys. 
The extraordinary Hall coefficient is found to 
obey none of the variously proposed theoretical 
relations. 


DETERMINATION OF THE FERMI VELOCITY 
BY THE HIGH-FIELD TILT EFFECT. Harold N. 
Spector, Institute for the Study of Metals, Uni- 
versity of Chicago, Chicago, Illinois (Received 
July 14, 1960). 


The attenuation of a sound wave in a semi- 
metal in the presence of a magnetic field is 
examined when the magnetic field is tilted from 
the direction perpendicular to the direction of 
propagation of the sound wave. It is found that 
in the limit of high magnetic fields (w,7 2 10q/) 
the behavior of the attenuation with the angle of 
tilt gives a point by point determination of the 
velocity on the Fermi surface. A detailed treat- 
ment of a two-spherical-band model is given and 
the case of a general Fermi surface is examined 
qualitatively. 


COLORATION OF PURE AND DOPED CALCIUM 
FLUORIDE CRYSTALS AT 20°C AND -190°C. 

W. J. Scouler* and A. Smakula, Laboratory for 
Insulation Research, Massachusetts Institute 

of Technology, Cambridge, Massachusetts (Re- 
ceived July 14, 1960). 


Coloration of pure CaF, crystals with 2.5-Mev 
electrons at room temperature results in bands 
at 580, 400, 335, and 225 mu. In YF,-doped 
crystals the bands are in the same spectral 
positions as in pure crystals. They are, how- 
ever, enhanced and their intensity ratios are 
significantly changed, the 400-my band being 
predominant. NaF-doped crystals show an even 
more drastic change: Coloration at room tem- 
perature produces bands at 605, 385, and 330 mu 
and is deeper than in either pure or YF,-doped 
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crystals. Since the 400-my band is strongly 
enhanced by YF, addition which forms F™ inter- 
stitials, it is correlated to neutral fluorine atoms 
in interstitial positions. Similarly, the 605-my 
band in crystals doped with NaF, which creates 
F~ vacancies, is correlated to electrons trapped 
in F~ ion vacancies (F centers). The enhance- 
ment of the 385-my band by NaF is explained by 
the reduction of the activation energy required 
for formation of interstitials and vacancies. Of 
several possible correlations the 330-my band 
is probably connected with a hole trapped in a 
Ca?* vacancy and that at 225 my with an electron 
trapped by a Ca”* interstitial. The bands at 440 
and 200 my which appear in NaF-doped crystals 
only must be connected with Na ions. Coloration 
at -190°C produces strong changes in the absorp- 
tion spectrum. In pure and YF,-doped crystals 
bands appear at 550, 320, and 270 mu, while in 
NaF-doped crystals bands are formed at 440, 
390, 315, and 200 my in addition to the 550-my 
band. The primary process of coloration at low 
temperature cannot differ from that at room 
temperature but the secondary processes are 
strongly influenced by temperature. 


*Present address: Lincoln Laboratory, Massachu- 


setts Institute of Technology, Lexington, Massachusetts. 


COLOR CENTERS IN ALKALINE EARTH 
FLUORIDES. D. Messner” and A. Smakula, 
Laboratory for Insulation Research, Massachu- 
setts Institute of Technology, Cambridge, Massa- 
chusetts (Received July 14, 1960). 


Single crystals of CaF,, SrF,, and BaF, have 
been colored by x rays (150 kv) and high-energy 
electrons (2.5 Mev). Five color-center bands 
in CaF, and SrF, and nine in BaF, have been ob- 
served in the spectral region between 0.175 and 
0.8 u. At room temperature the bands in CaF, 
are well resolved, while those in SrF, and BaF, 
overlap. At -190° and -255°C the bands sharpen 
but generally do not shift, except the 580-myu 
band in CaF, which shifts toward longer wave- 
lengths. Photochemical bleaching of the strong- 
est band shows a certain influence on the others. 
Three bands could be correlated to the Mollwo 
relation with the exponent 3.85. A possible cor- 
relation of the absorption bands to certain models 
is discussed. 

*Present address: Fa. Kalle A. G., Wiesbaden- 
Biebrich, Germany. 
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DISTRIBUTION OF ELECTRONS IN CADMIUM 
SULFIDE CRYSTALS. Ganesh P. Mohanty and 
Leonid V. Azaroff, Metallurgical Engineering 
Department, Illinois Institute of Technology, 
Chicago, Illinois (Received July 14, 1960). 


Very accurate x-ray diffraction intensities 
from a “pure” CdS crystal and one doped with 
chlorine were used to calculate electron density 
maps of their crystal structures projected along 
the c axis. To eliminate series termination 
errors, a difference map obtained by subtrac- 
ting the electron density of the “pure” crystal 
from that of the doped crystal was also prepared. 
According to this map, the electron distribution 
in the doped crystal is smeared out in the vicinity 
of the atom centers and an excess of electrons 
occurs in interstitial sites. This is contrary to 
the recently reported findings of Shuvalov who 
observed the formation of “electron bridges” 
along lines joining adjacent atoms. It is con- 
cluded that our observations are consistent with 
a partial disorder introduced in the doped crystal 
by the substitution of chlorine into sites normally 
occupied by sulfur atoms. Observed variations 
between corresponding x-ray intensities of the 
doped and undoped crystal can be similarly ex- 
plained by changes in the extinction effect due to 
this disorder. 


SUBLATTICE MAGNETIZATION AND RESO- 
NANCE FREQUENCY OF ANTIFERROMAG- 
NETS WITH LARGE UNIAXIAL ANISOTROPY. 
Bernard R. Cooper, Department of Physics, 
University of California, Berkeley, California 
(Received July 5, 1960). 


The effect of large anisotropy on the thermal 
behavior of the sublattice magnetization and 
antiferromagnetic resonance frequency is dis- 
cussed. It is shown that the curves for M(T)/M(0) 
and w(T)/w(0) vs T/T, are shifted upward from 
those predicted by a Brillouin function magneti- 
zation curve. The shifts are due to a combination 
of the effect of the large anisotropy energy with 
that of electron exchange correlation. Correc- 
tions to the molecular field approximation arising 
from pair correlations are taken into account by 
an extension of the Oguchi method. These cor- 
rections are expected to be roughly 1/Z as large 
as the more accurate corrections which would 
arise from correlation with all Z neighbors. 
This expectation is supported by a comparison 
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of the values of kT ./J and S(0) for the Oguchi 
method with those obtained by the Bethe-Peierls- 
Weiss and spin-wave methods, respectively. 

The results are related to the experiments on 
FeF, where the anisotropy is comparable to the 
exchange energy. 


RADIOACTIVE PREPARATION OF DEFECTS IN 
SOLIDS. John Lambe, Scientific Laboratory, 
Ford Motor Company, Dearborn, Michigan (Re- 
ceived July 5, 1960). 


An experimental study has been made of the 
use of radioactive materials in preparing defects 
in solids. Paramagnetic resonance was used to 
study such defects. The systems studied utilized 
tritium and krypton-85. Solid tritium was studied 
at 4.2°K. It was found that atomic tritium was 
formed and showed a characteristic hyperfine 
splitting of 543 gauss. The krypton-85 was also 
studied at 4.2°K. This should yield atomic rubidi- 
um but no paramagnetic resonance characteristic 
of this atom was found. 


DIFFUSION OF F CENTERS INTO POTASSIUM 
CHLORIDE SINGLE CRYSTALS. Hiroyuki 
Mizuno and Morio Inoue, Matsushita Electronics 
Corporation, Takatsuki, Osaka, Japan (Re- 
ceived June 13, 1960). 


Diffusion of F centers into potassium chloride 
crystals during an additive coloring process has 
been measured at various temperatures. A dif- 
fusion process which can be explained by a 
simple theory of volume diffusion seems to be 
predominant at higher temperatures. At lower 
temperatures, however, the occurrence of a dis- 
ordered interface diffusion with an F-center con- 
centration decreasing exponentially with the first 
power of penetration depth, is confirmed. 


ELECTRON LEVELS IN A ONE-DIMENSIONAL 
RANDOM LATTICE. H. L. Frisch and S. P. 
Lloyd, Bell Telephone Laboratories, Murray 
Hill, New Jersey (Received July 7, 1960). 


Let the potential of a one-dimensional scalar 
particle be V(x) =V,)- oo Ox -% 5), -o< x < 00, 
where V,<0, and where the sequence (x,;) is ran- 
dom, with a Poisson distribution. The quantity 
of interest is a certain limiting level distribu- 


tion, equal numerically to the node density of 
real solutions 7(x) of the Schrédinger equation. 
The random variables z j= v(x; - 0)/W(x5), 
-%0<j<oo, constitute an ergodic stationary 
Markov process. The stationary density T(z) 
of the (z,;) satisfies a first-order linear differ- - 
ential difference equation, and the node density 
is given (with probability 1) by lim,_.2?T(z) 
(Rice’s formula). Numerical results are ob- 
tained by integrating the second-order linear 
differential equation satisfied by the Fourier 
transform of T(z). 


THERMAL CONDUCTIVITY OF PURE INDIUM. 
R. E. Jones and A. M. Toxen, Research Center, 
International Business Machines Corporation, 
Poughkeepsie, New York (Received April 19, 
1960). 


The thermal conductivity of a pure indium 
specimen was measured in the normal and super- 
conducting states in the range of temperatures 
1.3 - 4.2°K. In the normal state the specimen 
showed sizable magnetothermal resistivity ef- 
fects which were not in agreement with Kohler’s 
rule for thermal conductivity. Near the transi- 
tion temperature the ratio of conductivities 
K,/Ky, exhibited the finite slope with temperature 
characteristic of electronic conduction limited 
by phonon scattering. The results were com- 
pared with a simple model proposed by Kadanoff 
and Martin and the agreement was found to be 
good. 


DEMBER EFFECT IN SILVER CHLORIDE. 
Alvin M. Goodman* and George Warfield, De- 
partment of Electrical Engineering, Princeton 
University, Princeton, New Jersey (Received 
June 20, 1960; revised manuscript received 
August 11, 1960). 


The Dember effect has been investigated in 
silver chloride using a chopped light-ac ampli- 
fication system of measurement. The effect has 
been resolved into a self-sustaining component 
which reaches a steady state (for particular 
wavelength and intensity of incident light) and a 
non-self-sustaining component which decays with 
time. The steady-state self-sustaining effect has 
been found to be of only one polarity (normal 
polarity) when the illuminated face of the sample 
is coated with silver, but of both polarities when 
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the illuminated face is not coated. The normal 
or positive polarity self-sustaining effect shows 
up at a shorter wavelength and the inverse self- 
sustaining effect at longer wavelengths, both in 
the vicinity of the optical absorption edge. A 
possible explanation for this behavior based on 
the existence of surface states is suggested. In 
the coated sample, the normal polarity effect was 
found to exhibit a peak at a wavelength of about 
0.383 yu at T=88°K. 

*Present address: RCA Laboratories, Inc., Princeton, 
New Jersey. 


OPTICAL MEASUREMENT OF THE PLASMA 
FREQUENCY AND M, , BAND OF CHROMIUM. 
Norman N. Axelrod* and M. Parker Givens, In- 
stitute of Optics, University of Rochester, 
Rochester, New York (Received July 7, 1960). 


The transmission of thin chromium films has 
been measured in the extreme vacuum ultra- 
violet. The onset of transmission at 500 ang- 
stroms or 24.8 electron volts agrees well with 
the theoretical prediction by Pines of 25 elec- 
tron volts; it also agrees well with the experi- 
mental electron eigenlosses measured by 
Watanabe, and Marton and Leder, of 26 elec- 
tron volts and 22 electron volts, respectively. 
The M, , band was also measured. 


*Now at National Aeronautics and Space Administra- 
tion, Goddard Space Flight Center, Washington, D. C. 


CALCULATION OF THE CRYSTALLINE FIELD 
STRENGTH: CHROME ALUM. A. J. Freeman, 
Materials Research Laboratory, Ordnance 
Materials Research Office, Watertown, Massa- 
chusetts, and Solid State and Molecular Theory 
Group, Massachusetts Institute of Technology, 
Cambridge, Massachusetts, and R. E. Watson,* 
Solid State and Molecular Theory Group, Massa- 
chusetts Institute of Technology, Cambridge, 
Massachusetts (Received July 1, 1960). 


A calculation of the crystalline field strength, 
Dq, is reported for the case of chrome alum, 
using the same model and crystal field potential 
employed by Kleiner and recently determined 
Hartree-Fock wave functions for the Cr** ion. 
On this basis, earlier theoretical attempts at 
determining, within the framework of crystal 
field theory, accurate Dq values are reviewed 
and analyzed. Particular emphasis is placed 
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on a consideration of reported point-charge cal- 
culations and Phillips’ method of including the 
effects of orthogonalization of ligand to metal- 
ion wave functions. Our results indicate that the 
point-charge model estimates for Dq gave good 
results mostly because they were based on the 
use of improper 3d wave functions for the transi- 
tion metal cation. Kleiner’s result is significantly 
improved—his wrong sign for Dq is reversed and 
a small positive Dq is obtained—but in such a way 
as to contradict Phillips’ conclusions. A dis- 
cussion is given of the various evidence for the 
inadequacy of the electrostatic potential theory 
and some of the necessary modifications are in- 
dicated. 


*Now at Avco, RAD, Wilmington, Massachusetts. 


EXPERIMENTS ON CHARGE-CHANGING COLLI- 
SIONS OF LITHIUM IONIC AND ATOMIC BEAMS. 
S. K. Allison, J. Cuevas, and M. Garcia- Munoz, 

Enrico Fermi Institute for Nuclear Studies, Uni- 
versity of Chicago, Chicago, Illinois (Received 

July 5, 1960). 


Beams of Li* ions accelerated to kinetic ener- 
gies in the range 10 - 475 kev were brought into 
charge equilibrium in the gases H,, He, and N, 
and the fractional amounts of Li’, Li®, Lit, Li**, 
Li*** measured. Because of the relatively large 
amount of Li in an equilibrated beam below 40 
kev in energy in propane and nitrous oxide, lim- 
ited investigations were carried out in these 
gases. The anomalously large Li** yield in heli- 
um at Li’ energies below 100 kev was noted. 

By holding each separated charged constituent 
successively in an orbit in a magnetic field, ad- 
mitting gas, and observing the beam attenuation, 
the total cross section for all charge- changing 
collisions was observed. Through application of 
the differential equations for growth or decay of 
a charged component, certain individual charge- 
changing cross sections can be deduced from 
these observed sums. Values are given for (0, 
+031), Sor» 0» 2» OC,» and (05.+05,) in H, He, 
and N, target gases throughout the energy range. 
Upper limits can be assigned to 0,3, 0,7, 0% 3, and 
0,,. Use of data on o,, in H, from other sources 
provides values of 0,; in that gas. Values for 0, 
for Li** in helium gas, compared, at common 
velocities, with o,, for He* in helium, and 0,7 
for H° in helium, shows the exchange nature of 
the HetHe interaction compared to that of the 
other isoelectronic structures. 
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DISINTEGRATION OF Te’!*®. C. W. Kocher, 

Allan C. G. Mitchell, C. B. Creager, and T. D. 
Nainan, Physics Department, Indiana University, 
Bloomington, Indiana (Received July 5, 1960). 


The disintegration of Te*® has been studied 
with the help of magnetic spectrometers and scin- 
tillation counters. Two half-lives are observed 
-one of 4.7+ 0.3 days and one of ~16 hours. 
These represent the decay from two isomeric 
states. The isomeric transition is not observed 
because of the strong competition from the decay 
to states of Sb™"°. The 4.7-day disintegration 
occurs by electron capture and exhibits gamma 
rays of the following energies: 0.153, 0.270, 0.930, 
1.10, 1.22, and 2.12 Mev. The following coinci- 
dences are observed: (2.12, 0.270), (1.22, 0.930, 
0.153), (1.10, 0.270). The K/(L+™M) ratio for the 
line at 0.153 Mev is 7.85 and for the line at 0.270 
Mev is 5.77. The 16-hour isomer exhibits posi- 
trons of end-point energy 0.627 Mev; and two 
gamma rays of energy 1.76 and 0.645 Mev. A 
disintegration scheme is given. 


DECAY OF Si**. E. L. Robinson and O. E. John- 
son, Physics Department, Purdue University, 
Lafayette, Indiana (Received April 8, 1960). 


A (2.14 0.3)-sec activity was observed when 
vacuum -distilled magnesium was bombarded with 
8-Mev He® ions. Half-life studies using Nal(T1) 
scintillation counters yielded evidence that this 
activity was due to the decay of a positron- 
emitting isotope with a maximum kinetic energy 
greater than 3.5 Mev. The features of the gamma- 
ray spectrum with the exception of a weak line 
at 824+ 15 kev could be understood in terms of 
the decay characteristics of known radioisotopes. 
An internally consistent argument based on the 
known decay characteristics of reaction products 
that may be expected from energy considerations, 
the results of half-life studies, experimental 
gamma spectra, and nuclear systematics can be 
made to support the conclusion that the (2.1+0.3)- 
sec half-life is that of Si*® produced in the re- 
action Mg“(He®, 7)Si?*, and a consistent decay 
Scheme can be proposed. The ground state of 
Si** (0+) decays by the emission of two positron 
groups to excited states of Al?®. The most in- 
tense transition, E,=3.76 Mev, is to the 0.228- 
Mev state (0+) of Al°®. The second transition, 
E,=2.94 Mev, is to the 1.05-Mev state (1+) of 
Al’, The 1.05-Mev and 0.228-Mev states are 
then connected by an (824+ 15)-kev gamma transi- 


tion. The energies of the positron transitions 
are derived from the known levels of Al** and 
the Si**-Al”® mass difference. 


ANGULAR DISTRIBUTIONS OF THE Be*(d,n)B*° 
NEUTRONS. Richard Bardes and George E. 
Owen, The Johns Hopkins University, Baltimore, 
Maryland (Received July 7, 1960). 


The angular distributions of neutrons to the 
ground state and to the first excited state of B’° 
in the reaction Be®(d,n)B’° have been studied at 
incident deuteron energies of 1.41 Mev, 1.88 
Mev, and 2.35 Mev. A proton recoil spectrom- 
eter utilizing a xenon gas scintillation trigger 
was developed to study this problem. This 
spectrometer operated with a resolution of the 
order of 7% at the neutron energies involved in 
this reaction. 

The analysis of the data in terms of a direct- 
interaction mechanism indicates that the results 
are consistent with an interpretation which in- 
dicates that the majority of the neutrons emitted 
at bombarding energies below the Coulomb bar- 
rier originate in the Be® target nucleus. 


MOTT-SCATTERING ANALYSIS OF LONGITU- 
DINAL POLARIZATION OF ELECTRONS FROM 
Co®. J. S. Greenberg, D. P. Malone, * R. L. 
Gluckstern, and V. W. Hughes, Yale University, 
New Haven, Connecticut (Received July 15, 
1960). 


Mott scattering has been used to analyze the 
degree of longitudinal polarization of beta par- 
ticles emitted from radioactive nuclei. The 
reliability of this method and the influences of 
the various systematic errors associated with 
this method on the accuracy of the measurement 
have been investigated in detail and are discussed. 
On the basis of a linear extrapolation of the 
inverse of the Mott asymmetry to zero scatterer 
thickness, the polarization of 194-kev electrons 
from Co® was found to be -(0.994+ 0.057)v/c 
with all known corrections applied. The effects 
of atomic screening and finite nuclear size have 
not been included. Using the quoted value for 
the polarization measured in the pure Gamow- 
Teller transition in Co yields C,4’=(0.7 to 1.45) 
XC. 

*Cornell Aeronautical Laboratory, Buffalo, New 
York. 
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PHOTONEUTRON CROSS SECTIONS OF COBALT 
AND MANGANESE. P. A. Flournoy, R. S. 
Tickle, and W. D. Whitehead, University of Vir- 
ginia, Charlottesville, Virginia (Received June 
15, 1960). 


The total photoneutron yields of Mn® and Co™® 
have been carefully measured from threshold to 
approximately 30 Mev. Analysis of these data 
using the Leiss- Penfold matrix indicates that 
the cross sections for both elements show a 
splitting in the giant resonance region in accord 
with the predictions of the classical hydrodynamic 
model. The Mn*® peaks occur at energies of 16.8 
+ 0.25 Mev and 19.75+ 0.25 Mev corresponding to 
cross sections of 90 mb and 77 mb, respectively. 
Co**® maxima occur at 16.75+ 0.25 Mev and 18.75 
+ 0.25 Mev with cross sections of 109 mb and 92 
mb. The cross sections o(y,) + o(y, 2) + o(y, np) 
+... integrated to 25 Mev are 627 Mev-mb for 
Mn** and 709 Mev-mb for Co™®. Breit-Wigner 
resonance lines were fitted to both cross sections 
and the intrinsic quadrupole moments determined 
from these fits are +0.78+ 0.10 barn for manga- 
nese and +0.76+ 0.11 barn for cobalt. 


CREATION OF THE HEAVY ELEMENTS BY 
NEUTRON CAPTURE ON A FAST TIME SCALE. 
Peter Fong, Department of Physics, Utica College 
of Syracuse University, Utica, New York, and 
Laboratory for Nuclear Science, Massachusetts 
Institute of Technology, Cambridge, Massachu- 
setts (Received February 10, 1960). 


The mathematical problem concerning the 
creation of the heavy elements by neutron capture 
on a fast time scale involves the solution of more 
than a thousand simultaneous differential equa- 
tions. This set of equations is solved by re- 
placing them with a mathematical equivalent — 

a single partial differential equation which turns 
out to assume the form of the equation of con- 
tinuity. The stream lines of the analogous hydro- 
dynamical problem are determined, from which 
the problem of nucleogenesis may be solved once 
the initial condition is given. By assuming a 
neutron capture mean lifetime of about 0.1 sec 
for all nonmagic-number nuclides and by making 
use of the fact that the neutron magic-number 
nuclides have smaller neutron capture cross 
sections, the position, width, and height of the 
abundance peaks at Br, Xe, and Pt may be de- 
termined; these agree reasonably well with the 
experimental results. 
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INELASTIC DIFFRACTION SCATTERING. 

E. Rost and N. Austern, University of Pittsburgh, 
Pittsburgh, Pennsylvania (Received July 12, 
1960). 


The method of distorted-wave Born approxi- 
mation is shown to be equivalent to the first- 
order adiabatic method for the calculation of 
inelastic scattering. The flexibility of the dis- 
torted-wave method is indicated and tested by 
numerical calculation for a simplified model. 
The calculations treat the inelastic scattering of 
alpha particles from strongly absorbing nuclei, 
and are directly comparable to the adiabatic cal- 
culations considered by Blair in Fraunhofer 
approximation. Good agreement is found at for- 
ward angles. Coulomb wave functions are used, 
and the energy difference between initial and 
final states is taken into account. Comparison is 
made to a few representative experiments and 
good agreement is found. 


GAMMA TRANSITIONS FOLLOWING THE DE- 
CAY OF NICKEL-65. R. Jambunathan, M. R. 
Gunye, and Babulal Saraf, Atomic Energy Estab- 
lishment, Trombay, Bombay, India (Received 
April 25, 1960; revised manuscript received 
August 8, 1960). 


The gamma transitions following the decay of 
2.6-hour Ni®, have been studied by coincidence 
scintillation spectrometry. Gamma rays of en- 
ergies 0.37, 0.51, 0.61, 0.77, 0.85, 0.95, 1.11, 
1.48, 1.62, and 1.72 Mev having intensities of 
55, 1.5, 0.6, ~1, ~0.6, ~0.3, 140, 290, 10, and 
6 per thousand disintegrations, respectively, 
are emitted in this decay. The gamma-gamma 
coincidence study is consistent with energy levels 
of Cu® at 0.77, 1.11, 1.48, 1.62, and 1.72 Mev. 
The angular correlation of 0.37- and 1.11-Mev 
gamma rays suggests the spin assignment of 
5/2 for the 1.11-Mev level and 7/2 or 3/2 for 
the 1.48-Mev level. Both the gamma transitions 
are mixtures of E2 and M1. 


DECAY SCHEME OF Dy’®*. R. C. Greenwood 
and E. Brannen, University of Western Ontario, 
London, Ontario, Canada (Received July 5, 1960). 


The electron capture decay of Dy’® has been 
shown to occur appreciably only to the 58-kev 
state and ground state of Tb'®®, with branching 
ratios of 37% and 63%, respectively. The X in- 
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ternal conversion coefficient of the 58-kev tran- 
sition has been measured to be ax = 8.5*9:3, 

which is consistent with a predominantly M1 tran- 
sition, with the possibility of some E2 mixing. 

The orbital electron capture ratio to the 58-kev 
state of Tb**® was obtained as Py y.../Px =0.17 
+0.15. From this value, a minimum decay energy 
to the ground state of Tb*®® was assigned as 230 
kev. The upper limit on transitions to the 138-kev 
and 364-kev states of Tb’®® was placed at 0.1%. 

A spin and parity of 3/2- was assigned to the 
ground state of Dy’®® from a consideration of 

logft values. An upper limit of 450 kev was placed 
on the decay energy to the ground state of Tb’®’. 


NUCLEAR SPINS OF SIX NEUTRON-DEFICIENT 
GOLD ISOTOPES. W. Bruce Ewbank, Lawrence 
L. Marino,* William A. Nierenberg, Howard A. 
Shugart, and Henry B. Silsbee, f Lawrence Radia- 
tion Laboratory and Department of Physics, 
University of California, Berkeley, California 
(Received July 5, 1960). 


The nuclear spins of six neutron-deficient 
isotopes of gold have been measured by atomic- 
beam magnetic-resonance methods. The results 
of these measurements are as follows: Au’* 
(3.0 hr), 3/2; Au’® (4.8 hr), 1; Au’® (17.5 hr), 
3/2; Au'™ (39 hr), 1; Au’® (190 day), 3/2; Au’®® 
(5.6 day), 2. 


*Present address: Convair, San Diego, California. 
tPresent address: University of Washington, Seattle, 
Washington. 


FINE STRUCTURE IN THE ENERGY SPECTRA 
OF PHOTOPROTONS FROM He*. Carmelo 
Milone, Istituto di Fisica dell’Universita, Centro 
Siciliano di Fisica Nucleare, Catania, Italy 
(Received June 27, 1960). 


The energy spectra of the photoprotons from He* 
irradiated with 31- and 32-Mev bremsstrahlung 
were studied by means of nuclear photoemulsions 
in AE, steps of 0.4/3 Mev. In the region of the 
giant resonance several peaks have been dis- 
tinguished. The more evident of these have 
widths less than 0.5 Mev and are located at He* 
excitation energies of 24.7+ 0.2 and 26.1+0.1 
Mev. Previous experiments on the photoproton 
spectra from He* could not give evidence of 
these narrow resonances because the spectra 
were examined in energy steps much larger than 


the width of these resonances. 

Other types of experiments which have been 
performed in order to obtain information on the 
existence of excited states in He* are discussed. 
The fact that narrow resonances are distinguished 
only in the present experiment is attributable to 
the selective nature of the (y, p) process, the 
most important contribution in the He*(y, p) re- 
action coming only from the states of He* having 
J=1~ and T=1. 


INTERPRETATION OF ISOMERIC CROSS- 
SECTION RATIOS FOR (n, y) AND (y,”) REAC- 
TIONS. J. R. Huizenga and R. Vandenbosch, 
Argonne National Laboratory, Argonne, Illinois 
(Received July 11, 1960). 


The relative probability of forming each mem- 
ber of a pair of nuclear isomeric states has been 
compared with theoretical predictions in order 
to learn which nuclear parameters can be deter- 
mined from these data. For thermal neutron 
capture reactions, the observed ratios do not give 
much information about the dependence of the 
nuclear level density on spin, but they are con- 
sistent with a spin cutoff factor, exp[-(J+4)?/20’], 
where o <5. The calculations are sufficiently 
consistent with experiment to make their pre- 
dictions usable as a guide for assigning spins to 
the compound states formed in thermal or reso- 
nant energy neutron capture. For (y,”) reactions, 
the calculations reproduce the energy dependence 
of the experimentally observed isomeric cross- 
section ratios. In order to obtain quantitative 
information about the spin dependence of the nu- 
clear level density, it is necessary to consider 
reactions where particles are emitted which can 
carry off enough angular momentum to reach 
many spin states of the residual nucleus. 


ISOMERIC CROSS-SECTION RATIOS FOR RE- 
ACTIONS PRODUCING THE ISOMERIC PAIR 
Hg'*’,197"_ R. Vandenbosch and J. R. Huizenga, 
Argonne National Laboratory, Argonne, Illinois 
(Received July 11, 1960). 


Excitation functions and isomeric cross-section 
ratios were determined for the Au'®( p, n) and 
Au’®"(d, 2n) reactions. Isomeric cross-section 
ratios were also determined for the Hg'™(n, y), 
Hg'(d, p), Hg'™(n, 2n), Hg’*(a, an), and 
Pt(a,xm) reactions. Those reactions for which 
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compound nucleus formation predominates are 
treated by a statistical model. In this treatment 
the compound state is characterized by angular 
momentum and excitation energy. The analysis 
yields a value of 4+1 for the parameter o which 
characterizes the dependence of the nuclear level 
density on angular momentum. Relatively small 
amounts of angular momentum are transferred 
in reactions which proceed predominantly by a 
direct-interaction mechanism. Such reactions, 
therefore, give a larger yield of the isomer with 
spin closer to that of the target nucleus. 


Ni*®(n, p)Co® CROSS SECTION FOR NEUTRONS 
OF ENERGIES BETWEEN 2.2 AND 3.6 Mev. 

L. Gonzdlez, J. Rapaport, and J. J. van Loef, 
Instituto de Fisica y Matematicas, Universidad 


de Chile, Santiago, Chile (Received July 1, 1960). 


The relative cross section for Ni™®(n, p)Co™ 
has been measured by an activation method for 
neutron energies of 2.2 to 3.6 Mev. The absolute 
cross section at 3.55-Mev neutron energy is 
found to be 167+ 27 mb by comparing the inten- 
sity of the 810-kev gamma ray of Co™ with the 
amount of Si* formed in the P™"(n, p)Si* reaction, 
which has a previously known cross section of 
96.2+9.0 mb. The cross section rises from 120 
mb at 2.2 Mev to a maximum value of 185 mb at 
3.0 Mev. Observed cross sections are compared 
with the predictions of the statistical theory of 
nuclear reactions by using a square well poten- 
tial and a diffuse edge potential. In neither case 
could a satisfactory agreement with experimental 
data be obtained. 


a+ 
FISSION OF GOLD BY CARBON IONS. Glen E. 
Gordon,* Almon E. Larsh, Torbjgrn Sikkeland, 
and Glenn T. Seaborg, Lawrence Radiation 
Laboratory and Department of Chemistry, Uni- 
versity of California, Berkeley, California (Re- 
ceived July 14, 1960). 


Angular distribution and kinetic-energy spectra 
of fragments, and cross sections for fission of 
gold with 68- to 124-Mev C”™ ions have been ob- 
tained by observation of the fragments in two 
types of detectors, gas scintillation chambers 
and silicon p-m junctions. From the parameters 
used to fit the angular distributions to the theo- 
retical curves of Halpern and Strutinski, we 
have obtained the average excitation energy of 
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the fissioning nucleus at the time of fission. This 
quantity is approximately 25 Mev, which is near- 
ly independent of bombarding energy, suggesting 
that fission is preceded by the emission of sev- 
eral particles from the compound nucleus. The 
fission cross section increases from a value of 
100 mb at 68 Mev to 1.28 b at 124 Mev. Over 
this range of bombarding energies, the total 
fragment kinetic-energy release rises from 142 
+6 to 146+6 Mev. At all bombarding energies, 
the variation of laboratory-system kinetic energy 
of the fragments with laboratory-system angle 
indicates full momentum transfer by the bombard- 
ing particle to the fissioning system. 


*Present address: Department of Chemistry, Massa- 
chusetts Institute of Technology, Cambridge, Massa- 
chusetts. 


SINGLE AND MULTIPLE STRIPPING OF a 
PARTICLES. M. el Nadi,* Sloane Physics La- 
boratory, Yale University, New Haven, Connect- 
icut (Received June 9, 1960). 


Assuming the stripping process to take place 
mainly at the nuclear surface, expressions are 
derived for the differential cross sections for 
the processes in which one, two, or three nu- 
cleons are stripped from an incident a particle 
by the target nucleus. Comparison between the 
derived formulas and experimental data was car- 
ried out for the (He*,) and (a,) reactions and 
shows fairly good agreement. 


*Present address: Faculty of Science, Cairo Uni- 
versity, Cairo, Egypt. 


HYPERFINE STRUCTURE OF AMERICIUM-241. 
Richard Marrus, William A. Nierenberg, and 
Joseph Winocur, Lawrence Radiation Laboratory, 
University of California, Berkeley, California 
(Received July 8, 1960). 


Atomic-beam measurements on Am™? have 
confirmed that the nuclear spin is /=5/2, and 
have established that the electronic ground state 
is characterized by an angular momentum J=7/2 
and a Landé splitting factor g7 =- 1.9371(10). In 
addition, the magnetic-dipole and electric- 
quadrupole hyperfine-structure coupling con- 
stants have been found to be, respectively, A 
=+17.144(8) Mc/sec and B = 123.82(10) Mc/sec. 
It is hypothesized that these values arise from 
a state that is primarily formed from the Hund’s 
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rule term °S of the configuration (5f)"(7s)*. How- 
ever, important contributions to the measured 
gy A, and B values are shown to come from the 
admixture of other terms in the ground state by 
means of the spin-orbit interaction, and also 
from the excitation of s electrons in the core to 
higher s states. . 


RADIOACTIVE DECAY OF Yb'®’. R. G. Wilson 
and M. L. Pool, Department of Physics and 
Astronomy, Ohio State University, Columbus, 
Ohio (Received May 27, 1960). 


Erbium oxide enriched to 35.1% in the mass 
number 164 was irradiated with 17- and 24-Mev 
alpha particles. An activity decaying by electron 
capture with a half-life of (17.7+ 0.2) minutes 
was produced and its assignment to Yb’®’ con- 
firmed by the identification of the thulium K 
x ray and the well-known daughter activity of 
Tm'*’, Gamma rays with energies of 106, 113, 
and 176 kev were observed in the spectrum of 
Yb'*’. Gamma-gamma coincidence measure- 


ments were performed on the observed radiations. 


Transition probabilities, multipole order admix- 
tures, and electron capture branching ratios 
have been calculated. Approximately 89% of the 
electron capture transitions occur to the 292.7- 
kev (7/2-) level of Tm?®”. 


LIFETIMES OF ROTATIONAL LEVELS IN 

SOME ODD-A NUCLEI MEASURED BY A MICRO- 
WAVE METHOD. A. E. Blaugrund, Y. Dar, and 
G. Goldring, Department of Physics, The Weiz- 
mann Institute of Science, Rehovoth, Israel (Re- 
ceived June 1, 1960; revised manuscript received 
August 15, 1960). 


Measurements of mean lives of the first ex- 
cited states of a number of odd-A nuclei have 
been carried out with the intention of obtaining 
more accurate information on magnetic transi- 
tion probabilities than is at present available. 
Most of these mean lives are in the range of 
2-20x10""' sec, and for this measurement a 
special timing apparatus has been devised which 
consists of a microwave beam pulsing system 
and a beta-ray spectrometer and cavity combina- 
tion which acts as a microwave shutter. In four 
cases the exponential decay could be determined 
directly; these were the first excited states of 
Lu'’, Hf!’®, and Ta‘! and the 118-kev level of 


Tm’®. The results for these levels are: 7(Tm*®*) 
= (9.0+ 0.4) x10™™ sec, 7(Lu’”’) = (14.64 1) x10-™ 
sec, T(Hf!”*) = (5.44 0.4) x10"™ sec, 7(Ta‘®) = (6.0 
+0.5)x10°"' sec. Lifetime measurements were 
also carried out for the first excited states of Hf'”, 
Re'®>, and Re*®’. The limitations and scope of this 
method of measuring lifetimes are discussed. 


SURVEY OF TRITIUM-PRODUCING NUCLEAR 
REACTIONS. Jose Gonzalez-Vidal and William H. 
Wade,* Lawrence Radiation Laboratory, Univer- 
sity of California, Berkeley, California (Received 
July 8, 1960). 


(p,t), (d,t) and (a, t) reactions have been in- 
vestigated throughout the periodic table by bom- 
barding stacked metal foils and determining 
directly the tritium produced in the reaction. In 
the (a, ¢) reactions, there is conclusive evidence 
that most of the tritons are produced with high 
energies, thus indicating the presence of direct- 
interaction processes. The curve representing 
the integrated cross section vs Z of the target 
rises with decreasing Z. This, and the appear - 
ance of low-energy peaks in the individual ex- 
citation functions of low-Z targets indicate that 
at low and intermediate values of Z the relative 
number of low-energy tritons increases. These 
tritons are probably the product of a compound- 
nucleus mechanism. For the (p, ¢) and (d, t) re- 
actions the same compound-nucleus and direct- 
interaction effects are noticed. 


*Chemistry Department, University of Texas, Austin, 
Texas. 


ENERGY LEVELS OF Ne” AND Ne*. Joan M. 
Freeman,* Laboratory for Nuclear Science, 
Massachusetts Institute of Technology, Cam- 
bridge, Massachusetts (Received June 28, 1960). 


The energies of twenty-seven excited states of 
Ne” and seventeen of Ne”* have been established 
by a study of the (d, p) reactions in Ne”° and in 
Ne**. Seven deuteron bombarding energies be- 
tween 4.75 and 7.5 Mev were used, and the target 
was a cell with thin windows, containing either 
natural or Ne**-enriched neon gas. The proton 
groups emitted at 90 degrees were analyzed with 
a broad-range magnetic spectrograph. The es- 
tablished calibration of the spectrograph, in 
terms of alpha particles from polonium, together 
with some well-known reactions, provided the 
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data necessary for precise energy determina- 
tions. The Q values obtained for the ground- 
state reactions Ne*°(d, p)Ne™ and Ne”*(d, p)Ne”* 
were 4.534+ 0.009 Mev and 2.791+0.009 Mev, 
respectively. Deuteron groups inelastically 
scattered from Ne”° and Ne” were also observed. 
The higher energy groups were used in the cali- 
bration procedure; the lower energy groups gave 
the results 4.250+ 0.008 Mev and 3.356 + 0.008 
Mev for the second excited states of Ne”° and 
Ne”, respectively. The interpretation of the 
level schemes of Ne” and Ne” in terms of the 
collective model has been considered and is dis- 
cussed here briefly. 


*Now returned to Atomic Energy Research Estab- 
lishment, Harwell, Berkshire, England. 


DECAY OF 16-MINUTE Ta’®?”. A. W. Sunyar, 
Brookhaven National Laboratory, Upton, New 

York, and P. Axel, Physics Department, Uni- 
versity of Illinois, Urbana, Illinois (Received 

July 5, 1960). 


A revised decay scheme of the 16-minute iso- 
mer Ta’*?” involves three states of Ta’®? at 
excitation energies of 147 kev, 319 kev, and 503 
kev. Transition multipolarities were classified 
by measuring both gamma-ray coincidences and 
internal conversion electrons. The 503-kev iso- 
meric state decays mainly (98%) by a 184-kev 
E3 transition to the 319-kev state. This 319-kev 
state decays mainly (94%) to the 147-kev state 
by means of a 172-kev transition that is pre- 
dominantly M1. The 147-kev transition to the 
ground state is also predominantly M1. Two of 
the three possible crossover transitions were 
observed; a 356-kev M4 transition originates at 
the isomeric level and a 319-kev E2 transition 
connects the second excited state with the ground 
state. 

The 147-kev and the 319-kev states are probably 
the first and second excited rotational states of 
the ground-state configuration. The relevant ro- 
tational parameters are of particular interest 
because very little is known about moments of 
inertia for odd-odd nuclei. 

The relative probabilities were determined for 
pile neutron activation of the 16-minute isomer, 
Ta'®?™ and the 112-day ground state, Ta’®’. If 
the ground-state formation cross section is taken 
as 21 barns, the corresponding value for the iso- 
mer is only 9 mb. 


294 


LEVEL STRUCTURE OF Eu’. R. E. Sund and 
M. L. Wiedenbeck, Harrison M. Randall Labo- 
ratory of Physics, The University of Michigan, 
Ann Arbor, Michigan (Received July 1, 1960). 


Gamma rays in Eu’ following the decay of 
Sm*** and Gd'* have been studied using coinci- 
dence and angular correlation methods. Results 
for the strong transitions are in agreement with 
the decay scheme given by McCutchen. Measure- 
ments were made of spectra of gamma rays in 
coincidence with the x-ray, 70-kev, and the 97- 
kev and 103-kev transitions in the decay of Gd, 
and with eight energy regions in the decay of 
Sm***. A number of new, weak transitions were 
observed in the decay of Sm'®’, and a consistent 
decay scheme is proposed. Directional correla- 
tion measurements were made on the 70-kev— 
103-kev cascade from the decay of Sm**® and 
from the decay of Gd'**. Possible spin assign- 
ments are discussed. 


PSEUDOSCALAR INTERACTION IN NUCLEAR 
BETA DECAY. C. P. Bhalla,* University of 
Tennessee, Knoxville, Tennessee, and M. E. 
Rose, Oak Ridge National Laboratory, Oak Ridge, 
Tennessee (Received July 5, 1960). 


The experiments on the allowed beta transitions, 
which lead almost uniquely to the V - 1.2A inter- 
action, do not have any bearing on a possible con- 
tribution from the pseudoscalar interaction. To 
determine whether or not any contribution from 
the pseudoscalar interaction is really needed, an 
examination has been made of the £ longitudinal 
polarization and the 8 shape factor in the 0+0 
(yes) beta transitions. The theoretical polari- 
zation for the mixture of the pseudoscalar and the 
axial vector interactions has been developed. In 
this work, the formulation of the pseudoscalar 
interaction as given by Rose and Osborn has been 
used. The numerical results on the 8 longitudinal 
polarization and the shape factor depend on two 
parameters, namely, the coupling constant ratio, 
Cp/MC 4, and x, the ratio of the two relevant 
nuclear matrix elements. M is the nucleon mass 
in units of the electron mass. The electronic 
functions occuring in the theoretical formulas 
for these effects are tabulated for Pr'** (0° - 0°) 
and Ho’®* (0° + 0*), All the electronic radial func- 
tions were computed considering the nucleus as 
a sphere of a uniform charge distribution with a 
nuclear radius of 1.2A¥“°x10"'% cm, and taking 
into account the finite deBroglie wavelength effect. 
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The results of extensive numerical analysis are 
presented. We conclude that the absence of the 
pseudoscalar interaction is consistent with the 
existing experimental data. The value of Cp/MC A? 
which also gives a satisfactory fit to the experi- 
mental data, depends on A. The upper limit of 

the value of |Cp/MC,| is found to be 0.05 for 
\~{=200. In this work, time-reversal invariance 
is assumed valid for the weak as well as the 

strong interactions, and the two-component theory 
of the neutrino has been used. 


*Present address: Westinghouse Electric Corporation, 
Atomic Power Department, Pittsburgh, Pennsylvania. 


DEPOLARIZATION DUE TO TARGET MOTION. 
Albert E. Werbrouck,* Palmer Physical Labora- 
tory, Princeton University, Princeton, New 
Jersey (Received July 14, 1960). 


The projection upon both of the lab scattering 
planes of a spin component along the center -of- 
mass angular momentum direction is computed 
for arbitrary target motion. The resulting de- 
polarization of A hyperons produced by 1080- 
Mev/c pions on nuclear matter and on deuterons 
is presented as a function of production angle. 
At this energy, the polarization, if constant for 
all production angles, would be reduced to 85% 
and 95% of the original value in heavy nuclei 
with respect to the 7x A and K XA planes, respec- 
tively. In deuterium, the reductions would be to 
98.4% and 99.5% of the initial polarization. 


*Present address: Istituto di Fisica, Turin, Italy. 


SLOW-NEUTRON INTENSITY AT HIGH BAL- 
LOON ALTITUDES. R. C. Haymes and S. A. 
Korff, New York University, University Heights, 
New York, New York (Received June 22, 1960). 


In order to extend measurements of cosmic-ray 
neutrons to top balloon altitudes, a flight was 
conducted on August 24, 1959, in which our sys- 
tem of identical counters filled with enriched and 
natural B*°F, was used. An altitude of 122 000 
feet at a geomagnetic latitude of 55°N was at- 
tained, corresponding to 4-millibar pressure, 
which is considerably higher than previous meas- 
wements. The counting rate dropped to about 1% 
of the maximum rate, which occurs at 60000 
feet. If one extrapolates this curve to the top of 
the atmosphere, one obtains an upper bound for 
the slow-neutron albedo of 0.03 neutron per 


cm’-sec. The background due to highly ionizing 
events in the counters increased exponentially 
to about 63 000 feet, with a derived mean free 
path of 152 g/cm*. The background continued 
an exponential increase up to 122000 feet, but 
with a mean free path of 15 g/cm’. 


METHOD FOR DETERMINING THE K° SPIN. 
Philippe Eberhard, Laboratoire de Physique 
Atomique et Moleculaire, College de France, 
Paris, and M. L. Good, University of Wisconsin, 
Madison, Wisconsin (Received July 11, 1960). 


A method is presented for determining the K°® 
spin, using only angular momentum conservation 
and the rules of quantum mechanics. The method 
is based on a proof that, in the reaction 71> +p—> 
A+K°®, the K,° decay intensity associated with a 
given direction of the A spin cannot be isotropic 
at any production angle, if the K° spin is equal 
to or greater than two. 

Account is taken of the particle-mixture aspect 
of the K,° and of possible magnetic moment pre- 
cession. 


DECAY OF BOUND MUONS. Victor Gilinsky 
and Jon Mathews, California Institute of Tech- 
nology, Pasadena, California (Received July 
11, 1960). 


The decay rate of bound muons from the K shell 
is calculated for elements up to lead. The muon 
is represented by a relativistic wave function for 
a point nucleus with parameters adjusted for 
finite-nuclear-size effects. The outcoming elec- 
tron is represented by a Sommerfeld-Maue wave 
function. The decay rate is found to be a mono- 
tonically decreasing function of the atomic num- 
ber. An estimate of the error is obtained by 
comparing with an exact calculation for the low- 
est electron angular state. The spectrum and 
angular distribution of the electrons are also 
presented. 


LOW-ENERGY NEUTRAL PION PHOTOPRODUC- 
TION IN H, AND He. A. Odian,* G. Stoppini,t and 
T. Yamagata,! University of Illinois, Urbana, 
Illinois (Received July 11, 1960). 


The results of an experiment on the photopro- 
duction of neutral pions from H, and He in the 
very low energy region are compared with the 






295 








VoLUME 5, NUMBER 6 


PHYSICAL REVIEW LETTERS 


SEPTEMBER 15, 1960 





dispersion relation calculations of Chew, Gold- 
berger, Low, and Nambu (C.G.L.N.). For H,, 
agreement was found for incident y-ray energies 
above 155 Mev for N*=+0.065 and disagreement 
below for all constant values of N*. The He re- 
sults give support to the C.G.L.N. effective- 
range formulas for the small p-wave phase shifts. 


*Now at the Istituto Superiore di Sanita, Rome, Italy. 

+Now at the Istituto di Fisica dell’Universita, Rome, 
Italy. 

tNow at CERN, Geneva, Switzerland. 


LEPTON DECAYS OF HYPERONS. D. R. Har- 
rington, Carnegie Institute of Technology, Pitts- 
burgh, Pennsylvania (Received July 5, 1960). 


An expression is given for the lepton energy 
spectrum and asymmetry for the lepton decay 
modes of hyperons in terms of six form factors. 
A weak interaction having the form of the usual 
lepton current coupled to an unspecified strong 
current is assumed, and the momentum transfer 
dependence of the form factors is neglected. 
Some applications to the decays A~p+e+v and 
A~p+u+v are given. 


COMPLEX SINGULARITIES OF PARTIAL-WAVE 
AMPLITUDES IN PERTURBATION THEORY. 

J. G. Taylor and A. E. A. Warburton, Depart- 
ment of Applied Mathematics and Theoretical 
Physics, Cambridge University, Cambridge, 
England (Received July 11, 1960). 


We show that the complex singularities which 
invalidate Mandelstam’s representation do not 
cause complex singularities of partial wave 
amplitudes for two-particle scattering processes, 
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except for the expected “kinematic” complex 
branch points. 


OSCILLATORY CHARACTER OF REISSNER- 
NORDSTROM METRIC FOR AN IDEAL CHARGED 
WORMHOLE. John C. Graves and Dieter R. 
Brill, Palmer Physical Laboratory, Princeton, 
New Jersey (Received July 15, 1960). 


A transformation is presented to remove co- 
ordinate (“pseudo”) singularities from metrics 
of a certain class, a special case of which is the 
transformation of Kruskal, extending the Schwarzs- 
child metric beyond its pseudosingularity. The 
transformation is applied to the Reissner-Nord- 
strom metric, which describes a concentration 
of charge and mass in general relativity. On an 
initial surface this metric shows the same gen- 
eral behavior as the Schwarzschild metric, de- 
scribing a “wormhole,” or bridge, between two 
asymptotically flat spaces, but with electric flux 
flowing through the wormhole. It is found that 
the region of minimum radius, the so-called 
“throat” of the wormhole, begins to contract, 
but reaches a minimum and re-expands after a 
finite proper time, rather than pinching off as 
in the Schwarzschild- Kruskal case: the radius of 
the throat pulsates periodically in time, “cush- 
ioned” by Maxwell pressure of the electric field 
through the throat. The motion of charged par- 
ticles in this metric is investigated, and it is 
shown that no particle can hit the geometric 
singularity at y=0, (1) quite in general, provided 
only that the mass of the test particle exceeds 
the value associated in general relativity with 
its charge, and (2) in particular when the test 
particle has no charge at all; but (3) such col- 
lisions are not avoided when the throat itself is 
not endowed with an electric flux. 





